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PREFACE 


Properties of dielectric materials are of interest to scientists 
in various branches: physicists, chemists, electrical engineers, 
biologists. Their interests concern different aspects; thus an 
electrical engineer will require the dependence of dielectric loss 
on frequency and temperature in order to find a substance 
which is nearly loss free in a certain range. The chemist on the 
other hand can use this knowledge to draw conclusions on the 
properties of molecules. For these, and for many other pur- 
poses, it is imperative to have a theory of dielectrics. 

The present book is intended to give a systematic account of 
the theory of the dielectric constant and of dielectric loss. This 
it is hoj)ed will satisfy the requirements of the various branches 
of research interested in dielectrics. In writing this account I 
found that the subject deserves interest also from a methodical 
point of view as an ap])lication of classical statistical mechanics. 
That this application is far from trivial is shown by some 
of the controversies in the literature which have lasted until 
very recent years ; also the general theorems derived in § 7 seem 
to be novel. 

It was my intention to write this monograph for the use of 
applied scientists. I hope, however, that the sections dealing 
with the general theory will also be of value to students. The 
required mathematical technique only occasionally exceeds 
acquaintance with calculus ; even so I have been told that its 
extensive use might be too heavy for biologists. The reader is 
assumed to have a certain elementary knowledge of atomic and 
molecular physics, statistical mechanics, and electrostatics. 
Quantum mechanics will not be required; its relation to the 
theory of dielectrics is discussed in van Vleck’s book {VI), 

Units unless stated otherwise refer to the electrostatic c.g.s. 
system. Vectors are represented by bold type. Unfortunately 
it was not always possible to avoid repetition of symbols. The 
meaning of the symbols cx:, ~ is respectively propor- 
tional to, approximately equal to, order of magnitude of. As 
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usual, k is Boltzmann’s constant, h is Planck’s constant, 
% = hj27T. 

I am grateful to a great number of colleagues for help of 
various kinds. To T)r. F. C. Frank, Professor Willis Jackson, 
Dr. H. Pelzer for reading all, or parts of, the manuscript, and 
making valuable suggestions; to Dr. R. Sack for his help in 
reading the proofs ; to Mr. S. Zienau for making the index. My 
particular thanks, however, are due to Dr. B. Szigeti who helped 
to collect the experimental material, and to Dr. J. H. Simpson 
who read the first draft of the manuscript and made many useful 
suggestions. 

I should like to use this opportunity to express my ap])recia- 
tion to the British Electrical and Allied Industries Research 
Association (E.R.A.) without whose support much of the work 
described here would have remained undone. 
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CHAPTER I 


MACROSCOPIC THEORY 

1. Static fields 

The electric properties of dielectric substances are usually 
described in terms of the dielectric constant. For most materials 
this quantity is independent of the strength of the electric field 
over a wide range of the latter, but in the case of alternating 
fields depends on the frequency. It also depends on parameters, 
such as the temperature, which define the state of the material. 
In the macroscopic (phenomenological) theory, which will be 
summarized in the present chapter, the dielectric constant is 
supposed to be known empirically. The purpose of the rest of 
the book will be to derive the dielectric constant (and its varia- 
tion with tem])erature, frequency, etc.) from the atomic structure 
of the material. 

Throughout this book we shall be interested in homo- 
geneous substances only, and in electric fields which are inde- 
pendent of the space coordinates, although they may dej)end 
on time. 

Consider now a condenser consisting of two parallel plates in 
vacuum whose distance apart d is small compared with their 
linear dimensions, and suppose that they are charged electri- 
cally with charges -f cr^ and —gA respectively. A is the surface 
area of a plate, and hence +cr is the charge per unit area which 
will be denoted as surface density. The charges give rise to an 
electric field which inside the condenser is practically homo- 
geneous and directed perpendicular to the surface. Its amount is 

given by (vacuum). (1.1) 

In this equation the factor 47 t is due to the particular way in 
which the unit of electric charge is defined, leading to the electro- 
static e.g.s. system of units usually chosen in atomic physics. 
Equation 1 . 1 should be completed by the definition of the electric 
field as the force (in dynes) acting on the unit of electric charge 
which then leads to a unique definition of this unit. 
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Between the two plates of the condenser a voltage can be 
measured whose absolute value 0 is given by 

® = \E\d. (1.2) 

Now suppose the space between the plates to be filled with a 
homogeneous dielectric material while the charge on the plates 
remains unaltered. This will cause the voltage to drop to a 
smaller value, and the ratio of its former value to the present 
one is denoted by the static dielectric constant Since equation 
1.2 still holds, the electric field strength has also decreased, its 
present value being ^ ^ j 3^ 

From equations 1.1 and 1.3 it follows that the drop in the field 
strength accompanying the insertion of the dielectric might also 
be achieved by reducing the surface-density a of the electric 
charge by the amount 

P = (1.4) 


Therefore the influence of the electric field on the dielectric 
is equivalent to charging the two surfaces of the dielectric with 
charges of opposite sign in such a way that the positive condenser 
plate is faced by the negatively charged surface of the dielectric, 
and vice versa. The surface density of the charge is constant, 
and amounts to P. This behaviour of the dielectric is to be 
expected from the atomic point of view according to which any 
substance containing no net charge consists of an equal number 
of positive and negative elementary charges. In a dielectric in 
particular these charges cannot move freely through the medium 
(as in a conductor), but they can be displaced. Clearly negative 
charges will be displaced towards the positive plate and con- 
versely. The total charge passing through any unit of area 
within the dielectric, which is parallel to the condenser plates, 
is the same and its amount is equal to P. P, therefore, is called 
the polarization of the dielectric. 

In macroscopic physics the introduction of the polarization 
P through the displacement of charges is rather fictitious 
because these charges cannot be removed from the dielectric. 
They compensate charges of opposite sign but equal absolute 
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value of the condenser plates. These latter charges are some- 
times denoted as ‘bound’ charges, whereas the original charges 
at the condenser plates are called ‘true’ charges (cf. Fig. 1). It 
is customary then to introduce a new field-quantity D, which is 
described in terms of the true charge. This is the electric displace- 
ment, defined by 


In vacuum, therefore, D 
1.3 and 1.5 


Ittct. (1-^) 

” E. In dielectrics, however, using 

D^e,E, ( 1 . 6 ) 

Clearly from equations 1.4, 1.5, and 1.6 

(1.7) 



Fig. 1. Dielectric material with dielectric constant Cg — 2 between condenser 
plates, (a) Macroscopic description. The left-hand side of the figure shows 
true charges 4- > ~ as sources of D. The right-hand side shows that the true 
charges can be considered as composed of bound charges 0, 0 and of free 

charges 1 4- L I being sources of 4ttF and E respectively, and loading to 

D = EF^rrP. 

(b) Atomic description. Only true charges exist, and the field is described 
by E only; the polarization of the dielectric (indicated by -\ h ••• — )» how- 

ever, leads to surface charges wliich compensate some of the charges on the 

condenser plates. 


Thus in macroscopic physics the electric field in a dielectric 
must be described by two field quantities. The electric field- 
strength E and the electric displacement D are usually chosen, 
and the polarization P can then be derived with the help of 
equation 1.7. D is defined by the (true) charge according to 
equation 1.5, and E can be derived from D with the help of a 
special relation 1.6 which is characteristic for the particular 
dielectric material. 

To link up this description with atomic physics it should be 
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noted that the surface charges :^PA give rise to an electric dipole 
moment M of the dielectric given by 

M - PAd - PF, (1.8) 

where V — Ad is its volume. On the other hand, as shown in § 4, 
this electric moment M can be calculated from the configuration 
of the positive and negative elementary charges which form the 
constituents of the substance. Thus from 1.8 with the use of 
equations 1.6 and 1.7 it follows that 

= 47tMIVE = 4.7TPIE. (1.9) 

This equation provides the required link between macroscopic 
and atomic theory. 


2. Time -dependent fields 

Now consider that the charges on the condenser plates, and 
hence the electric field, depend on time. As in the static case, 
a dielectric which is placed between the plates will be polarized 
by the field. The displacement of charges connected with this 
polarization usually shows some inertia. Thus if a constant 
field is suddenly applied the polarization will not reach its static 
value immediately, but will approach it gradually (cf. Fig. 2). 

As in the static case, two field quantities are required to de- 
scribe the electric field inside the dielectric and one usually 
chooses the electric field-strength E and the electric displace- 
ment D. The latter is still defined by equation 1.5, and the 
correlation between E, D, and P given by equation 1.7 still 
holds. Equation 1.6, however, correlating E and TJ, is no longer 
valid in the present case, but has to be replaced by a more general 
relation. 

Consider first the important case of a periodic field, e.g. 

E = Eq cos cat, (2.1) 

where Eq is independent of time and a>/27T is the frequency in 
cycles per second. If a field of this type has persisted for a 
sufficient length of time, D too must be periodic in time. In 
general, however, D will not necessarily be in phase with E, but 
will show a phase-shift <f>, i.e. 

D = Z>ocos(a>f— = DiQOscot+D2Bmu)t, 


(2.2) 
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where according to elementary trigonometry 

™ DqCos^, ~ Do sin <^. (2.3) 

For most dielectrics IJq will be proportional to Dq, but the 
ratio Dq/Eq usually depends on frequency. Therefore, two 



Fig. 2. Time-dei)ondence of the polarization P of a dielectric when a constant 
electric field E is suddenly applied to it. 

different dielectric constants, e^iw) and €.^{co), both frequency 
dependent, can be introduced by 

Dj ^ €^Eq and D^ e^Do. (2.4) 

Thus from 2.3 and 2.4 

tan^ == (2.5) 

^1 

It will be shown in § 3 that is proportional to the energy loss 
in dielectrics. 

As the frequency approaches zero the present description 
must become identical with that given in § 1 . Thus (assuming 
that there exists no dielectric loss in static fields) 

€ 2 (^ 0 ) -> 0 , 6 j( 6 t>) — > €gj as oj —> 0 . ( 2 . 6 ) 

We shall add the further relation (of. also § 10 ) 

ei(ca) 600 , as o) -> 00 , (2.7) 

which is to be understood in such a way that is the value which 
€ 1 ( 0 )) approaches at the highest frequencies contemplated in the 
present book. They correspond to wave-lengths in the infra-red 
region. 
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The above equations can be written in a condensed form by 
introducing a complex dielectric constant 

€ = (2.8) 

and replacing equation 2.1 by 

E (2.9) 

considering, however, only the real part of this equation (which 

is identical with 2.1). Then the real part of the equation 

D eE (2.10) 

is identical with equations 2.2 and 2.4. 

The two dielectric constants and if considered as functions 
of the frequency co are not entirely independent if the relation- 
ship between E and D is a linear one.f This linear relation- 
ship is usually expressed by the principle of superposition and is 
best explained with the help of a more general time dependence 
of the field than has been considered above. Assume that during 
the time interval between axidu-\~du an electric field of strength 
E{u) has been applied to the dielectric and that the electric 
field vanishes outside this time interval. A displacement D will 
result which in view of the inertia of the polarization P will 
persist at times t > u-\-du, but which will gradually vanish. J 
D is thus a function oit—u, i.e. 

D{t~-u) = E{u)oL{t-~u)du if t > u-]~du, 

where (x(t—u) is the decay function describing the gradual 
decrease of P, in paiticular 

(x(t—u)-^0 if ^->00. (2.11) 

The displacement D contains a part which can follow the field 
practically immediately, and which in view of the meaning 
given to will be assumed to be equal to E(u), 

Thus 

D{t--u) == E{u)+ E(u)(x(0) du if u <t < u+du^ 

t Cf. B. Gross [G4] and S. Whitehead [TF5], where further references are 
given. 

X Note that according to 1.7 = 47rP if = 0. 
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where a may be considered to remain at the value a(0) during 
the short interval du. 

Suppose now that at a later time interval between u' and 
u' -\-du^ another field E{u') is applied. Then by the principle of 
superposition it will be assumed that the corresponding dis- 
placement is superposed linearly on the former one. 

This principle of superposition if applied to a continuous time- 
dependent field E{u) initiated at the time u requires that 
the displacement D{t) at the time t is given by 

D{t) = J E{u)oc{t-u) du. (2.12) 

0 

This equation will now be apphed to periodic fields. Intro- 
ducing E from equation 2.1 into 2.12 thus leads to 

t 

D(^) — €ooi?ocosa>^ — Eq j oc{t~u)GO&coudu 

0 

t 

— Eq J a(a;)cosaj(^— x) dx 
0 

ifx = t—u is introduced. It should be noted that in all integra- 
tions t is to be considered as a parameter. Again it will be 
assumed that the field has persisted sufficiently long to make D 
a periodic function of time. This means that t is larger than the 
time at which oL{i) practically vanishes. Then in view of 2.11 
the above integration over the coordinate x can be extended to 
infinity without appreciably altering the value of the integral, i.e. 

00 

D{t)—e„EQCOsiot — Eq j (x{x)coao){t—x) dx, 

0 

or applying a simple trigonometric formula, 

00 

D{t)—€a,EQCosajt = E Q cos cot j oL{x)coaciJx dx + 

0 

00 

-^-^?Q8inco^ J oL{x)^in(x}xdx, (2.13) 
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Comparing this expression with equations 2.2 and 2.4 it follows 
that 

CX3 

€i(a;) — e^o = / a(a:)cosa>:r dx (2.14) 

0 
00 

e 2 (a>) = J a{x)^m(x)X dx. (2.15) 

0 

Both functions ei(w) — foo and e^ico) can thus be derived from the 
same function a:(a:) and, therefore, cannot be independent. 

A calculation carried out in the appendix (A l.iii) shows that 


00 



€i(co) — “ J 
0 


(2.16) 

and 

€2(40) = 1 J 

00 

f {ei(M)— foo} 2 2 

a> — iJL 

(2.17) 


0 

where ja is a variable of integration. Both integrals are principal 
values. 

Equation 2.16 can be used to calculate the static dielectric 
constant from e 2 (^)- Ir ^his case clearly 

00 

= «i(0) = fco + - f Hil^) (2.18) 

TT J fl 

0 

which indicates that substances for which is very small 

cannot show appreciable dielectric losses (which are proportional 
to 62 )- 

Finally, it should be mentioned that the correlation between 
macroscopic and atomic theory is provided by an equation 
similar in nature to equation 1.9. It is based on the fact that 
equation 1.7 holds in the time-dependent case as well as in the 
static one. With a similar argument to that in § 1 it follows that 
the polarization P is equal to the electric moment M per unit 
volume, leading to equation 1 . 8 . Introducing the complex 
dielectric constant c with the help of 2.10, using 1.7, it follows 

(e-l)E = iirMIV, (2.19) 
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in which only the real part of the expression on the left-hand side 
is considered. Thus if 

M — (2.20) 

then €i— 1 4:TrMJVEQ, €2 irrMJVEQ. (2.21) 

3. Energy and Entropy 

A. Static fields'f 

In many text -books the expression 



is stated to represent the electric energy per unit volume of a 
dielectric material with a static dielectric constant in the 
presence of an electric field E, This statement is misleading 
whenever depends on temperature. In fact it suggests that 
the energy difference per unit volume of the dielectric, first in 
the presence and then in the absence of an electric field, is always 
given by the above expression. This energy difference ought, 
however, to depend on the state in which the dielectric is kept 
while the electric field is applied; this might, for instance, be 
done isothermally or adiabatically. A more accurate discussion 
given below shows that the above expression actually is the 
change of the free energy of the dielectric. 

Before this discussion is commenced the reader will briefly 
bo reminded of the two fundamental laws of thermodynamics. 
Consider as a simple example a gas of volume v, pressure p, 
and temperature T. By expanding it by the small volume dv, 
work amounting to pdv will be done. Conservation of energy, 
therefore, requires that 

dU - dQ-pdv (3.1) 

is the change in the energy content of the gas if dQ is the influx 
of heat during the expansion. Equation 3.1 represents the first 
law of thermodynamics in this simple case. An analysis shows 
that the quantity dQ is not a total differential, i.e. that no 
unique function Q of the variables exists such that dQ is the 

t Cf. Abraham-Becker [Alt Chapter XI]. 



MACROSCOPIC THEORY 


10 


I, §3 


difference between two neighbouring values and 
reversible process, the expression 



For a 
(3.2) 


does, however, represent a total differential of a function S, 
the entropy. S is of fundamental importance in connexion with 
the second law of thermodynamics (cf. text-books). 

With the help of S the Helmholtz free energy F can be derived 

F=--U~TS. (3.3) 

It represents the maximum amount of work which the system 
can be made to do in an isothermal (i.e. constant temperature) 
process. 

In the case of dielectric material in an electric field, electro- 
magnetic theory (cf. Appendix A l.i) shows that the quantity 

~EdD (3.4) 

477 


represents the influx of energy into the dielectric (per unit 
volume) if the displacement D is increased by the small amount 
dD. 

Assume now that the volume of the dielectric is always kept 
constant, and that the temperature T is the only parameter 
besides the electric field E considered to be varied. Then the 
increase dU oi the energy U per unit volume of the dielectric 
in a process in which either T or E, or both, are varied slightly 
is given by ^ 

dU ==dQ+~dD (3.5) 


if dQ is the influx of heat per unit volume. 

This equation is similar in structure to equation 3.1 for a gas 
if E and D are replaced by —p and v. The relation existing in 
gases between p, v, and T (the equation of state) is, however, 
different from the relation between E, D, and T. For the latter, 
equation 1.6 will be supposed to hold with a dielectric constant 
€g which may depend on T but is independent of E, Thus 

dD = d{t,E) = €,dE-^Ede, = €,dE-\-E^^dT, 
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which means that a variation of D may be due to a change in 
the field-strength E at constant temperature, and to a change 
in temperature at constant E. For the following it will be useful 
to take T and E^ as the independent variables. Equation 3.5 
representing the first law of thermodynamics then becomes 






A further relation will now be obtained from the entropy law 
according to which dS^ given by equation 3.2, is a total differen- 
tial. This means that a unique function S{T,E‘^) must exist 
such that 


Thus if it is found that 

dS - A(T, E^) dT+B(T, E^) d(E^), (3.8) 

where A and B are both functions of the two variables T and 
the condition that dS is a total differential requires that 

“ = ^ (3.9) 

ST 8(B*) ' ' 

because both sides of the equation 3.9 are equal to d'^SjdTd[E^), 
Now according to 3.2, inserting dQ from 3.6, 




^ndTj 




d{E^). (3.10) 


T his equation is of the type 3.8, and equation 3.9, therefore, 
becomes 

8 ji( dU e,\i _ a jiidU 
8T[T\d(E^) »7rji d(E^}\T\dT 47r 8T)j' 

Carrying out the differentiations one finds 

^ r 'T 

8(E^) ~ 8irl * dTj 


Integrating with respect to E’^ yields the energy density 

U-F, (!■) + (..+ (! 


(3.11) 
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where Uq(T) is independent of but depends on T and thus 
represents the energy of the dielectric in the absence of a field. 

The entropy S can now easily be calculated because on com- 
paring equation 3.10 with 3.7, both dSjdT and dSld(E^) are 
known if U is introduced from equation 3.11. Thus 

dS ^ 1 dU.E^d^e, dS __ 1 d€, 
dT ““ T "^ 877 d{E^) 877 dT' 

or integrating 8 = (3.12) 

where ^o(^) entropy in the absence of a field. From 

equation 3.3 one finally finds for the free energy 

F^F^{T) + ^-^, (3.13) 

077 

where Fq{T) is the free energy in the absence of a field. This 
proves our original contention. 

The above expressions for f/, S, and F are very instructive. 
Thus equation 3.13 for the free energy shows (in analogy to the 
meaning of F in gases) that the amount of electric energy avail- 
able in an isothermal reversible process is E^lSn. 

From equation 3.11 for the energy it follows that for 
substances such as dilute dipolar gases for which (cq is inde- 

pendent of T) ^ .„+constant/r, (3.14) 

the energy change due to the field is given by E^/Stt. Thus in 
this case the temperature-dependent part of does not make 
any contribution to the energy. The remaining free energy 
{€g—€Q)E^/H7T available besides CqE^IStt is thus entirely due to a 
change in entropy, t 

Finally, equation 3.12 shows that the entropy is increased by 
the field if deJdT is positive, and decreased if this quantity is 
negative. Since the entropy is a measure of the molecular dis- 
order, an external field creates order in dipolar liquids and gases 
for which decreases with increasing T, This may be expected 

t Near the absolute zero of temperature this would contradict the third 
law of thermodynamics ; a temperature dependence of the type 3.14 can, 
therefore, not be valid near T = 0, 
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because the field will orientate some of the dipoles which in the 
absence of a field are at random. In some dipolar solids, on the 
other hand, increases with T, which means that an external 
field increases the disorder. This too is understandable if one 
assumes that in the absence of a field the dipoles are in a well- 
ordered state as may be expected in solids. The field by turning 
some of the dipoles into different directions can thus only 
decrease the existing order. 



Fig. 3. Schematic temperature -dependence of the dielectric constant and 
of the entropy change S oc deJdT due to the polarization by a field. If /S' > 0 the 
field creates disorder, if < 0 it creates order. Near the absolute zero of 
temperature the substance is already perfectly ordered. Hence deJdT cannot 
bo negative near T = 0. 


B. Periodic fields 

The corresponding calculations in the general case of a time- 
dependent field become very involved. It is easy, however, to 
consider the isothermal case in the presence of a periodic field, 
and to calculate on an average over one period the amount of 
electric energy which is transformed into heat. 

On an average over a period the energy U of the dielectric 
cannot change because the temperature is kept constant, and 
the field E is periodic. Thus in equation 3,5 dU = 0, and hence 
integrating over one period 

27r/a> 

0 

The heat produced per second and per unit volume, i.e. the 
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rate of loss L of energy from the electric field, is therefore given 

2njo) 

r OJ r ^dD 

^” 5 ? J 


E — dt. 
8t 


Inserting here E from equation 2.1 and D from 2.2 and 2.4, 


one finds 


wEl 

8n^ 


27r/cu 

I J ooswt^e 


3 COS cot 


3 sin (oi 


dt 


3t 




or after integration, 


L ^ - 


€2 El CO 


877 


(3.15) 


This can be expressed in terms of the phase-shift (/> (equation 
2.5) by 

(3.16) 


L = tan (}>. 


877 


For this reason (f> is usually described as the loss angle. 

Another derivation of equation 3.15 is given in the appendix 
(A l.ii). The connexion of and with the optical constants 
of the material is also discussed in the appendix (A. l.iv). 



CHAPTER II 

STATIC DIELECTRIC CONSTANT 

4. Survey 

In the present chapter the intention is to calculate the electric 
dipole moment induced by an external field in a dielectric from 
its atomic and molecular structure. The dielectric constant 
can then be obtained with the help of equation 1.9. 

In § 7 a formula of general validity will be derived which 
connects the static dielectric constant 6^ with structural proper- 
ties of the substance. The calculation of explicit values of €g 
and its temperature dependence is usually, however, beset with 
great formal difficulties and two types of approximations are, 
therefore, introduced. 

Firstly, a simple model is in general chosen to represent a 
material of much greater complexity. Secondly, mathematical 
approximations are often introduced which hold only within a 
certain range of some parameter such as the temperature. 
Consequently, in order to decide whether an approximate 
formula is applicable for a given material, one has to judge (a) 
whether the basic model can be chosen to represent the actual 
material, and (b) whether the mathematical approximation 
holds for the given range of the parameter concerned. 

A dielectric substance can be considered as consisting of 
elementary charges e^, and 

2^,-0 (4.1) 

if it contains no net charge. The electric dipole moment f of a 
charge relative to a fixed point is defined by the vector if 
is the radius vector from the fixed point to The total moment 

of the whole system is the vector sum of all the individual dipoles, 
^ 1^*. This quantity is independent of the position of the fixed 

point if the sum of the charges is zero. For using this fact (cf. 4.1) 
the dipole moment relative to a point at a distance b from the 
original one is 

t Also referred to as electric moment or dipole moment or moment. 
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It will now be assumed that in the lowest energy state (ground 
level) of the substance its dipole moment vanishes. Then if is 
the positional vector of the charge in the ground level, 

(4.2) 

i 

Therefore, if is the displacement of the charge from its 
equilibrium position in the ground state, 1^. = Then 




M(X) = 2 


(4.3) 


is the electric dipole moment of the substance for a given set of 

di.placeme„te X ^ (4.4) 


which has been expressed in abbreviated form by X. Many such 
sets X may lead, of course, to the same moment M. 

Very often it is useful to collect some of the elementary charges 
into a group forming an atom, a molecule, a unit cell of a crystal, 
or some larger unit. Let the jth unit of this type contain the 
8 elementary charges 

Xj r^2j’**5 

be an abbreviation for the set of all their displacements r^x,..., 


Then m(x,) == 2 (4-6) 

is the electric moment of this jth group of charges, and 

M(Z) = 2 (4-7) 

3 

where the sum extends over all the groups. The vector sum of 
their individual moments tn(xj) thus forms the total moment 
M(X). Our task is to find the average displacements, and hence 
the average electric moment under the influence of an external 
electric field. 

In order to obtain a preliminary idea about the average con- 
tributions of certain displacements to the electric moment we 
shall consider two cases, each of which has its characteristic 
type of displacement: 

Case (i) the displaced charge is bound elastically to an equili- 
brium position; 
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Case (ii) a charge has several equilibrium positions, each of 
which it occupies with a probability which depends on the 
strength of an external field. , 

The interpretation of case (i) is that on displacing the charge 
e, carried by a particle of mass m a distance r, a restoring force 
proportional to — r acts on the particle in a direction oj^posite 
to the displacement (hence the ~ «ign). Thus, if a constant 
external field f is applied {i ^ time) 


dt^ 


— a>o r -] f, 

m 


(4.8) 


where wqI27t denotes the frequency of oscillation, and ~-7ncoQr 
is the restoring force. Equation 4.8 can be written 




(r-r) 


--cog(r-r), 


(4.9) 


where ? (4.10) 

mcDg 

i.e. drjdt — 0. The charge e, therefore, carries out harmonic 
oscillations about the position r which thus represents the time 
average of its displacement, i.e. if C and 8 are constant 

r = r+Ccos(a>o^+S). 

The average electric moment is, therefore. 



(4.11) 


As an example of case (ii) consider a particle with charge e 
possessing two equilibrium positions A and B, separated by a 
distance b. In the absence of an electric field the particle has 
the same energy in each position. Thus it may be assumed to 
move in a potential field of the type shown in Fig. 4. If in 
equilibrium with its surroundings it will oscillate with an energy 
of order kT about either of the equilibrium positions, say about 
A. Occasionally, however, through a fluctuation, it will acquire 
sufficient energy to jump over the potential wall separating it 
from B. On a time average, therefore, it will stay in A as long 
as in B, i.e. the probability of finding it in either A or B is 

4960.11 ci 
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The presence of a field f will afifect this in two ways. Firstly, 
as in case (i), the equilibrium positions will be shifted by an 
amount r which for simplicity will be assumed to be the same 
in A and B. Secondly, the potential energies of the par- 
ticle in the two equilibrium positions will be altered because 



Fig. 4. Potential energy of a charged particle with two equilibrium positions. 
The dotted curve holds in the presence of an external field f. 

its interaction energy with the external field differs by e(bf), i.e. 

V^-Vb = e(hf). (4.12) 

The particle should, therefore, on the average spend more 
time near B than near A. Actually, since according to statistical 
mechanics, the probability of finding a particle with an energy 
V is proportional to 

_ ^-VBlkT 

— ^-VAlkTj^^-VBlkT 

are the probabilities for positions A and B respectively. They 
have been normalized in such a way as to make 

Pa+Pb == 1 (4.14) 

in agreement with the physical condition that the particle must 
be in one of the two positions. Thus from 4.12 and 4.11 

P^’^Pb = ^ (4.15) 
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It follows from the definition of the probabilities jpj and 
that if the condition of the system over a long time i-^ is con- 
sidered, the particle will spend a time (use 4.14) 


PAh = [I-UPb-Pa)]*'! 


in position A, and a time pj^t^ = position 

B. It has thus been displaced by the distance b from A to B 
during the fraction \{Pb-~Pa) The average 

moment induced by the field is thus 


leb{pB-PA)- 


(4.16) 


Hence if 6 is the angle between b and f, the projection of the 
induced moment into the field direction is, using 4.16 and 4.15, 
given ,ebfcos9ikT_i 

In most cases it is permissible to assume 

€bf<^H\ (4.18) 


for putting e = electronic charge, / — 300 volts/cm. — 1 e.s.u., 
b = 10"® cm. distance between neighbouring atoms in a 
molecule, and T — 300° (— room temperature) one finds 


ehf _ 4*8x10-1^X10-® XI 
kT “ " I*4xl0-i®x30() 


(4.19) 


Developing 4,17 in terms of ebfjkT, the average induced 
moment in the field direction is found to be 


(|e6)2 flkT-\-er, (4.20) 

where er is a term similar to those considered in case (i) which 
has been added to account for the elastic displacement. 

Often two charges +e and e are strongly bound, forming an 
electric dipole p. — ed where d is the distance between them. 
The above case (ii) then leads to the same result as that of a 
dipole (1 having two equilibrium positions with opposite dipole 
direction, but with equal energy in the absence of a field. In 
a field £ the energy of interaction between field and dipole is 

-M. (4.21) 

SO that 2p/cos 0 is the energy difference between the two positions 
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if 6 is the angle between jjl and f . This is equivalent to equation 

4.12 if It.. 

jx = |eb. (4.22) 

Actually putting an immobile charge — e half-way between A 

and B turns case (ii) into the present case. Clearly the induced 

moment must be the same for both cases because the charge 

— e is immobile, and its distance from A and B is leading to a 

dipole moment /x, 4.22. Introducing this into 4.20 yields for the 

induced moment in the field direction 


fjL^ gohW 


f+er. 


(4.23) 


In contrast to case (i) the electric moment now depends on tem- 
perature. In view of equation 1.9, a substance consisting of a 
great number of such dipoles will have a temperature-dependent 
dielectric constant in contrast to a substance in which all charges 
are bound elastically. According to equation 3.12 this means 
that in the dipolar case (ii) the entropy of the substance is 
decreased by the field. This is evident because the field causes 
the fraction of dipoles with components in the field direction 
to be larger than the fraction of dipoles with components in 
the opposite direction, thus leading to a state which is disordered 
to a smaller degree (i.e. having a lower entropy) than the state of 
complete disorder in which p^ ~ pj^. 

The difference between the action of the field in the two cases 
(i) and (ii) should be well understood because this is essential for 
the whole theory of dielectric constant. In case (i) the field 
exerts a force on an elastically bound charge, thus shifting its 
equilibrium position. In case (ii) this force of the field on the 
charge again leads to contributions of type (i) denoted by er in 
equations 4.20 and 4.23. It would be wrong, however, to assume 
that the field by this force turns a dipole from one equilibrium 
position into another. This is effected in a more indirect way 
because the field slightly alters the probabilities of a jump of a 
dipole from one equilibrium position to another. This will be 
described in greater detail in § 9, where the dynamic properties 
of the present model will be investigated. 

It should also be realized that though every charge is displaced 
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elastically (case (i)) the fraction of dipoles turned by a field of 
reasonable strength (case (ii)) is very small. This fraction is 
given by which, according to 4.15, 4.18, and 4.19, is of 

the order 1()~^ in a field of 300 volts/cm. Even in a field of 
100,000 volts/cm. only about 2 per cent, of all dipoles are 
orientated. 

5. Dipolar interaction 

In a dielectric two essentially different types of interaction 
forces should be distinguished. Forces due to chemical bonds, 
van der Waals attraction, repulsion forces, and others have all 
such short ranges that usually interaction between nearest 
neighbours only need be considered. Compared with these 
forces dipolar interaction forces have a very long range. This 
can be readily shown as follows. 

As indicated previously (cf. § 4) a polarized dielectric can be 
considered as composed of small regions each having a certain 
dipole moment, and the total dipole moment of the body is the 
vector sum of the moments of these regions. Now it is well 
known from macroscopic theory that the energy per unit volume 
of a macroscopic specimen depends on its shape (cf. Appendix 
A 2.iii). This implies that interaction between dipoles must be 
taken into account even at macroscopic distances and illustrates 
the great importance of the dipolar interaction forces. 

Due to the long range of the dipolar forces an accurate cal- 
culation of the interaction of a particular dipole with all the 
other dipoles of a specimen would be very complicated. How- 
ever, very good approximation can be made by considering that 
the dipoles beyond a certain distance, sayf can be replaced 
by a continuous medium, having the macroscopic dielectric 
properties of the specimen. Thus the dipole whose interaction 
with the rest of the specimen we are calculating may be con- 
sidered as surrounded by a sphere of radius containing a 
discrete number of particles, beyond which there is a continuous 
medium. To make this a good approximation the dielectric 
properties of the whole region within the sphere should be equal 
t The suffix m stands for macroscopic. 
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to those of a macroscopic specimen, i.e. it should contain 
a suflScient number of molecules to make fluctuations very 
small. 

From this point of view we are thus led to H. A. Lorentz’s 
[Z/3] method for the treatment of dipolar interaction: from a 
macroscopic specimen select a microscopic spherical region 
which is sufficiently large to have the same dielectric properties 
as a macroscopic specimen. The interaction between the dipoles 
inside the spherical region will then be calculated in an exact way, 
but for the calculation of their interaction with the rest of the 
specimen the latter is considered as a continuous medium. To 
demonstrate this method we shall make use of a very simple 
model. It consists of a cubic lattice of 'atoms’ whose linear 
dimensions are very small compared with the lattice distance. 
Each atom consists of a positive charge -fe which is rigidly 
bound to the lattice point, and of a negative charge — e which is 
bound elastically to it. The force acting on a negative charge 
when it is displaced while all others are at rest at their respective 
equilibrium positions will be denoted as the restoring force. When 
other charges are displaced as well, an additional force will act 
on the charge due to a change in the interaction between the 
charges. It will be assumed that electrostatic interaction only 
exists, and in particular that there are no short-range forces. 
Furthermore, the temperature will be assumed to be so close to 
the absolute zero, T ~ 0, that thermal oscillations can be dis- 
regarded. 

On these assumptions when a macroscopic electric field E is 
set up in the substance all negative charges are displaced by the 
same amount, say r, each forming a dipole 

m — (—e)r. (5.1) 

The field f acting on one charge and displacing it against the 
restoring force is often referred to as the local field, or the inner 
field; it has to be distinguished from the macroscopic field E. 
Thus if the restoring force is denoted by —c^r. 
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and hence according to 5.1 

m = (5.2) 

To calculate f according to Lorentz’s method, separate the 
sources of f into those inside the spherical region and those out- 
side it, making contributions and respectively, 

f - (5.3) 

Since all induced dipoles are identical, the field f must be the 
same at every lattice point. To find consider therefore the 
dipole at the centre of the spherical region and calculate its 
interaction energy I with all the other dipoles of the region. 
This will depend on f , and 

(— e)fi = —grad /(f). (5.4) 

To calculate this interaction energy I assume the field to be 
sufficiently weak to make f very small compared with the lattice 
distance a^. Each dipole m can then be considered as a point 
dipole. The electrostatic interaction energy between two parallel 
dipoles of equal moment is given by 

^(l-3cosV), (5-5) 

where I is the distance between them and 0 is the angle between 
1 and m. For simplicity let m be parallel to a crystal axis, say 
in the ^-direction. Then the three components of 1 are 
X = TOq, Y = Z — qa^ 

where n,p, q are positive or negative integers. Since cos 0 = Zfl, 
the total energy of interaction follows from 5.5 by summation 
over all lattice points, i.e. 

T _m^X 

Now to each set of three values ?^, q, say n == Uq, p Pq, 
q ^ qQ, two others can be coordinated by cyclic permutation, 
namely n ^ p^, p = q^, g = and n q^, p ^ q p^. 
These three terms just cancel in the above sum, and hence / = 0. 
It thus follows, using 5.4, that 


(5.6) 
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The external contribution/^ has to be calculated macroscopi- 
cally. It represents the electric field inside the spherical region 
due to all sources except the polarization inside this region. To 
understand this clearly it should be remembered that the 
macroscopic electric field E is partly due to true charges outside 
the specimen (or at its surface) and partly due to the polariza- 
tion P of the dielectric which acts in the opposite direction. 
To obtain the contribution to E due to the spherical region 
should be omitted. Therefore if E^ is this contribution, 

f, E-E,. (5.7) 

By definition, E^., the self-field, is the field inside a spherical 
specimen with a permanent polarization P. Its direction is 
opposite to P. A simple electrostatic calculation (Appendix 
A 2.21) shows that 

E«=-yP. (5.8) 

Therefore, using 5.3, 5.6, 5.7, 5.8, and 1.0, 

477 ^ 2 

f = = E + = (5.9) 

o o 

It should be noted that this expression is independent of the size 
of the spherical region, which means that the contribution of any 
homogeneously polarized spherical shell to the field inside it 
vanishes. In other words, the interaction between such a shell 
and a dipole inside it vanishes. This is in agreement with the 
above result 5.6 obtained by considering dipole-dipole inter- 
action in detail. The present model, therefore, has been chosen 
in such a way as to make macroscopic and microscopic treat- 
ment identical. Such a result, of course, could only be obtained 
by assuming short-range forces to be entirely absent. The exis- 
tence of such a force would modify f^, but it would not alter our 
result for f^. 

The fact that the interaction energy of the dipoles of a spherical 
region vanishes must mean that if a spherical specimen is brought 
into a homogeneous field Eq in vacuum this field represents the 
local field acting on each dipole. Actually from electrostatics 
it follows that the field inside a spherical specimen of dielectric 
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constant is E “ 3Eo/(eg+2) (cf. Appendix A 2.16), and hence 
with equation 5.9, f = Eq. 

The above result 5.9 can be derived by an alternative method 
due to Onsager |Oi] which is capable of a generalization and 
which will be of use at a later stage of the development of a 
general theory. 

The reader should realize, however, that in Onsager ’s paper 
the spherical region contains one single dipolar molecule only. 
This is, of course, open to the objection that the spherical region 
is certainly not large enough to have the same properties as a 
macroscopic specimen. Tn our use of his analysis, however, the 
spherical region is sufficiently large to make it certain that it is 
correct to treat the outside macroscopically. Now consider, 
in the absence of a macroscopic field, a spherical region which is 
polarized homogeneously, thus having a dipole moment M — PF 
where V is its volume. If considered separately without the rest 
of the specimen the electric field inside the spherical region is the 
self-field E^, equation 5.8. If, on the other hand, the spherical 
region is considered inside the specimen there will be a certain 
interaction between the polarized sphere and the surroundings, 
giving rise to an altered field inside the spherical region. Its 
difference from E^. is denoted as the reaction-field R and is given 
by (Appendix A 2.18; — radius of the spherical region) 

^ ^ 2(e,~l) M ^ 2 e 47 tM ^ 2 e,-l 
26^+1 3 2eg-|-l V 3 26^+1 

(5.10) 

Thus R is the field produced inside the spherical region by the 
surroundings if polarized by the former. If a macroscopic field 
E is now produced in the specimen without altering the moment 
M of the sphere, the field inside it will be increased by an addi- 
tional field, the cavity-field G (cf. Appendix A 2.15), 



Thus the total field inside the sphere, due to outside sources, is 
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This expression holds whatever the value of P. In our particular 
case if irrP — (e^— 1)E (cf. § 1) is introduced, the above expression 
becomes identical with 5.9. 

Finally to calculate €g, insert f from 5.9 into 5.2. The polariza- 
tion P is then given by 


whore iVJj is the number of particles per unit volume, 
using 1.9 

fiZZ— 

^ 8+2 


3 


Hence 

(5.13) 


which is usually known as Clausius-Mossotti formula (cf. 
Clausius, 01; Mossotti, M5). The above derivation of this 
formula is exact except for the assumption that r is a small 
quantity. It will be noted, however, that the change in as r 
increases is essentially a manifestation of the field dependence 
of the dielectric constant (since the change in magnitude of f 
depends only on field strength) and this docs not concern us here. 
That is, formula 5.13 is correct for the limiting case in which the 
dielectric constant is independent of the field. The reader should 
realize, however, that this derivation holds for the above model 
only; there are few substances which it could claim to represent, 
though for simple non-polar substances it will often be a useful 
approximation. 


6. Dipolar molecules in gases and dilute solutions 

Molecules can be divided into two classes, polar and non-polar, 
according to whether or not they possess an electric dipole 
moment in their lowest energy -level (ground state). In general 
it is easy to recognize the class to which a molecule belongs 
because a non-polar molecule must have a point of symmetry 
defined in such a way that the distribution of charges along (or 
near) any straight line passing through it must be symmetrical 
with respect to this point. Thus a diatomic molecule is polar 
unless its two atoms are equal (e.g. the polar molecules HCl, 
CO, and the non-polar molecules Hg, Og). Triatomic molecules 
of the type ABg, where A and B represent different atoms, are 
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polar unless their nuclei lie on a straight line with A half-way 
between the two B atoms. Examples are the triangular polar 
H 2 O and the straight non-polar COg molecules. CgHg is a more 
complicated non-polar molecule forming a plane hexagon with 
the centre as point of symmetry. On replacing one H atom by 
another type of atom, say by Cl, the resultant molecule C^HgCl 
becomes polar. 

The magnitude of molecular dipoles is usually of the order 
of one electronic charge e.s.u.) displaced by 

J X 10~® cm., i.e. about 10“^® c.g.s. units. Often dipole moments 
are expressed in Debye units, one such unit being c.g.s. 

units. Debye was the first to recognize the importance of the 
investigation of dipole moments for a study of the constitution 
of molecules (cf. his book, reference D2). A detailed discussion 
of molecular structure is, however, not intended in the present 
book. Our purpose in this section is to determine the properties 
of dielectric substances consisting of molecules having per- 
manent dipole moments and our model for such a molecule (in 
vacuum and free from perturbing influences) will, therefore, 
have a dipole moment 

In addition to its translational motion, a free molecule can 
carry out oscillations and rotations. However, unless stated 
otherwise, it will be assumed that these do not alter the average 
value of the dipole moment. In a complicated molecule intra- 
molecular rotation of dipolar groups may lead to considerable 
changes in the dipole moment of the molecule as a whole (cf. 
the book by Le Ffevre, L2 ) . Such molecules will not be considered 
at present. 

A constant electric field f will influence the molecule in two 
ways. Firstly, it will perturb the free rotation of the dipole, 
and secondly, it will induce a further dipole moment, say af, 
by elastic displacement of the atomic electrons relative to their 
respective nuclei, and to a smaller extent by elastic displacement 
of the nuclei relative to each other. The total moment of the 
molecule is thus 

m = (6.1) 

The quantity a has the dimensions of a volume and is called the 



28 


STATIC DIELECTRIC CONSTANT 


II, §6 


polarizability of the molecule. In anisotropic molecules (that is, 
molecules having different polarizabilities along different axes) 
the induced moment need not always have the same direction 
as the field f ; the polarizability a is in this case not a scalar but 
a tensor quantity. Its average value, a, obtained by allowing 
the molecule to have all possible directions relative to the field, 
is, however, a scalar. Assuming that the polarizability of a 
molecule is entirely due to electronic contributions, a is deter- 
mined by the optical refractive index n. For, according to Max- 
well’s relation, is the dielectric constant at optical frequencies 
at which there are no dipolar contributions (because the time 
required by the dipoles to attain equilibrium with the field is 
much longer than the period of the field). Assuming the mole- 
cules to be isotropic, it will be shown below that the Clausius- 
Mossotti formula holds approximately in this case, i.e. 


OL ~ 





(6.2) 


where a is the radius of a sphere which on an average contains 
one molecule. 


Oases 

We shall now proceed to calculate the dielectric constant 
of a gas of dipolar molecules. In order to simplify calculations, 
it will be assumed that the density of the gas is so low that the 
dipolar interaction energy is small compared with the thermal 
energy (c^ kT per molecule) and can therefore be neglected. 

According to § 5 the dipolar interaction energy is of the order 
/x^/Z^ if is the number of molecules per unit volume. 

It will thus be assumed that 

lxlN^<^kT. (6.3) 

This simplifies the calculation of considerably because it 
implies that the local field f, acting on a dipole, is entirely due to 
external sources, i.e. f = D = e^E (cf. § 1). Also if Nq is small 
€g— Tmust be small as well, so that the further assumption 

6,-1 <1 (6.4) 

leads to f == E. (6.5) 

It will be found below that 6.4 follows from 6.3. 
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In order to calculate the electric moment M of the gas we make 
use of the fact (§ 4) that M is equal to the vector sum of the 
moments m of all N molecules. Since, on the other hand, the 
time average in of m is the same for all molecules, 

M “ Nin. ( 6 . 6 ) 

Thus from 1.9 ^ 5 — 1 = ^TrNQmjE, (6.7) 

where according to 6.1 and 6.5 

m == (6.8) 

In this expression the first term represents the average value of 
the intrinsic moment, and the second one is the average induced 
moment per molecule. The calculation of m is similar to that 
carried out in case (ii) of § 4. In the present case, however, 
instead of having only two possible directions as described in 
§ 4, a dipole has a whole continuum of possible directions. As 
in § 4, the behaviour of the dipole may be considered on a 
statistical basis, without inquiring into its dynamic properties 
under the influence of the field E. Thus since —EfjL^coBd is 
the energy of the dipole in the field, the probability of finding 
in a direction forming an angle between 6 and d-\-d6 with E 
is given, according to statistical mechanics, by {2rrBinddd is the 
element of the solid angle between 6 and 6-\-d9) 

cos Ojk T g d dd j { ^ sin d dd, ( fi • ^ ) 

' 0 

provided that thermal equilibrium has been attained. At 
present in considering static properties this will always be 
assumed to be the case. 

In evaluating 6.9 we shall assume, as in the case of equation 
4.18, that the field is sufficiently weak to make 

^< 1 - ( 6 . 10 ) 

Then developing 6.9 in terms of EjkT, and keeping terms up 
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cos 6, 

'Tt . Tt 

cos ^ — j 0^Efi^coB0lkT gjjj 6 dd ^ ^EfL^coBdlkT 6 dO = . 

0 ' 0 

( 6 . 11 ) 

Similarly the components of perpendicular to E are found 
to vanish, so that the average moment of a molecule is a vector 
with the same direction as E, amounting to — fjL^ cos 6, i.e. 

iI. = 4,E. (6.12) 

On using 6.7 and 6.8, therefore, 

+ if ^,-1 < 1 (6.13) 

or introducing Coo = i + ^TraiVJ) (6.14) 

as dielectric constant at frequencies so high that the dipoles have 
no time to attain equilibrium, 

= (6.15) 


In agreement with the conclusions drawn in § 4 the dipolar 
contribution to is temperature-dependent in contrast to the 
non-polar contribution €oo. Measurement of the temperature- 
dependence of thus makes possible the separation of the dipolar 
contribution, and hence the determination (examples 

in § 16). 


Dilute solutions 

The above derivation suggests that a formula similar to 6.16 
should hold for dilute solutions of dipolar molecules in a non- 
polar liquid having the temperature-independent dielectric 
constant For at sufficiently low densities the interaction 
between dipoles can again be neglected. Also, replacing 6.4 by 

€,-eo<l, (6.16) 


t For with x = cos 0 and y ~ ii^ EjkT 1, 

008 0 ==s J xsy^dxj j e’^^dxQt j {x-\-YX^)dx j j dx = 


0 “}- 
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equation 6.5 holds as before, leading to 6.11 because in a liquid, 
in the absence of an external field, all directions of the dipole 
are equally probable. Furthermore, if is again the non-polar 
contribution to e^, then {€^—\)El4:7r represents the non-dipolar 
contribution to the electric moment of the solution per unit 
volume. Since is the total moment, {e^—€^)El4:7T, 

in view of the additivity of moments, is, as before, the dipolar 
contribution. Assuming that the dipole is rigid, i.e. that a = 0, 
equation 6.15 would follow with eg. Actually a molecular 
dipole is not rigid, but can be polarized. In a solution this leads 
to a change in the effective dipole moment of the molecule, as 
was recognized by various authors (cf. Weigle, W 1\ Frank, Fl\ 
Higasi, H2\ Frank and Sutton, F4\ these authors also consider 
other effects which may lead to an alteration in the effective 
dipole moment when the molecule is non-spherical or has a 
large quadripole moment). The dipole polarizes its surroundings 
which in turn produce a reaction field at the position of the 
dipole. This field polarizes the molecule and thus alters its dipole 
moment. Thus if we define the resultant moment which the 
molecule has in the solution as the 'internal moment’ jx^, the re- 
action field will be proportional to say In the absence 
of an external field, therefore, from 6.1 (using f == g\Li and 


m = jjLj) 

Hi = 

(6.17) 

or 

1 — CX^ 

(6.18) 


A quantitative calculation of the reaction field, and hence 
of (X^, is very difficult, however, except in the case of spherical 
molecules. 

It follows from the above that, when a field E, produced by 
external sources, exists in the solution the resultant field acting 
upon a molecule in the solution will be different from E. Also 
the moment of the molecule in vsolution will be different from 
that it has in vacuum (namely In order to take account of 
these facts in calculating the dielectric constant of the solution 
we may consider the dipole in one of two ways. In the first of 
these we calculate the field acting upon the 'internal moment’ 
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This will, of course, depend on the shape chosen for our model 
of the molecule. In the second method we calculate the inter- 
action between the field of the dipolar molecule in the solution 
around it and the applied field E. The same dipolar field might 
alternatively be produced by a rigid dipole having in vacuum a 
moment fx^ which will be denoted as the ‘external moment’ of 
the molecule. f A solution of rigid dipoles fx^ should, therefore, 
have the same dielectric constant as the solution which we 
actually consider. 

The former method will be used below in deriving the Onsager 
formula. At present the latter method will be chosen because it 
enables us to make use of our previous derivation for gases to 
determine the dielectric constant of the solution, i.e. to employ 
equation 6.15 after replacing fx^ by [x^. Thus with the use of 6. 1 8 

_<r — 

* “ 3/cT ZlT \l-ixg) ’ 

if e^— and e^— €o<l. (6.19) 

Here the ratio of external to internal moment depends on 
the shape of the molecule. For spherical molecules, according 
to the definitions of and p,:, the internal moment p^; is identical 
with the moment of a sphere inside the dielectric containing p^ 
at its centre. An ordinary calculation in electrostatics shows 
that in this case (cf. Appendix A 2.31) 

and if a is the molecular radius (cf. Appendix A 2.19) 


2(e,-l) l 

26,+ 1 


(6.21) 


Furthermore, if the polarizability is isotropic and is mainly due 
to displacement of electrons, equation 6.2 holds, n being the 
refractive index of a pure liquid of the dipolar molecule. Then 
with 6.18, using 6.21 and 6.2, the internal moment of a spherical 

molecule is -2.,+ l«>+21 




/2«,+ l »2+2\ 


( 6 . 22 ) 


Compare Appendix A 2.ii for a discussion of dipole moments. 
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e.(»=>+2) 




(6.23) 


The dielectric constant in this case is (cf. 6. 1 9) given by (replacing 
€g by 6 q on the right-hand side according to 6.16) 


( €g(n^+2)Y 


./fo^ 
12^0 
rrpgiV'o / eo+2 \^( 
3/cT \ 3 / ( 


UT 

47rp*iV'o 


(6.24) 


2 ( eo - l )( 6 o -««))2 


(2eo+n^)(eo+2)| 

for spherical molecules and 


< 1. 


In some cases we might expect the approximation of spherical 
molecules to give a correct order of magnitude for the deviation 
from unity of the ratio of the dipole moments in vacuum and in 
solution, but we should not expect to obtain a more accurate 
result than that. 


Onsager formula 

In the case of spherical molecules Onsager \OT] has shown 
that it is possible to go one step farther in the approximate 
calculation of the dielectric constant. The interaction between 
molecules will no longer be entirely neglected, but one com- 
ponent of it, namely the long-range dipolar interaction, will be 
taken into account. The following assumptions will, therefore, 
be made: 

(а) A molecule occupies a sphere of radius a and its polariza- 
bility is isotropic; 

(б) The short range interaction energy is negligible (i.e. 
<^kT per molecule). 

Assumption (6) means that the surroundings of a molecule 
will be treated as a macroscopic continuum with dielectric 
constant because long-range forces only will be considered. 
Estimates on the range of validity of this method can be made 
(of. reference F12), but will be postponed until the more general 
Kirkwood formula is derived (§ 8). 

On the assumptions made above, the contribution of a single 
molecule to the dielectric constant can be calculated in the 
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same way for a pure liquid or for a mixture. In both cases the 
local field f acting on the molecule is the field inside a cavity of 
radius a within a continuous medium of dielectric constant 
and is composed of (i) the cavity field G due to external sources, 
and (ii) of the reaction field R due to the moment m of the mole- 
cule itself. Both G and R were discussed in § 5, although there 
the cavity of radius d contained a great number of molecules. 
To use the same values for G and R in the present case is only 
possible because of the assumptions (a) and (6). Thus, according 
to 5.10, 5.11, and 5.12 (M/a^ of § 5 corresponds now to mja^) 
and using 6.21 

f = G+R = -J^E+srm. (6.25) 

Inserting this value into 6.1 the moment m becomes 


m = fx. 


3e, 


2 ^ 5+1 




or solving with respect to m. 


m = 


. 3 €, oE 

l-‘Ocg'^2€^+l l — ocg^ 


(6.26) 


With this value for m the internal field f (6.25) becomes 


3c, E 
2c,+ l 1 — ocg 


g 

l—ocg 




(6.27) 


To calculate the average polarization m we shall require 
Considering that 




3eg cos e g 2 

2eg+l \—0Lg 1— 


(6.28) 


is the energy of a dipole in the field f , the probabihty of finding 
it in a direction forming an angle 6 with E is (similar to 6.9) given 

by 


sinOdO^ j sin 0 dd, 


(6.29) 


where should be introduced from 6.28. It will be seen that 
the second term in 6.28 is independent of 6 and therefore in 6.29 
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only the first term remains. Following the same procedure as 
in the derivation of 6.11 one finds 


COS0 = 






E 


and hence 


\^v =: 


2 €,+ 1 ^m\l-oLgy 

36, 




26jj+1 ‘ikT{\ — <xg) 


E. 


(6.30) 


Thus from 6. 26, using 6.30 (and a — oL for isotropic polarization), 
the average moment m of the molecule is given by 

3e, 


Let us consider first the application of 6.31 to the simple case of a 
pure liquid of non-polar molecules for which = 0. Since by 
definition a molecule occupies an average volume of 


47r . 1 

— a'^ — — , 

3 Ny 


equations 6.7 and 6.31 yield 


€„-l = 


3c, 


(6.32) 


(6.33) 


2€^.+ 1 1— 

Introducing g from 6.21 and solving with respect to a leads to 
the Clausius-Mossotti equation 


— 1 a 

€^+2 


(6.34) 


This equation can be considered to prove 6.2, which, together 
with 6.21, transforms 6.31 into 


3e,E f 2e,+ l , r^^-1 J 

2€,+n^\2e,+n^\ 3 ) ZkT ''3 /' 


(6.35) 


The Onsager formula for a pure dipolar liquid is then obtained 
by inserting 6.35 into 6.7, making use of 6.32: 


e*— W® = 


2e,+n^ 3kT \ 3 } ‘ 


(6.36) 


For very small densities this expression becomes identical 
with 6.15, as should be expected (because for Cg— 1 ^ 1 and 
Ti^—l < 1, (w2-j-2)/3 1 and 3cg/(26g+l) 1). 
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Consider next a mixture containing mole- 

z 

cules per c.c. of z different compounds. Then is the 

total moment per c.c. where is the average moment of a 
molecule of the type s obtained from 6.35 by inserting for 
n, and a the respective values and which these quantities 

have for the 6’th kind of molecule. From 1.9 the dielectric con- 
stant of the mixture thus follows: 


eg— 1 = 477 2 nhlE. (6.37) 

As a particular case take a mixture containing per c.c. 
polar (/x,yi, Uj) and Ag non-polar (fi^^ ^ ^ 2 ’ ^ 2 ) molecules. 

Inserting 7\ and mg from 6.35, equation 6.37 leads to 
/ 477^^1 iVi / 2eg+l n\-\-2 \ 

ZkT \2eg+wf 3 ) 


1 3eg (4 
* 2eg+ll 


+ 


2 €.+ 1 


(nf- 




2e^+nl 


(nl-1) 


faIN, . 


(6.38) 


For a dilute solution of polar molecules, i.e. when Ag, 
this expression should become identical with 6.24 ifn~ 
n\ = €q, and — Nq. Actually since represents the 

volume occupied by the non-polar molecules in 1 c.c., N-^ 
means that 477a| A 2/3 1 . Since the second term in 6.38 can be 

neglected, we find with the above substitutions for and 

, . 3eg 3eg 4,7^§i / 2^^+ 1 

^ 26g+€o^ « ^■^2€g+l %kT \2e,+n^ 3 ) ’ 

which, after a simple transformation, becomes identical with 
6.24 (since e^). 


7. General theorems [FIO] 

The formulae derived in the previous section for the dielectric 
constant have the great advantage of representing €g in a very 
simple way with the help of a few parameters. It should not be 
forgotten, however, that they only hold subject to certain con- 
ditions which in many practically important cases are not ful- 
filled. In the present section, therefore, expressions for the 
static dielectric constant will be derived which hold in a very 
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general way for any dielectric substance which is not per- 
manently polarized. Clearly the derivation of such expressions 
requires some mathematical abstractions, but in view of their 
importance it was thought desirable to give a detailed account 
of their deduction. 

As in § 5, select from an infinite homogeneous specimen a 
macroscopic spherical region of volume V which is large com- 
pared with a region just big enough to have the same dielectric 
properties as a macroscopic specimen. The surface of the 
spherical region need not be an exact geometrical sphere; 
deviations of a molecular magnitude are permitted and the 
exact surface will be laid in such a way that no molecules are cut 
by the surface. This has no influence on the field at distances 
which are large compared with atomic dimensions. We shall 
calculate the projection Mj^ of the average electric moment of 
the sphere in the direction of the macroscopic field E. For this 
purpose all particles inside the sphere will be treated according 
to the rules of classical statistical mechanics. The outside, 
however, will be considered as a continuous dielectric described 
by the macroscopic dielectric constant €^. It will be assumed 
throughout that the macroscopic field E is sufficiently weak to 
prevent saturation, so that is independent of E, 

The spherical region consists of a number of elementary 
charges e^, each of which can be described in terms of its displace- 
ment from the position it would have in the lowest energy -level 
(ground state) of the whole system. This displacement is a vector 
quantity and is denoted by r^. A set of all the displacement 
vectors will be collectively denoted by X, according to equation 
4.4. Except at the absolute zero of temperature, such a system 
of particles does not stay in the same configuration, even if 
macroscopically it is in equilibrium. Owing to thermal fluctua- 
tions there is a probability 

^-v(x.mikT I j ( 7 . 1 ) 

of finding it with any set of displacements in a space element 
between 

= (fi. Ta.-. r<.-) and X+dX = (ri+dr,,..., r<+dr<,.,.). 
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Here, U (X, E) is the potential energy of the system in the con- 
figuration X, in the presence of a field E, and 

dX = (7.2) 

is the product of the volume elements dr^ of all the displacements 
each being again a product of its components dr^^dr^ydr^. 
The integration must be carried out over all possible values of all 
the displacements. 

Each set of displacements X leads to a dipole moment M(X) 
as shown in equation 4.3. Therefore, if 0 is the angle between 
M(X) and E, Jf(X)cos0 represents the projection of M(X) in 
the direction of the macroscopic electric field E, and is 
the average value of ilf(X)cos0. Thus using 7.1, 

Mjs = ^ (Z)cos e e- uixMkT dxj^ e- mx.mkr gx. (7.3) 

The value U (X, 0) of the energy in the absence of a macro- 
scopic field will be denoted by U(X), and the zero of potential 
energy will be chosen so that U(X) vanishes in the ground state, 
i.e. when all displacements vanish. The energy U(X) can be 
considered as composed of (i) the energy of interaction ?7^(X) 
between the particles of the spherical region, and of (ii) their 
energy of interaction U^,{X) with the external region, 

U{X) = U,(X)+U,(X), (7.4) 

U^(X) will depend not only on X, but also on the dielectric 
constant because the external region is to be treated on a 
macroscopic basis. 

The fact that U^(X) contains the parameter which may 
depend on temperature requires some consideration. Our 
present investigation is a special case of a wider group of pro- 
blems — the statistical mechanics of systems containing tempera- 
ture-dependent parameters — ^which has been investigated by 
Gross and Halpern [Gl], Applied to our case their investigations 
show that €g has to be treated as constant parameter whenever 
differentiations with respect to temperature are required; 
U^(X) is then the energy required to establish at constant 
temperature the equilibrium polarization in the outside region. 
From the point of view of the whole system (spherical region 
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+outside region) U^{X) is thus the corresponding free energy 
composed of the energy of interaction of the charges inside the 
sphere with the polarization induced by them in the outside 
region, and of the free energy of this region (so far as it depends 
on X), as shown in the subsequent example (cf. 7.17 and 7.18). 
This is reasonable because ?7(X) is the energy required to estab- 
lish the displacement X at constant temperature. 

Now for a given configuration X, assume the macroscopic 
field E to be applied (which, of course, alters the probability 
of this configuration). Inside the spherical region this leads to 
an additional homogeneous field G (the cavity field) given by 
equation 5.11. It follows from electrostatics that its interaction 
with the charges of the spherical region is given by 

-M(Z)G = --^^M(X)EcoHe. (7.5) 

Thus in the presence of a field E, 

U{X,E)=U{X)--^^M{X)Ecose. (7.6) 

This expression for f7(X, E) must be inserted into equation 7.3. 
It should be remembered now that E was to be considered as 
sufficiently weak to prevent saturation. Thus if the right-hand 
side of 7.3 is developed into a power series in E, only the first 
term need be considered, i.e. 


^-U(X,E)lkT „ 


Q-mX)}kTll 


Se, M(X)E 
26,+ 1 kT 


COS0+... 


)• 


(7.7) 

Now Jilf(X)cos0e-t^<^)'*^eiX = 0 (7.8) 

because this integral is proportional to the average moment in 
the E-direction in the absence of a macroscopic field. Therefore, 
inserting from 7.7 into 7.3, making use of 7.8, 

= 2^ % J e-uixmT^x, (7.9) 

where 1/J = J (7.10) 

Assume now that the field E may have any direction relative to 
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a given M(Z), and average over all these directions. Then 
coqW will be replaced by its average value, i.e. by Inserting 
this value into equation 7.9, and making use of 1.9, yields 




47 t 3e„ 


where 


3F 26,+ 1 kT^ 


(7.11) 


(7.12) 


is the average value of M^(X) in the absence of a macroscopic 
field. Equation 7.11 represents our first general result. It 
shows that the dielectric constant can be expressed in terms of 
the mean square of the spontaneous polarization of a sphere 
of the dielectric embedded in a large specimen of the same 
material. 

Before developing the general theory further, equation 7.11 
will be shown to be self-consistent. This means that 7.11 must 
be fulfilled identically if the spherical region is treated on 
a macroscopic basis. For this purpose tlie following theorem 
will be employed (cf. reference Tl), Suppose the free energy 
F((Xi, aav) of a system depends on a number of macroscopic 
parameters a^, ag,.--- Then the probability of finding the system 
in a range is given by 




Also it will be noted that the quantity {doL-^doL^^.,) in the above 
expression represents the volume of the 'a-space’ element 
between ((x^, ag...) and (ai+c^aj, ag+dag,...). In the present case, 
since the free energy depends on the absolute magnitude of M 
only, this element is the sheU between the two spheres of radius 
M and M--\~dM in the Jf -space, and has the volume ^irM^dM, 
Since M may in theory have any value between zero and infinity, 
the probability that the spherical region will have a moment 
between M and M+dM is 

00 

e-PmkT^z^Mj f e-^(^^^M^dM. (7.13) 

0 

00 00 

Hence ^ = J dM j j M^dM. (7.14) 

0 ^ 0 
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Comparing this with equation 7.12, and remembering that 
according to 3.3 the entropy 8 is given by 8 = (U—F)IT, it 
follows that dXr>.e^l’‘M^clM (7.15) 


represents the number of states of the spherical region, in a range 
between M and M -{-dM, This equation might have been used as 
an alternative starting-point for the macroscopic considerations 
based on 7.12. 

Following the procedure used in the case of the energy U(X), 
the free energy F{M) can be considered as composed of an in- 
ternal free energy F^{31), and an external one, F^{31), due to the 
interaction with the surroundings. The former is the self-energy 
Fg calculated in the appendix (A 2.37), 




e,+ 2 

3V e,-l’ 


( 7 . 16 ) 


the latter can be obtained with the help of the reaction-field R 
(cf. equation 5.10). It is composed of the energy — MR of the 
dipole M in the field R, and of the free energy required to polarize 
the surroundings accordingly. This latter should be propor- 
tional to say Then 

F,{M) = -MR+|8i?2. (7.17) 


To find j8, Ff. can be considered to depend on the parameters 
R^, Ry, [for a fixed value of M\ and hence since for a system 
in equilibrium the free energy should be a minimum 


or 


dF,, 

8R^ 


0 , 


8F, 

dR.. 


0 , 


2i3R = M. 


8F, 

dR, 


= 0 


Inserting this into 7.17, using 5.10, 
F,{M) = -MR+^MR = -|MR 


Thus from 7.16 and 7.18 


47rif2 e,-l 
3F 2e,-M- 

(7.18) 


F = F,{M)+F,{M) 


V (2e.-fl)(e,-l)’ 


(7.19) 
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and hence withf 7.14, 

IP = (7.20) 

Ztt oeg 

which shows that equation 7.11 actually is fulfilled identically. 
The derivation of this result demonstrates the importance of the 
long-range interaction terms — included in FJ^M ) — in equation 
7.11. This means that in view of this interaction the fluctuations 
of a sphere embedded in its own medium are different from 
the corresponding fluctuations for a sphere in vacuum. In the 
latter case F ~ F^, and 7.14 leads to 

(7.21) 

Ztt €g+2 

Comparing 7.20 with 7.21 shows that becomes relatively 
large for substances with high dielectric constant in contrast to 
which does not depend appreciably on if ^ 1. 

In returning now to the further development of equation 
7.11 a method will be used which forms a generalization of the 
method employed by Kirkwood {K4) in the case of liquids con- 
sisting of rigid dipoles. Assume that the spherical region is 
composed of N molecules or other groups of atoms in such a way 
that each group has the same average polarization in an external 
field. In a pure liquid, for instance, each group contains one 
molecule; in a crystal it contains all the particles of a unit cell. 
The spherical region can then be divided into N ‘units’ each of 
which makes the same contribution to In view of equations 
7.9-7.12 it follows that is also composed of N equal terms. 
This finds its mathematical expression in equations 7.31 and 
7.32. A more direct proof that all these terms are equal is given 
below after equation 7.30. Each such unit contains the same 
number, say k, of elementary charges which in each case can be 
arranged in a similar way relative to each other, and relative to 
their surroundings. Let denote the displacements 

of the k elementary charges of the jth unit, and let 


xH'~^^x^dx — f 


* dx. 


t Using 
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denote the whole set of displacements of the jth unit. Then the 
set of displacements X of all the charges in the whole spherical 
region is composed of the sets of displacements 0 : 3 , .--j Xp...,x^ 
of all the N units. Correspondingly the volume element dX 
(cf. 7.21) is the product of the volume elements 

dxj — drjj^drj2,---drjf^ 

of all the units 

dX = dxidx2...dxj...dx^. 

This means that in an integration over the displacements of all 
the elementary charges of the spherical region one can start 
integrating over the displacements of the charges in the first 
unit, denoted by dx-^, and so on. Alternatively, the volume 
element dX may be written as 

dX = dX^dxj, (7.22) 

where dXj — dxidx 2 ->.dxj^idxj^i...dxjsj (7.23) 

refers to integration over the Avhole spherical region except the 
jth unit. 

Now let m(Xj) be the electric dipole moment of the jth unit. 
Then in view of 4.7 the moment M(X) of the whole region is com- 
posed as the vector sum of all the 111 ( 0 ::^), 

M(X) = I m(x,). (7.24) 

j = l 

Therefore M^X) - M{X)M(X) = f m(x^)M(X), (7.25) 

Inserting this into 7.12 leads to 

= f J J (7.26) 

The integrations in each term of this sum will now be carried 
out in two steps, as indicated by equation 7.22. In the jth 
term the integration will first be carried out over the whole 
spherical region except the jth unit, and subsequently over this 
unit. For the first step (volume element dXj) the quantity Tn{x^) 
must be treated as a constant because it depends only on the 
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displacements of the jth unit. The term of 7.26 corresponding 
to this unit can be written as 

= J J m(xj)^M{X)eri^<-^yik^'dX^)dxj 

= J ra{Xj)m.*{Xj)p(Xj)dXj, (7.27) 

where 

m*{Xj) = J M(Z)e-f'(-'')/‘-^dZ_,/J erViXMkr (7.28) 
and (cf. 7.10) 

J)(x.^) - J e-^'(''')/*-^’cZZyj e-W^'^')/'^^'(iZ (7.29) 

have been introduced. m*(Xj) represents the average moment 
of the whole sphere for a fixed set Xj of displacements of the jth 
unit leading to a moment m(Xj). p{Xj) is the probability of finding 
the jth unit with this particular set of displacements. Inserting 
from 7.27 into 7.26 yields 

= 2 ( i^{oi^j)Tn*{^j)p(^'j)dXj. (7.30) 

i=i 

Practically all terms of this sum are equal. For m* represents 
the moment of a large spherical region polarized by one (say 
the jth) of its units whose moment is kept at a value m. Accord- 
ing to electrostatics (cf. Appendix A 2.ii) the same moment m* 
is contained in any other sphere enclosing the jth unit (not 
necessarily co-centric), but with no further restriction of its 
position and radius so long as it is sufficiently large to be treated 
macroscopically. The actual value of m* is thus determined by 
short-range interaction and is independent of the position of the 
jth unit as long as its distance from the surface is sufficiently 
large to allow its interaction with the outside to be treated on 
a macroscopic basis. The number of units for which this does not 
hold can be made very small compared with the total number N 
if the spherical region is sufficiently large. 

The equivalence of all N terms in the sum 7.30 means, of 
course, that all units have the same average polarization. This 
fact might indeed have been used to show the equivalence of all 
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the above terms in a much simpler way. The above considera- 
tions give, however, some insight into the forces determining 
m*. It follows from these considerations that a region makes a 
contribution to m* only if the average moment induced in it 
by the jth unit cannot be obtained by treating the jth unit as 
either a point dipole or as a polarized sphere with moment m. 
Thus it is essentially the short-range forces and the deviation of 
the shape of molecules from a sphere which makes m* different 
from m. 

It now follows that 

mm* = J Tti{xy)Ta^{x^)'p{x^)dx^ 

^ ~ J m{x^)m^{xj)'p{x^)dx^ ■= ... (7.31) 

and hence with 7.30, 

= iVmm*. (7.32) 

Inserting this into equation 7.11 leads to our final expression 
for the static dielectric constant, 


where 


36^ ^ttNq mm* 
263+1 “3 kT~ 
N, = N/V 


(7.33) 

(7.34) 


is the number of units per unit volume. By its definition 7.28 
m* represents the average dipole moment of a spherical region 
embedded in its own medium , if one of its units is kept in a given 
configuration leading to a dipole moment m. mm* is the 
average value of the product mm* taking into account all 
possible configurations and weighing them according to the 
probability of finding the unit in such a configuration, m* 
differs from m because of the existence of short-range forces 
or because of the non-spherical shape of molecules. 

Equation 7.33 is perfectly general. It will now be specialized 
by separating the contribution to due to elastic displacements 
of electrons by treating them on a macroscopic basis. This con- 
tribution can be measured with the help of the optical refractive 
index n because at optical frequencies all other contributions 
have ceased to exist in view of the higher inertia of the heavier 
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particles. If is the electric moment due to electronic dis- 
placement by a field E, on the assumption that atomic nuclei 
are not displaced, then 

n^-l = 47 tMJVE. (7.35) 

Since according to Maxwell’s law the dielectric constant at 
high frequencies is equal to the square of the refractive index, 
this relation follows directly from equation 1.9. 

The spherical region now consists of charges which follow 
the laws of statistical mechanics and which are embedded in a 
continuous medium with dielectric constant n^. Thus if is 
the moment due to all other displacements, is the total 

moment of the substance. Hence equation 1.9 becomes 

6,-1 == 4n(Mj,+M,,)IVE, 

or using 7.35 = ^TrMj^jVE. (7.36) 

Instead of calculating the electronic polarization we may 
introduce it through 7.36 as an empirical quantity on the under- 
standing that n is the optical refractive index. Nearly all of the 
above developments remain unaltered if M(Z), m, and m* now 
refer to non-electronic displacements. Alterations have to be 
introduced only at the following two points: 

(i) In equation 7.5 the cavity field G must now be replaced by 
G', the field inside a spherical cavity with dielectric constant 
v? instead of in an empty cavity. Then according to Appendix 
A 2.14, 

i.e. the denominator (2eg+ 1 ) in 7.5 and in the following equations 
must now be replaced by (2eg+n*). 

(ii) In deriving 7.11 from 7.9 equation 7.36 should be used 
instead of 1.9. Hence 7.11 will be replaced by 


, 477 3eg 

~ ZV kT' 

All further developments remain unaltered, leading to 

, „2 3€g 477^^0 mm* 

* 2c.+n2 3 kT 


(7.38) 


(7.39) 
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instead of equation 7.33. It should be realized that in consider- 
ing the interaction energy V{X) required for the calculation of 
m and m* (cf. equations 7.27-7.31) the electronic polarization 
has to be included on a macroscopic basis. 

Finally, it should be noted that the difference between 7.39 
and 7.33 is of rather a trivial nature whenever ^ i.e. when 
the electronic contributions arc small. In this case 




3 


2e,+7^2 2e,+ l 2’ 




and 7.39 becomes 


mm** 


27tNq ^8 > 


(7.40) 

(7.41) 


while 7.33 would have €^—i on the left-hand side, which thus 
differs from 7.41 by the optical contributions n^—l. 


Mixtures 


Some substances contain a number of distinctly different 
groups of charges such as the different types of molecules in a 
mixture, or the positive and negative ions in ionic crystals. In 
these cases a unit cell would often be uncomfortably large or 
have an undesirable shape. Whenever these groups are well 
separated from each other it is possible to derive a formula 
for in which the contributions of these different groups are 
separated. Suppose that there are z different types of these 
groups. Then z different kinds of units will exist, each being 
representative for one type. If the substance contains per c.c. 

such units of the 1st, 2 nd,... 2 th kind, then 
7.24 can be written as 


M(X) 

V 


Ni Ni+Ni 


(7.42) 


For the whole further development each group can be treated 
separately, so that instead of 7,33 the final result is 


-1 = 


3gg 477 Y 
2€^+1 


mgtnfN^, 


(7.43) 


where m? is the average dipole moment of a spherical region of 
the substance embedded in its own medium if one unit of the 
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ivth kind is kept in a configuration corresponding to a moment 
Trig. Tiig m* is the average value of m^. m*. As before, the spheri- 
cal region must be large compared with a region just big enough 
to have the same dielectric properties as a macroscopic specimen. 


Summary 

A short summary of the results of § 7 is given below for the 
benefit of the reader who has not studied it in detail. 

(i) The static dielectric constant according to 7.1 1 and 7.21 
satisfies exactly the relations 


477- 36^ if 2 ^ 47T €g+2 if2^c 

3V2€,+ 1 kT 3F 3 kT ’ 


(7.44) 


where if 2 is the mean square of the spontaneous dipole moment 
of a sufficiently large sphere of dielectric material of volume V 
embedded in its own medium; if vac fhe corresponding quan- 
tity for a sphere in vacuum. 

(ii) If the sphere consists of components all of which on an 
average make the same contribution to the polarization, equa- 
tion 7.44 can be developed further into equation 7.33 which is 
also of a very general nature. 

By treating the electronic contribution on a macroscopic 
basis the less general equation 7.39 follows. 


8. Special cases 

The results of the previous section are contained in formula 
7.33 [also 7.11] for the dielectric constant which is valid in a 
very general way, and in the more specialized equations 7.39 
and 7.43. All these formulae represent in terms of quantities 
mm* or Jf 2 which refer to properties of the material in the ab- 
sence of a field. To calculate these quantities requires a detailed 
knowledge of the structure of the substance in question and of 
the interaction between the particles of which it is composed. 
In general such a calculation cannot be carried out without the 
use of approximations. The importance of the general formulae 
lies in the possibility of deriving from it approximate formulae 
of various types, each with a clearly defined range of validity. 
It also shows that even for a restricted class of substances, such 
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as dipolar liquids, one cannot hope to find a simple formula 
representing in terms of a few parameters (e.g. the Onsager 
formula), which is valid over the whole stable range of these 
substances. Thus it will be found, for instance, that for dipolar 
liquids the Onsager formula should hold asymptotically at high 
temperatures only, though the deviations from it at lower 
temperatures may be very different for different substances. 

As a first step to get acquainted with the handling of the above- 
mentioned formulae one may assume that all short-range forces 
can be neglected. For spherical molecules this should according 
to § 6 lead to either the Clausius-Mossotti or the Onsager formula, 
depending on whether on© considers non-polar or polar mole- 
cules. The calculations proving this are carried out in the 
appendix (A3). They demonstrate, as the main difference 
between these two cases, that for elastic displacement (non- 
polar molecules) the energy of a unit depends on its dipole 
moment m \ hence it is shown that (which for negligible short- 
range interaction is equal to mm*) increases proportional to 
the absolute temperature T , and the dielectric constant becomes 
independent of the temperature. For rigid dipoles, to take the 
other extreme, the moment of a unit is equal to the dipole 
moment p,, and hence = /x^ independent of temperature. 
Polar molecules, of course, show always both dipolar and elastic 
contributions to 

Polar liquids; Kirkwood's formula 

We shall now proceed to consider the general case of polar 
liquids consisting of molecules with an intrinsic dipole moment 
and with a polarizability cx. Liquids of this type have been 
already investigated in § 6 on the assumption that short-range 
interaction can be neglected; this led to the Onsager formula 
6.36. This assumption will no longer be made in the present 
calculation. 

It was found in § 6 that a molecule in the Hquid state has a 
dipole moment which is different from that in the gaseous state 
(fi.^) . For spherical molecules the difference is due to the polariza- 
tion of the molecule by the reaction field of the surroundings. 
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For no other case than for spherical molecules is it possible to 
calculate the ratio of the moments in any simple precise manner. 
In view of this difficulty it seems reasonable to treat the effect 
of the polarization of molecules in a macroscopic way, assuming 
(as in § 6) that the main contributions are due to electron dis- 
placement. The liquid will thus be considered as consisting of a 
continuous medium with dielectric constant [n = optical 
refractive index) containing dipoles with a moment {Ji. In this 
model a spherical molecule in vacuum consists of a sphere of the 
continuous medium containing a dipole pi at its centre. Accord- 
ing to the appendix (A 2.32) the moment of such a molecule is 


— 


3 


spherical molecules. 


( 8 . 1 ) 


For spherical molecules the present model permits us, there- 
fore, to express the moment pt. of the dipoles by the moment pi.y 
of a free molecule. 

On the basis of our model the general formula for the dielectric 
constant is given by equation 7.39. A unit contains just one 
dipole pi and its moment is, therefore, 


m = pi. (8.2) 

Since all contributions to m* are due to dipole orientation we 

introduce pi* by * * /o ov 

^ m* = pi.* (8.3) 


as the average moment due to the dipoles of a spherical region 
if one of its dipoles is kept in a fixed direction. Equation 7.39 
contains the quantity mm*, where the bar indicates averaging 
over all possible values of m. In the present model the only 
variable is the direction of the dipole pi. Since in a liquid all 
dipolar directions are equivalent, pipi* will have the same value 
for aU these directions. Thus 


mm* = pipi* == p^ji*. (8.4) 

Inserting 8.4 into 7.39 leads to the Kirkwood \K4\ formula 
. ^2 _ 47TiV^0Wi* 


ZkT 


( 8 . 6 ) 


In order to make use of this formula we must calculate pi*. 
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To simplify this calculation it seems permissible to assume that 
short-range interaction between nearest neighbours only need 
be considered (cf. first section of § 5). In this case fx* is the 
vector sum of the moment p. of the centtal dipole, kept in a 
fixed direction, and of the average of the sum of the moments of 
the nearest neighbours. Therefore if z is the average number of 
nearest neighbours, 

jjLp* ~ iJi^l-\-zco8y), (8.6) 

where cos y is the average of the cosine of the angle between neigh- 
bouring dipoles. This average value has to be calculated by 
applying equation 7.28 for m* to the present case. Since the 
directions of the dipoles are the only variables, this equation 
reduces to 

cosy = J coBye-^^l^^ da)^doj2 j J (8.7) 

where U is the part of the energy of interaction between neigh- 
bouring molecules in the liquid which depends on the angle 
between their dipoles; this energy may depend on the state of 
other molecules of the liquid, and U is then the energy averaged 
over all states of the other molecules (considering their proba- 
bilities). dcoj and dco^ are the surface elements of the solid 
angles of the directions of the two dipoles. 

The actual value of cos y depends on details of the interaction 
between the two molecules which may be very different from 
the electrostatic interaction between point dipoles as pointed 
out in § 5; repulsion forces, chemical bond, and other types of 
interaction must be considered. It should be noted, however, 
that a large value of U does not necessarily lead to large values 
of cosy; to do so U must also have the correct symmetry. Thus 
if, for instance, U is proportional to an even power of cosy 
e’^uikT jg even function of cosy, and cosy therefore 

vanishes, in contrast to the case in which U is an odd function 
of cos y. This means that an interaction tending to direct dipoles 
with equal probability either parallel or antiparallel does not 
play any role in the determination of (xp.*. Either of the two 
types of interaction leads, however, to a restriction of free 
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rotation of a molecule termed by Debye [D3] as hindered rota- 
tion. Hindered rotation, as we see, does not necessarily have an 
influence on the dielectric constant. 

The energy U need not only depend on the angle between the 
dipoles, but may also be a function of the direction of the vector 
joining the dipoles. Even so, we can always put 

U ~ t/oven+^)dd» (^•^) 

where reversal of the direction of one of the dipoles does not 
alter C^ven but changes the sign of 

Now let us consider that the temperature is sufficiently high 
to make ^ 

for all values of C/odd- Then 


^-VikT ^ ^-UevenlkT 


and hence, using 


t'-Tl?)' 


J COB yc'- dojidco^ 


= 0, 


one finds 


cosy 


£o 

'kT' 


(8.10) 


( 8 . 11 ) 


where 


Uq ~ j dcoj^dco.^! j (8.12) 

since J dco^ dw.^ ^ 0. 

The energy may be positive or negative depending on 
whether the interaction tends to orient neighbouring dipoles 
anti -parallel or parallel. 

Inserting 8.6 into 8.5, we find a specialized form of the Kirk- 
wood formula, 


47riV^ /x^ 


( 1 + 2 ; cosy). 


(8.13) 


" 2€,+n2 ZkT 

In the particular case of spherical molecules /x can be expressed 
by the vacuum moment jjl^ with the help of 8.1, and hence 

for spherical molecules. (8.14) 
This equation differs from the Onsager formula 6.36 by the 
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z COB y term. According to 8.11 this term tends towards zero at 
temperatures at which hT ^ IT/qI. Therefore it is found that 
the Onsager formula should hold in liquids at temperatures for 
which kT is large compared with the directional part of the 

interaction energy. This energy may have different orders of 
magnitude for different liquids; and this may lead to very 
different ranges of validity of the Onsager formula. In some 
liquids this range may fall into the whole liquid range at normal 



Fig. 5. Temperature -dependence of the dielectric constant in the high- 
temperature region. The full line represents Onsager’s formula. Short-range 
interaction tending to orient dipoles parallel (cos y > 0) leads to larger values 
of ( ) ; the opposite case (cos y < 0) loads to smaller values ( ). 

pressures, whereas for others zcosy may be appreciable in the 
normal range of the liquid and may have either sign (cf. Fig. 5) 
thus making the Onsager formula inapplicable in these cases. 
In gases, in view of the large distance between molecules, 
dipolar forces only need be considered. For very small densities 
interaction may be completely neglected, leading to equation 
6.15 for 6g. For higher densities Onsager’s formula should hold. 
This has been shown more directly by van Vleck [V2\ 

Dipolar solids 

The average potential energy of a polar molecule in a crystal- 
line solid in general depends on the direction of its dipole relative 
to the crystal axes. In liquids in contrast, though the dipoles of 
neighbouring molecules have a tendency to orient themselves 
in definite directions relative to each other, the average energy of 
a single dipole is the same in all directions because a liquid does 
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not possess any preferential direction. The average potential 
energy of a molecule in a solid is often described as being due 
to a crystalline field acting on the dipoles. f The crystalline 
field, being due to the interaction between molecules, normally 
depends on temperature. Usually there are several dipole direc- 
tions for which the average energy of a molecule has a relative 
minimum; they will be called equilibrium directions. The poten- 
tial barrier between these equilibrium directions is in general 
very high and prevents free rotation of the molecules even at 
temperatures near the melting-point. The hypothesis that above 
a critical temperature molecules should be able to rotate freely 
was introduced by Pauling [Pi] to explain the sudden change 
of the dielectric properties of many solids at a critical tempera- 
ture. This idea was a useful working hypothesis, but the present 
evidence seems to indicate that rotation does not occur in most 
solids (e.g. LI). Instead, the transition has to be considered as 
leading from an ordered arrangement of dipoles to disorder. 

This transition from order into disorder will be used as a guid- 
ing principle in discussing the general behaviour of the dielectric 
constant of polar solids. Its detailed properties, of course, 
depend on details of the crystalline field, but the main features 
can be discussed without such a detailed knowledge. J To de- 
scribe the main features of an order-disorder transition consider, 
a simple two-dimensional model consisting of dipoles arranged 
in a cubic face centred lattice. The crystalline field is assumed 
to lead for each molecule to two equilibrium positions with oppo- 
site dipole direction. At the absolute zero of temperature the 
solid will be in its lowest energy state. In this state the 
dipoles are arranged in an ordered way, but there are a number 
of possibilities for doing this. The energy of these ordered states 
depends on the particular interaction between the molecules 
which, as has been pointed out before, contains interactions 
of various types (e.g. dipolar, repulsive, etc.). These ordered 
arrangements can be classified according to whether or not they 

t The effect of various types of crystalline fields on the average dipole 
direction was considered by Bauer [BI], Frank [F5J, and others. 

J Various calculations of the behaviour of the dielectric constant in order- 
disorder transitions have been carried out for special models [K5t FS]. 
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lead to a residual dipole moment of the crystal. Thus, if all 
dipoles are parallel (Fig. 6 a), the solid has a permanent dipole 
moment, but if, for instance, dipoles in the corners have the 
opposite direction from those at the centres (Fig. Gb) the dipole 
moment of the crystal vanishes. In the former case the energy 



Fig. 6. Equilibrium positions of dipoles for a simple model of a crystalline 
solid. Ordered state : (a) in the case of permanent polarization ; (6) for vanishing 
polarization (two dipoles per miit coll); (c) disordered state, both directions 
have equal probabilities. 




Fig. 7. Average potential energy of a dipole (due to short-range interaction) 
as a function of its direction for the model of Fig. 6, assuming that all the 
other dipoles stay in the fixed directions indicated in Fig. 6: (a) for the central 
dipole of Fig. 6 6; ()3) for a corner dipole of Fig. 6 6 or for any dipole of Fig. 6 a ; 
(y) for the disordered case, Fig. 6 c. 


of a dipole is lower when it is directed to the right than for the 
opposite directions (cf. Fig. 7 J3). In the latter case the energy of 
the corner ‘dipoles behaves in a similar way, but for the centre - 
dipoles the two directions are exchanged (Fig. 7 C3^). The internal 
fields for these two types of dipoles are, therefore, different; 
they transform into each other by a reflection on a plane perpen- 
dicular to the dipole direction. In the former case (a) all lattice- 
points are equivalent and a unit cell contains one single dipole. 
In the latter case (6) two different types of sites, comers and 
centres, have to be distinguished, and a unit cell contains two 
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dipoles. Alternatively, the crystal may be considered as com- 
posed of two lattices of type (a), but with opposite directions of 
the dipoles. 

Clearly the division of ordered states according to whether 
or not they lead to a residual moment holds quite generally 
(i.e. for three-dimensional lattices as well), although there may 
be structures leading to more than two equilibrium directions 
for each dipole. In the following the case of two equilibrium 
directions only will be considered; most results hold qualita- 
tively, however, for more complicated structures as well. 

Now, starting from an ordered state at the absolute zero of 
temperature, imagine the temperature to increase gradually. 
This will cause some dipoles to turn into the second equilibrium 
position in which they have a higher energy, thus creating some 
disorder. This in turn will cause a decrease of the average energy 
difference between the two directions because the energy of 
interaction has its lowest value in the completely ordered state. 
The average energy difference of a molecule in the two equi- 
librium directions is thus a function of temperature, say 

V{T) ^ 0, 

which decreases with increasing temperature. Jt is found that a 
critical temperature Tq exists at which it vanishes (Fig. 7 y). 

Let us now define the direction which a dipole at a given 
lattice-point has in the completely ordered state (at T = o) 
as the 'right’ direction, and the opposite direction as the 'wrong’ 
direction. Thus in Fig. 6 a, ->is the right direction; in Fig. 
6 6, is the right direction for corner-dipoles and is the right 
direction for centre-dipoles. Let w be the probability of finding 
a dipole in the wrong direction, and (1— therefore, the pro- 
bability of finding it in the right direction. Since V{T) is the 
energy difference between the two positions, it follows from 
statistical mechanics that 

^ (8.16) 
1—w 

and hence 

^-nT)ikT I 

'ijij =— 1 <111 — . 


(8.16) 
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A calculation of the temperature dependence of F(r), and hence 
of w is very difficult, but approximations have been devised to 
simplify it. A summary of these methods can be found in 
reference [Ai]. We shall not give details of such calculations 



Fig. 8. Temperature -dependence of the average energy difference V{T) 
between opposite dipole directions, schematically. 



Fig. 9. Temperature -dependence of the probability w{T) of finding the dipole 
in the ‘wrong’ direction, schematically. 

here. It is sufficient for us to note (cf. Fig. 8) that near T = 0, 
V{T) is nearly constant, but as T approaches it decreases 
rapidly towards zero, which it reaches at T = JJj. Thus 

V{T) = 0 if T ^ Ji- 
lt follows (cf. Fig. 9) that for low temperatures 

w <^l iikT (8.17) 

whereas w ~ \ if T ^ Tq, (8.18) 

Thus above Tq both dipole directions have equal probabilities, 
i.e. the lattice is disordered (cf. Fig. 6 c). Calculations also show 
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that Tq is of the order of V(^)jk] the exact value depends on 
details of the structure. 

The quantity ( 1 — 2w) is a measure for the degree of order of 
the lattice because \~-2w ~ \ at T ™ 0, and l~2r(; = 0 if 
T > Tq. This type of order is often called long-distance order 
because it defines right and wrong directions for any lattice- 
point. In contrast, short -distance order is the order of neigh- 
bours relative to each other. It means that, in view of the 
interaction, the direction of a dipole is always influenced by the 
directions of its neighbours. Each dipole tends to orient itself 
in a certain direction relative to its neighbours. Long-distance 
order vanishes in the disordered state, or in liquids. Short - 
distance order persists, however, though it decreases with 
increasing temperature. In fact the absolute value of cosy 
introduced above for liquids is a measure for the short-distance 
order. 

Let us now discuss the dielectric properties of ordered solids 
near T ^ 0. If the order is of type (a), the solid is permanently 
polarized. This case will not be further considered here. For the 
other kind of order the only contribution of the dipoles to the 
dielectric constant is due to elastic displacement of the dipole 
direction by an external field. Thus near T = 0 the dielectric 
constant is larger than say where n is the optical 

refractive index due to elastic displacement of electrons (cf. §§ 6 
and 7). The difference is usually small. 

At higher temperatures, where dipoles are capable of changing 
their directions, €g increases with temperature, as will be shown 
presently. As in the case of liquids, contributions due to elastic 
displacements will be treated in a macroscopic way, and the 
following considerations are, therefore, based on equation 7.39. 
For the type of order which we consider at present a unit cell 
contains two dipoles as mentioned before, and at T = 0 these 
dipoles have opposite direction. For the two-dimensional case 
discussed above the unit cell contains a comer- and a centre- 
dipole as shown in Fig. 6 6. The moment m of the unit cell 
depends on the directions of the dipoles. Since the potential 
barriers between the equilibrium directions are very high, it 
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seems sufficient to consider equilibrium directions only. The 
moment m can thus have a discrete number of values 
nig,.- - If Pi is the probability that the unit cell has the moment 
equation 7.31 applied to the present case becomes 

mm* = 2 m;m? Pi, (8.19) 

where the sum extends over all values of m^. mf is the average 
moment of a sufficiently large spherical region around the unit 
cell if the moment of the latter is kept at nij-.. 

In our particular case of two equilibrium positions with 
opposite dipole direction the following four states exist: 


i 

Configuration 


Pi 

energy 

1 

< — > 

Wi — 0 

Pi = (1-U>)» 

0 

2 


~ — 2/a 

P 2 = w{l-~w) 

V{T) 

3 

-> — > 

m3 = 2/x 

Pq ^ w{l-w) 

V{T) 

4 


m4 = 0 

Pi = 

2ViT) 


Here the probability p^ was assumed to be the product of the 
probabilities of finding either of the two dipoles in its respective 
direction. Thus for i — 1 both dipoles are in ‘right’ positions, 
each having the probability (l—w). This way of calculating 
is not free from objections because it does not take into account 
the correlation between neighbouring dipoles giving rise to 
short -distance order. It is, however, a sufficiently good approxi- 
mation for our present qualitative treatment. 

With the above four states equation 8.19 becomes 


mm* = 8w(l— w?)pp.*, (8.20) 


where + 2^* is the moment of a spherical region around the unit 
cell if the moment of the latter is kept at di2|X. 

Inserting from 8.20 into 7.39 [it should be remembered that 
this equation holds for poly-crystalline material], we find by 
adding (see above) as contribution of elastic dipole- 

displacement 




_ 3€^ 47rJ\rQpp.* 

2€,+n^ UT 




( 8 . 21 ) 


Here is the number of dipoles per unit volume (not the number 
of unit cells as in 7.39). 
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For temperatures below the transition point, w{l—w) is 
probably the determining factor leading to an increase of 
with temperature. Above the transition point w ~ i.e. 
4^iv{l—w) “ 1. Equation 8.21 then becomes nearly identical 
with the Kirkwood equation 8.5 for liquids. A difference remains 
(i) because 8,21 contains too instead of at the left-hand side, 
and (ii) because /x* is defined differently. f 

A more exact calculation can be carried out, however, above 



Fig. 10. Temperature-depondenoe of the dielectric constant of a dipolar solid 
showing an order-disorder transition, schematically. 

jTq since the unit cell then contains only one dipole. This calcula- 
tion follows the lines given above for liquids and leads exactly 
to the Kirkwood formula 8.5. The difference from 8.21 should 
thus be due to the approximations made in its derivation. 

Thus the dielectric constant of a solid (cf. Fig. 10) which has no 
permanent polarization should start at T = 0 with a value 
which is slightly larger than and should rise very slowly at 
first and more rapidly as T approaches Tq, provided is below 
the melting-point. Above Tq, should decrease with T and 
there should be no appreciable variation at the melting-point. 
If, on the other hand, Tq cannot be reached below the melting- 
point, then €g increases up to the melting-point, and then de- 
creases. 

The similarity of the dielectric behaviour of disordered solids 

t Kirkwood: is the moment of a spherical region if one of its dipoles is 

kept in a fixed direction p//x. Equation 8.21 : 2ji* is the corresponding moment 
if two dipoles (of one unit cell) are kept parallel with a total moment 2}i. 
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and of liquids is not surprising. In fact we have seen before 
(§§ 4 and 6) that a dipole having two equilibrium positions with 
equal energy behaves similarly to a dipole with a continuous 
range of such positions. 

Finally, it should be remarked that as for liquids the quantity 
IfJL*— fji|//x is a measure for the short-distance order. Above the 
transition temperature Tq this quantity can be found experimen- 
tally if is known over a sufficiently large temperature range. 
Near it may be expected to depend sensitively on the structure. 

Intramolecular rotation 

Consider briefly the case of a liquid consisting of molecules 
in which intramolecular rotation is of importance. Again the 
unit cell contains one single molecule, but its dipole moment 
li(x) may depend on the set of coordinates x describing intra- 
molecular rotation. Thus we put m(a:) = fx(:r), and introduce 
correspondingly m*(x) = as moment of spherical region. 
In the Kirkwood formula 8.5 the term (ijx* must then be replaced 
by the average jjl(x)(x*(x). If in particular the molecule has a 
number of semi-stable states with moments (Xj, {X 25 ***> besides the 
ground state whose moment is jjLq, then similar to 8.19 

jA(;c)|A(a;)* = (8.22) 

i>0 

where represents the probability of flnding the molecule in 
the state i. 



CHAPTER III 


DYNAMIC PROPERTIES 

9. The establishment of equilibrium 

The present chapter will be devoted to an investigation of the 
dynamic properties of dielectrics. This will be found to be a 
much more difficult task than the development of the theory 
of the static properties, and in fact it will be seen that only for 
dilute solutions of dipolar molecules has it been possible so far to 
carry out quantitative calculations. The reason for the greater 
difficulties in dealing with the dynamic properties is evident if 
we remember that for the static case it was not necessary to 
investigate the kinetic properties of molecules. In the present 
section a qualitative discussion of the problems involved will be 
given. 

It will be remembered that in § 2 it was found that in alternat- 
ing electric fields the electric displacement shows a phase-shift 
with respect to the electric field. This leads to the introduction 
of two dielectric constants and both of which depend on the 
frequency co/27r and which according to 2.8 can be considered as 
real and imaginary components of a complex dielectric con- 
stant €(ca). The two quantities and eg are not independent of 
each other, but can both be derived from a function oc{t) of time 
with the help of equations 2.14 and 2.15. Making use of 2.8, 
these equations can be written as a single complex equation, 

C30 

e{o}) = j dx, (9.1) 

0 

where a: is a variable of integration. 

The function a(^), according to § 2, describes the decay of the 
polarization of a dielectric with time if the external field is 
suddenly removed. Alternatively the gradual increase of polari- 
zation with time to its equilibrium value if the dielectric is 
brought into a constant field can also be described with the help 
of the decay function a(<). It is thus seen that with the use of 
the macroscopic equation 9.1 the frequency dependence of the 
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complex dielectric constant (and hence of dielectric loss, cf. 3.15) 
is uniquely connected with the way in which equilibrium is estab- 
lished in a dielectric brought into a constant field. We shall 
choose this latter process for the qualitative considerations 
which follow because it is very suitable for the purpose of showing 
in a simple way the difficulties involved. 

As in § 4, consider the two characteristic types of displace- 
ment of charges, namely (i) elastic displacement, and (ii) dis- 
placement to another equilibrium position, and assume that all 
interaction forces between particles can be neglected. From § 4 
we then know the average position of a charge in the presence 
and in the absence of a field, or rather its average displacement 
by the field. It is now our task to consider the influence of the 
field in more detail by assuming that the dielectric consists of 
such an assembly of (non-interacting) charges which, in the 
absence of a field, are in thermal equilibrium. 


Case (i). Elastic displacement 

Each particle of mass m and charge e is bound elastically to its 
equilibrium position and will carry out harmonic oscillations 
with frequency coq about it. If r is the displacement, then in the 
absence of a field 

rt 

(9.2) 


dt^ 


-6ugr, 


and hence r Cocos(a>o^+8o)» {^•^) 

where the maximum amplitude Cq and the phase Sq are indepen- 
dent of time, and the energy is given by ^mcol (7§. If the motion 
is not perturbed, the oscillator will keep this same energy in- 
definitely. It should be realized that this is in contradiction to 
the postulates of statistical mechanics, according to which 
within a sufficiently long interval of time the oscillator should 
be found with various energies according to the Boltzmann 
theorem, and its average energy at a temperature T should be 
kT. To achieve this some kind of interaction with other particles 
or with the surrounding medium — permitting an exchange of 
energy — must be assumed. Frequently it is assumed that this 
interaction takes the form of collisions of extremely short dura- 
tion. This means that the equation of motion 9.2 holds except 
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during a collision. Therefore the solution 9.3 holds between two 
collisions, but at a collision both the amplitude Cq and the phase 
Sq usually change their value. The assumption of collisions of 
very short duration means that each oscillator actually satisfies 
the equation of motion obtained by neglecting all interaction 
terms, but as a result of a collision it makes a transition into a 
state with different energy (and phase). This, as may be shown 
by statistical mechanics, leads to the correct average energy if 
averaged over a time which is long compared with the time 
between two colJisions. It is of great importance that in equili- 
brium this average is independent of the nature of the collisions. 
Equilibrium properties can, therefore, be derived without con- 
sideration of collisions — quite in contrast to the question of 
attaining equilibrium which we intend to investigate presently. 
It should be remarked first that the collisions need not occur 
between the oscillators themselves as one would expect in a gas 
of such oscillators. They might be considered as dissolved in 
another medium (liquid or solid) which is in thermal equilibrium, 
and the collisions to be considered occur between the oscillators 
and the particles of the medium. 

Assume now a constant field f to be applied at the time t = 0. 
Then the displacements satisfy the new equation of motion 4.8 
instead of 9.2. A particle now oscillates about a new position of 
equilibrium displaced by a vector r relative to the former one, 
but the motion is still harmonic (cf. 4.9 and 4.10). Both the 
maximum amplitude C and the phase 8 usually change, however, 
at ^ ~ 0, although r remains continuous, as is shown in Eig. 11. 
This assumes, of course, that no collisions occur during the 
represented time interval. 

To obtain the polarization P of a dielectric substance con- 
sisting of a laige number of oscillators the vector sum of all 
displacements r has to be formed according to 4.6. At a given 
instant of time ^ < 0 in the absence of a field the phase-angle 8 q 
can be assumed to have any value between 0 and 27r with equal 
probabilities. The vector sum of all displacements, and hence 
the total electric moment of the dielectric, will, therefore, vanish 
at any instant of time. Now let us assume as a first approach 
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that no collisions occur after the field has been applied at if = 0. 
Then since all individual displacements vary continuously the 
same must hold of the polarization P. On the other hand, the 



(a) 



Fig. 11 . Time dependence of the displacement r of a harmonic oscillator if a 
constant field / is applied at the time ^ — 0. The displacement and its time 
derivative remain continuous at < = 0, but the average displacement is 
altered. Two characteristic cases are shown: (a) both phases before and after 
application of the field vanish ; the maximum amplitude C for < > 0 is smaller 
than it was for t < 0; {b) both phases are equal to — tt and the maximum 
amplitude is larger for ^ > 0 than for ^ < 0. 


time average of the displacement of any single oscillator is equal 
to r so that the average polarization too is different from zero. 
It will therefore oscillate about its average value with the 
frequency cuo/27r. This means that after application of the field 
the oscillators are no longer distributed uniformly over all values 
of the phase 8, as can also be shown by direct calculation. 

Collisions, as has been pointed out above, alter the phases. 

4980.11 P 
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They will tend to suppress the oscillations of the polarization 
about its equilibrium value. One should expect, therefore, that 
after application of a field the equilibrium polarization is reached 
steadily if the average time between two collisions of an oscillator 





(c) 

Fig. 12. Time dependence of the polarization of a dielectric if a field is applied 
at the time f ~ 0. (a) and (6), the material consists of dipolar oscillators, 
(a) without, (6) with collisions, leading to oscillations of the polarization about 
its average, (c) The material consists of rigid dipoles with several equilibrium 

positions. 


is shorter than its period. Otherwise transient oscillations of 
the polarization will be excited which last for a time of the order 
of the time between two collisions (cf. Fig. 12 a, 6, and c). 

Ca^e (ii). Displacement into another equilibrium position 
As in § 4 (cf. Fig. 4), a charged particle is assumed to possess 
two equilibrium positions A and B at distance b from each other 
and separated by a potential barrier. In the absence of an 
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external field the potential energy of the particle is the same in 
A and B. In the presence of a field f according to 4. 12 there is a 
difference in the potential energies at A and B. In an assembly 
of a great number of particles of this kind in the absence of a 
field an equal number will oscillate about each of the two 
equilibrium positions with an average energy k7\ if thermal 
equilibrium is assumed to exist, and if the height of the potential 
barrier is large so that the potential energy H of a particle on top 
of it is large compared with hT, 

H:p>kT. (9.4) 

This means that the fraction of particles with sufficient energy 
to go over the top of the potential barrier is extremely small. 
Its order of magnitude is given by the Boltzmann factor 

/kT). 

Assume again that at the time / — 0 an electric field f is 
applied. Then in the absence of collisions the number of particles 
oscillating about A and about B is not altered, because the only 
action of the field is to alter slightly the equilibrium position — 
as pointed out above, case (i). The field is not able, however, to 
lift a particle over the potential barrier (cf. § 4). On the other 
hand, in equilibrium the number of particles near A is larger by 
a fraction of the order ebfjkT (independent of H) than the num- 
ber of particles near B, as was shown in § 4. In the absence of 
collisions, therefore, equilibrium cannot be established. This 
was also found to be so in case (i), but in the present case lack of 
equilibrium has a much more serious effect. For in adding all 
the displacements in order to calculate the polarization, lack 
of collisions meant in case (i) that the polarization was 
oscillating about its equilibrium value. In the present case, 
however, the required displacements cannot be carried out 
without collisions, which means complete lack of polarization of 
the required type. 

Collisions will tend to establish equilibrium. The time re- 
quired for this process depends on details of the model. It is often 
assumed that collisions on either side of the potential barrier 
are very frequent, so that the particles on each side may always 
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be considered to be in equilibrium amongst themselves. By this it 
is meant tliat a particle having sufficient energy to flow over the 
potential barrier to the other side will suffer so many collisions 
that it is very unlikely to oscillate back to the side from which 
it came originally because on an average over many colhsions 
a ])article has an energy kT H, 

Starting now from equilibrium in the absence of a field, the 
immediate effect of the application of a field in the A-B direc- 
tion is to lift the potential near A by the amount efb. As a con- 
sequence the fraction of particles with sufficient energy to move 
over tlie potential barrier is approximately given by 

for tlie A B and the B A directions respectively,']' because 
tlie height of the potential barrier measured from A is now only 
H~ehi. Thus if is the frequency of oscillation of a particle, 
the probability per second for the transfer of a particle from B 
to A is given by 

«’2i = (9.5) 


and for the transition A B considering 4.18 it is 
.... . ^ 

(9.6) 

Thus if at any instant of time there is a number of particles N^{t) 
at A and a number K^{t) at B, a number "'^dll flow per 

second from Ato B and a number from JB to ^ . Therefore 
the rates of change of N-^ and respectively are given by 

(9.7) 

(9.8) 

t By introducing a different normalization of the potential of the external 
field the two exponents can be written in a more symmetrical way as H — Jc(bf ) 
and leaving the results unaltered. 
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It follows that the total number N of particles, 

N = 

is independent of time, as required, because 
dt dt dt 

Subtracting 9.7 from 9.8 with the use of 9.9 yields 

j^{N^—Ni)= —(Wi^+W2i)iN^—N^)+(Wi2-Wzj)N. (9.11) 


Now according to 9.6 and 4.18, 



(9.12) 

efb 

W12—W21 — 

(9.13) 

Using these two equations, 9.11 becomes 



(9.14) 

If we assume 


II 

o' 

II 


at t = 0, 9.14 is solved by 


N efh 

(9.15) 


The induced polarization is proportional to and thus 

approaches its equilibrium value exponentially, as shown in 
Fig. 12 c. 

In contrast to case (i), equilibrium is reached in the present 
case by the influence of the field on the probability of transfer 
of a particle from A to B and J5 to ^ respectively, according to 
equations 9.5 and 9.6 or 9.12 and 9.13. This model can be 
generalized by introducing transition probabilities Wi 2 and W 21 
measuring the probabilities per second for the transition of a 
particle from A to B and from B to A respectively, without 
specifying them by equations 9.5 and 9.6. This means that 
equation 9.11 holds as before. Instead of calculating W 12 and 


(9.9) 

(9.10) 
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W 21 from a model, and tlien deriving 9.12 and 9.13, these equa- 
tions can also be obtained from 9.11 by considering solutions 
corresponding to equilibrium which can be obtained from § 4 
without the knowledge of the transition probabilities. Thus 
since in equilibrium d{Ny^~N^jdt — 0, in the absence of a field 
™ iVg, and hence from 9.11 In the presence of a 

field, using 4.13, 4.12, and 4.18, 

N^-N, = N(p^-p^) = f (9.16) 

This must be a solution of 9.11 if N^~N^ is independent of time, 
which leads to 9.12 and 9.13. 

The model discussed above can be furtlier generalized by 
considering dipoles which can have several equilibrium positions 
with certain transition probabilities between them. In equi- 
librium in the absence of a field the number of transitions from 
a given position is just balanced by the number of transitions 
from other positions to it. An electric field alters the transition 
probabilities, and hence, in equilibrium, the distribution of 
dipoles over the various positions. The time required to attain 
equilibrium depends on the transition probabilities. 

The preceding considerations have thus shown that in the case 
of elastic binding the field displaces the charges which then 
oscillate about their equilibrium positions. In the case in which 
charged particles or dipoles possess several equilibrium positions 
the field does not, by immediate action on the charges, transfer 
themi to their new positions, but it alters the transition pro- 
babilities between them. This in turn leads to the establishment 
of equilibrium. 

10. The Debye equations 

In the present section it is intended to derive equations for the 
frequency dependence of the complex dielectric constant €(oj) 
which should hold for the case of dilute solutions of dipoles in 
liquids and solids and for a few other cases. These equations 
were first established by Debye [D2] and subsequently have 
been applied to many substances, not always, unfortunately, 
with the necessary discrimination with respect to the intended 
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rangt'. of validity. We shall base the considerations of this 
section on the hypothesis that in constant external fields equi- 
librium is attained exponentially with time, as corresponds to 
case (ii) of § 9. We thus assume for the decay function a(^), 

(x(t)oc (10.1) 

where r is independent of time but may depend on temperature. 
That this assumption leads to the required properties can easily 
be shown from the relationship 2.12 between the electric field 
E(t) and the electric displacement both of which may de- 
pend on time. In 2.12 it was assumed that both JJ and E vanish 
for times t < 0. If this is no longer tlie case, then 2.12 has to be 
replaced by 

D(t) = j E{u)<x{t-u)du. (10.2) 

— 00 

This integral equation can easily be transformed into a differen- 
tial equation. For on differentiating 10.2 with respect to the 
time and making use of 

--«(<). (J0.3) 

at T 

which follows from 10.1, we find after multiplication by r 

t 

T^-^ = €*T^:5^ + Ta(0)i5{<)- f E{u)cc(t-u} du. 
at at j 

— 00 

(10.4) 

Adding 10.2 and 10.4 yields 

r~{D-^^E) + {D-^^E) = raH))E. (10.5) 

To determine the constant a(0) consider the special case of equi- 
librium in a constant field. This means that 

|(Z>-e„^) = 0, D = ,,E, 

and hence from 10.6 

Ta:(0) = €^—€ 00 * (10.6) 
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Therefore inserting 10.6 into 10.5, 
d 


JS?) + (D-eooi?) - (es-~e^)E 


(10.7) 


is found as the differential equation connecting D(t) with E(t) 
on the assumption that the decay function «(<) is given by (cf. 
10.1 and 10.6) 

( 10 . 8 ) 


a(<) = 


We shall now use equation 10.7 in the investigation of the 
approach to equilibrium of a condenser. The following two 
cases are to be considered: 

(a) Constant charge on the condenser plates. Then 
dD ^ ^ j. 


dt 


= 0 , D 


and hence, using 10.7, 
dt 


= i.e. Zj,-e,i?oc 


where T'==~r. 

(6) Constant voltage at the condenser plates, i.e. 

d^ 

dt 

It follows with 10.7 that 

i.e. D~eJJ,ace-<lr. 


(10.9) 

( 10 . 10 ) 


0 , 




( 10 . 11 ) 


Both cases thus lead to exponential approach to equilibrium. 

In periodic fields assume E to be represented by equation 
2.9, i.e. E oc ex'p{—icut). Then, introducing the complex di- 
electric constant e (cf. 2.8), we find, using 2.10, 

~ = -icoE, D = €{a))E, ^ = -i<o€(co)E. 

at dt 

( 10 . 12 ) 

Introducing this into 10.7 leads to 

e(a>)-e„ = i2i;^. (10.13) 

1 — ^caT 
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An alternative way of deriving this equation is to insert 10.8 
into 9.1, 

00 

€(a>)— = (eg— Coo)^ J dx, (10.14) 
0 

which after integration results in equation 10.13. 


Separating the real and imaginary parts in 10.13 according to 
2.8 we find 



(10.15) 

, (f8-eoo)‘t>T 

(10.16) 

and for the loss-angle (j>, using 2.5, 


^1 €^+€^aj^T^ 

(10.17) 


Equations 10.15-10.17 (also 10.13) will be denoted as the 
Debye equations and the constant r will be called the relaxation 
time. They describe the properties of a dielectric substance in 
alternating fields on the assumption of an exponential decay 
function 0 L{t) (cf. 10.8). Some models leading to such a decay 
function will be studied in § 11. Most of them require 

< 1, (10.18) 
a condition which normally is fulfilled in dilute solutions only. 

In discussing the properties of the Debye equations it should 
be noted that the dielectric constants cj and depend on at least 
two parameters, the angular frequency w and the temperature 
7\ The frequency dependence is expressed explicitly, but the 
temperature appears implicitly through and r, both of 

which usually depend on T, They may depend on other para- 
meters as well, variations of which will not be considered here. 
For the following it will be assumed that both and are known 
as functions of T. If r were known as well, a new variable 

z = logWT = logCU+logT (10.19) 

could be introduced in terms of which 10.15 and 10,16 become 


^00 I ^ ^ 

^8 — ^00 1 + e ^^ 6^+6“^* 


€g~€ao 


( 10 . 20 ) 
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Fig. 13 shows these functions. It will be noted that 
is a symmetrical function of z. 

Actually r cannot be assumed to bo known but must be 
determined from measurements of and €3 at various frequencies 
and temperatures. Assuming that the Debye equations are 
fulfilled t{T) can, however, easily bo found from the frequency 



Fig. 13. The Debye functions € 2 /(€^ — Coo) (tull line), and (ei — €oo)/(eg~eoo) 
(dotted line), according to equations 10.20. 

at which €3 maximum. In fact at constant temperature 

the angular frequency of this maximum is determined by 

— = 0 , if o) — and T — constant. ( 10 . 21 ) 

Hence using 10.16, 

( 10 . 22 ) 

Inserting this value into 10.15-10.17 we find for the dielectric 
constants and the loss angle 

«i = ^2 == f=o). tan^ = 

^ 8+^00 

if a> (10.23) 

Thus from the frequency at which eg has its maximum we find 
T, and from its height we ought to obtain e^o if the Debye 
equations are fulfilled. 

Instead of considering the maximum of eg it is often preferred 
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to determine the angular frequency co^ at which the loss-angle 
<f> has its maximum. We then require 

= 0 , if w ^ OJA and T ^ constant. (10.24) 
doj ^ 


Using 10.17, this yields 



(10.25) 


In view of 10.18, oj^ and oj^ will be nearly equal for most sub- 
stances for which the Debye equations can be expected to hold. 
Inserting 10.25 into 10.15-10.17 we find 






^{^8 ^00 )> 


tan<^= /-; if a> = co^. (10.26) 

It is an interesting feature of equations 10.23 and 10.26 
that the values of and Cg at the fre(juency oj„^ or co^ at which 
either eg or tan (f> has its maximum are independent of this fre- 
quency and of the relaxation time, and are expressible by the 
static dielectric constant and the high-frequency dielectric 
constant e^o. 

We shall use this occasion to give a more precise definition of 
than that given previously. According to the Debye equa-- 
tions decreases from €, to eoo in the frequency region in which 
€2 has relatively large values. Therefore 600 is the value which is 
asymptotically approached by 6 i(a>) at frequencies sufficiently 
larger than to make € 2 (oj) relatively small. This value 
of €00 need not coincide with the optical dielectric constant 
because most substances absorb in the infra-red region (see 
also § 13). 

If £00 is very small, as is required by 10.18, then e^, Coo, and 
€i are nearly equal, 

£, - £^ - ,,{ojh (10.27) 

and very precise measurements would be required to obtain 
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€'oo- This quantity can, however, be obtained with the help 
of the phenomenological relation 2.18 if 

€^(co) == cj^tan^ e^taiie^ 

is known for all frequencies for which it has appreciable values. 
For on inserting eg from 10.16, equation 2.18 becomes 


, .doj , r 

eglo;) — =- (e^— eoo)- - 

CO TT J i 

0 


cor dio 
-{-CO^T^ CO ’ 


(10.28) 


which is an identity because 


00 


J 


cor 

l+co^r^ 


dco 

CO 


TT 

2 ’ 


In the particular case of dilute solutions of dipolar molecules 
in non -polar liquids e^—eoo can be inserted from 6.24. In this 
case the dielectric constant cq of the solvent is nearly equal to 
Coo, i.e. with 10.27, 

€o e® ^ Cs- (10.29) 

Thus from 6.24, 10.16, 10.17, and 10.29 



SfcT \ 3 ; I (2e,+»2)(e,+ 2) 


(or 

l-foiV’ 


(10.30) 

for spherical molecules and 1. 


Here n is the refractive index of a pure liquid of the dissolved 
molecules. The factor containing n is not of great importance 
because { usually differs from unity by a few per cent, 

only. 

In the above discussion of the Debye equations we have been 
mainly concerned with the frequency dependence of the dielec- 
tric constants using the temperature as a parameter. In practice 
€i and eg are often measured as functions of the temperature T 
with the frequency as a parameter. This does not permit such an 
immediate comparison with the Debye equations because the 
latter do not depend on T explicitly. Assuming that e,— 
is known as a function of T, it is advisable in this case to plot the 
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function € 2 /( 6 ,— Coo) against T, Its maximum is determined by 
(cf. 10.16) 

n ^ / ^2 \ ^ ^ I • 1 

^ ^ T+^ ^ 

as before 10 . 22 . This means tliat the maximum for a given 
frequency occurs at a temperature at which 

r{TJ = l/co. (10.31) 

In the case that and €2 are measured as functions of cd for a 
given temperature T, the Debye equations can be checked 
immediately once t(T) has been found from the position of the 
maximum of according to 10 . 22 . For r does not depend on 
cDj and its value can thus be insei*ted into equations 10.15-10.17 
w^liich then contain no unknown quantities. If, on the othei* 
hand, and measured for a given frequency as functions 
of T, the value found for r(TjJ must not be introduced into 
10.15-10.17 because r varies with T, Now if 10.18 is fulfilled 
and thus an exact measurement of €g—€ao is difficult, knowledge 
of the temperature dependence of eg at a given frequency does 
not permit a check of the Debye equations. If, however, e^— Coo 
and € 1—600 are known as functions of the temperature, ojt can 
be derived from 10.15, i.e. 

(a,T)2 = (10.32) 

and then inserted into 10.16, resulting in 

€2 = V{(e,— €,)(€, —e„)}. (10.33) 

Equations 10.32 and 10.33 are equivalent to the Debye equations 
10.15 and 10.16 in their original form. They have the advantage 
of showing clearly that the relaxation time r is a quantity which 
can be calculated from the measurable quantities according to 
10.32, whereas 10.33 is a relation between measurable quantities 
only. Their disadvantage is that they require measurements of 
€ 5—61 and € 1—600 which in the cases where the Debye equations 
are supposed to hold are very small (cf. 10.18 and 10.27). When- 
ever this is possible, however, 10.33 affords a very simple test 
of the Debye equations whether temperature, frequency, or both 
are varied. 
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11. Models for the Debye equations 

In the previous section it was shown that an exponential 
decay function leads to the Debye equations. In the present 
section, models for which these equations hold will be studied 
and expressions for the relaxation time will be obtained. A very 
simple model has already been discussed in case (ii) of § 9 and 
§ 4. In this model the dielectric material contains an assembly 
of charged particles whose interaction can be neglected. Each 
charge has two equilibrium positions separated by a high 
potential barrier. It is then assumed that each particle collides 
with the surrounding medium and that the number of collisions 
per second is such that the average time between two collisions 
is small compared with the average time r which a particle 
spends near one of its equilibrium positions before jumping to 
the other one. This leads to the linear diffei ential equation 9.1 1 
for the difference of the number of particles occupying the two 
positions. Its solution 9.15 approaches its equilibrium value 
exponentially with a relaxation time (cf. 9.5) 

T = ^ = -e"/«’, H>kT, T„<T, (11.1) 

21^21 ^0 

where H is the height of the potential barrier and the 

frequency of oscillation about either of the equilibrium positions. 
For this model, according to 9.5 and 9.6, the transition pro- 
babilities w ^2 ^^21 ^ particle between the two equilibrium 

positions are equal in the absence of a field, but are altered 
slightly when an external field is applied. 

It should be pointed out that from the assumptions made it 
follows that the derived exponential law can be true only on an 
average over time intervals which are large compared with Tq, 
Therefore the derivation of the Debye equations from tliis law 
can hold only if the period of the field is large compared with Tq. 
It must, therefore, be concluded that this model leads to the 
Debye equations only when the following conditions hold: 

To<T, To<l/a>. (11.2) 

Thus, in order to ascertain whether or not the Debye equations 
hold in the main region of absorption (near cdt ^ 1), it would 
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be necessary to determine and this in turn would require a 
detailed knowledge of the interaction between the charges and 
their surroundings. It will also be seen from conditions 11.2 
that even if the Debye equations do hold in this main region 
deviations may be expected at higher frequencies. 

It is important to realize that an exponential decay function 
(and hence the Debye equations) can only be obtained when the 
particles can be considered as independent of each other. Inter- 
action between particles (which is to be distinguished from that 
between a particle and its surroundings) would mean tliat the 
transition probabilities are not constant, but depend on the 
positions of the neighbours. Hence equations 9.7 and 9.8 would 
be no longer linear and, tlierefore, could not be solved by ex- 
ponential functions. 

Dipolar solids 

Models more directly applicable to solids, but leading to a 
mathematical treatment similar to the above, can easily be 
devised. Consider, for instance, the model for dipolar solids 
used in § 8. It consists of dipolar molecules, each of which, owing 
to the crystalline field, has a number of equilibrium positions 
with different dipole directions, which are separated by potential 
barriers. In the simplest case only two equilibrium positions 
with opposite dipole directions exist. In such a model at low 
temperatures the dipoles, because of interaction with each 
other, form an ordered arrangement. At a temperature Tq an 
order-disorder transition occurs and for T > Tq long-distance 
order vanishes. Short -distance order, i.e. order relative to 
neighbours persists, however, and can be neglected only at still 
higher temperatures. It will be shown that at these temperatures 
the Debye equations hold. They should be invalid at tempera- 
tures near Tq, but become valid again if T < Tq. 

Let us consider first the high -temperature region. Then in 
the absence of a field the lowest energy -level of a dipole is the 
same for both equilibrium directions. To carry out a transition 
between them a minimum energy, say J7, is required to lift 
the molecule over the potential barrier separating the two 
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equilibrium positions. In the calculation of the probability of 
a transition between tlie two positions equation 9.5 cannot 
now be applied because in general a molecule is also capable of 
internal excitation, a fact which was not considered in § 9. We 
shall take account of this by introducing a factor A which only 
varies slowly with temperature. Then if 1/2 t is the transition 
probability, 11.1 will be replaced by 

T = — (11.3) 

where is the average time required by an excited molecule 
to turn from one ecpiilibrium direction to the other. In the 
presence of a field f the molecule has an additional energy -[-(xf 
and “ fjif for the two equilibrium directions of its dipole respec- 
tively. Tlie minimum energy required by tlie molecule to enable 
it to carry out a transition in the presence of a field is, therefore, 
Thus, assuming 

], (11.4) 

the two transition probabilities are given by 

respectively. Hence if and are the numbers of dipoles in 
the two directions, equations 9.7, 9.8, and 9.1 1 still hold, provided 
these new values of ^’21 used. Therefore inserting 

11.5 and 11.6 into 9.11 yields 

(U.7) 

For constant fields this equation leads to exponential approach 
to equilibrium and hence according to § 10 to the Debye equa- 
tions with a relaxation time r. 

Assuming that a molecule interchanges energy very quickly 
with its surroundings, it can be considered as passing through 
a number of levels above the energy H during the transition 
between the two equilibrium positions. A quasi-thermal 
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equilibrium of the excited molecule with its surroundings is 
therefore established. Under these conditions the probability 
of finding a molecule with an energy above H is approximately 
(Dj^lDQ)ex:p{~HlkT), where is the number of energy-levels 
in a range of the order kT near the ground level and is the 
corresponding number in the excited state above H. Now the 
transition probability 1/2t must be equal to this quantity 
multiplied by 7r/a;„, which is the average time required by an 
excited molecule to turn by an angle 77, and divided by 2 because 
only half of the molecules move in the required direction. Hence 
by comparison with ecpiation 1 1.3 it is found that ap])roximately 

( 11 . 8 ) 

Expr^ission 11.3 for the relaxation time is sometimes (e.g. 
Frank, F2\ Kauzmarin, Kl) related to similar expressions for the 
rate of unimolecular chemical reactions. In general a formula of 
this type will be obtained for any process which requires excita- 
tion to an energy //, but calculation of the absolute value of a 
rate of reaction (requiring the factor A jeo^) is difficult. The first 
calculation of this kind is due to Pelzer and Wigner [P4 1. Usually 
can be assumed to be of the order 


o)^ ~ 10^2-10^^ per second, (H-^) 

so that measurement of the temperature -dependence of r leads to 
semi-empirical values for the order of magnitude of A (11.8). 
Sometimes (e.g. Eyring, El, E2) it is assumed that the trans- 
lational or rotational motion of the excited molecule can be 
separated from the other types of motion. In this case if 
kT, there are about kTjUa)^ energy -levels, in an interval 
kT, connected with the rotation of the molecule as a whole. 
Denoting by the number of energy -levels due to internal 
excitation of the molecule, we then obtain, using 11.8, 


D, 


'H 


)+ 


i.e. 


A _ h 


( 11 . 10 ) 


In view of the assumptions that have been made in the deriva- 
tion of this expression one should not expect it to give more than 

4980.11 a 
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an order of magnitude. And if '> kT, result might be 
obtained which is not even of the correct order (cf. Pelzer, 
PS), 

Let us now consider our model at lower temperatures. In the 
neighbourhood of the transition temperature interaction 
between dipoles becomes important. The energy of a molecule 
then depends on the dipole directions of its neighbours which 
invahdates the assumptions which led to the Debye equations. 
It will be shown presently, however, that this is no longer the 
case when the temperature is well below the transition point 
and the dipoles are therefore in an ordered state. As pointed 
out in § 8, the lattice then contains two types of sites; and if the 
state is completely ordered, these sites are occupied by molecules 
with opposite dipole directions. They may form, for instance, 
the centres and the corners of a body-centred cubic lattice. Each 
molecule has a second equilibrium position with opposite dipole 
direction in which its average energy is higher by an amount 
V{T) than in the original position. If the temperature is suffi- 
ciently low (well below T^), it is possible to replace V{T) by F(0), 
the value of F when T — 0. It is in this region that we can 
assume the energy of a dipole to depend on its own direction 
only, and not on the directions of its neighbours, because the 
latter are nearly always in the state of lowest energy (ground 
state). Hence we expect the Debye equations to hold. 

To prove this it is sufficient to show that equilibrium is 
approached exponentially with time since, as shown in § 10, 
the Debye equations then obtain. It will simplify the calcula- 
tions if we assume at a given time the existence of a dipole 
moment of the substance in the absence of a field and calculate 
its variation with time. Let us distinguish quantities referring 
to the two types of site by the use of + and — indices. Thus, 
for instance, probability per second for a transition of 

a dipole on a + site from the ‘ 1 ’ into the ‘ 2 ’ direction. Similarly 
BXidNp are the numbers of dipoles in the ‘L -direction for the 
two types of site. The total dipole moment is proportional to the 
quantity 


^N = N^-Nt+Np^N^. 


( 11 . 11 ) 
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Also if N is the total number of dipoles 

N^:+N^ - Nr+N^ - i-A. (11.12) 

Furthermore, in the absence of a field the 'T ('2’)-direction 
plays the same role for +sites as the '2’ (‘ l’)-direction for 
—sites. Thus in equilibrium at T “ o, — iV^ = and 
Af - ^ 2 ^ - 0. Also 

'^21 ^ (11.13) 

Thus by the same method as that used in deriving equations 9.7 
and 9.8 we now find 


dN^ 

dt ““ dt 




(11.14) 


from which, with the help of 11.12, we obtain 

±(N+-N+) - -2w+N+ + 2w+Nt 

= — '^^^ 12 )^' (11.15) 

The same equation holds for the — sites, if we exchange 1 for 
— indices. Using equation 11.13 this means that 

— (iVi A2 ) " + ^2 )““2(^^’21 ^2)^'* 

(11.16) 

Therefore by adding 11.15 and 11.16 and making use of 11.11 
we find 


~ AA, i.e. AA oc 

at 

(11.17) 

whichprovesthat AA approaches its equilibrium value (AA — 0) 
exponentially. 

Dipolar liquids 

It was shown in § 8 that when a disordered solid melts there 
should be no appreciable change in the static dielectric constant 
because in both the liquid and solid phases the average energy 
of a dipole is the same for all equilibrium directions. The fact 
that — ^in contrast to liquids — there is only a discrete number of 
equilibrium directions in solids has no influence on the static 
dielectric constant. 
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There is, however, an essential difference in the dynamic 
properties of solids and liquids. In solids in view of the interaction 
of a molecule with its Neighbours a dipole has a number of equi- 
librium directions. They are separated by potential barriers 
over which the dipole must pass in turning from one such direc- 
tion to another. In liquids the average distance between neigh- 
bours, and hence the interaction, is about the same as in solids. 
Therefore, if the positions of all molecules but one could suddenly 
be fixed, this selected molecule would behave similarly to a 
dipolar molecule in a solid, i.e. it would probably possess a 
number of equilibrium directions separated by potential barriers. 
It is an essential property of a liquid, however, that its mole- 
cules have no fixed positions. Thus, if we imagine one molecule 
to be turned out of its momentary equilibrium directions, its 
neighbours will tend to rearrange themselves in such a way as to 
make this new direction an equilibrium position. So we may 
have two conceptions of the way in whicli a dipole in a liquid 
alters its direction. It may jump into a different direction in a 
similar way as in solids. This requires — at least for the duration 
of this jump — that the arrangement of its neighbours should 
remain unaltered. The second possibility is that such jumps may 
occur very rarely and that a dipole may alter its direction only 
in conjunction with a rearrangement of the positions of its 
neighbours. In fact a dipole might be considered as fixed fairly 
rigidly relative to its neighbours. The rotation of a dipole would 
then affect the motion of molecules at some distance from it. 

The average motion of these neighbouring molecules might be 
described by replacing them by a continuous medium with the 
properties of a macroscopic viscous fluid. This possibility leads 
to the model used by Debye \D2] in which a dipolar molecule is 
considered to be a sphere of radius a moving in a continuous 
viscous fluid with viscosity rj and obeying the macroscopic 
equations of flow. The fluid is considered to adhere to the sur- 
face of the molecule. On these assumptions the frictional 
constant ^ of the sphere is given by Stokes’s law, 

I = STTTja^. 


( 11 . 18 ) 
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This means that when an electric field f, making an angle d with 
the direction of the dipole ft, is applied, 

= -~iu/‘sin(9, (11.19) 

if no other forces act on the molecule. In this equation it is also 
assumed that inertia effects of the dipole — which macroscopi- 
cally would lead to an additional term IdWIdf^ [1 = moment of 
inertia) on the left-hand side — can be omitted. This omission is 
justified in so far as this form of the term is concerned since it 
would hardly be expected that such effects could be described 
adequately by this macroscopic expression. 

The macroscopic equation of flow 11.19 can yield correct 
values for 6 only if account is taken of thermal fluctuations due 
to the interchange of energy between the molecule and its 
surroundings. For the above model the energy of a molecule is 
the kinetic energy of rotation of the sphere and has an average 
value of order kT. This rotation — corresponding to Brownian 
motion in translation — changes its direction and magnitude very 
frequently due to collisions with the molecules of the surround- 
ing liquid. However, if no external forces act on the dipole its 
average displacement from a given direction vanishes because 
displacements in all directions have equal probability. The 
mean square of these displacements, on the other hand, will 
increase steadily with time. 

Equation 11.19 may, therefore, be considered to hold for the 
average value of d (for a single dipole) when a field is applied. 
The actual value of 6 may, however, only be expected to be near 
its average value as long as the mean square displacement due 
to thermal fluctuations remains a small quantity. For a macro- 
scopic dipole this should practically always be so. According 
to 11.19 such a dipole will, therefore, gradually approach the 
direction 0 0 parallel to the field. In equilibrium it will remain 

in this direction apart from small fluctuations. 

It is to be noted that the value 6 used in 11.19, and assumed 
above to be the time average of the 0-values of a single 
dipole, may also be considered as the average of such angular 
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displacements of a large number of (non-interacting) dipoles all 
having the same 6^- value within a small range. 

In contrast to a macroscopic dipole the fluctuations of a 
molecular dipole are very large. This follows at once from con- 
sideration of the equilibrium distribution of an assembly of 
dipoles which, according to 6.9 (replacing E by/), is proportional 
to exp(ju/cos0//:T). Since ixf<^kT for all practicable fields 
(cf. § 4), the field has only a slight influence on the distribution 
function. In other words, a dipole has only a slightly higher 
probability of being near 6^0 than in the opposite direction. 

In order to derive now the macroscopic dielectric properties 
of dilute solutions of dipolar molecules in non-polar liquids, 
consider an assembly of non-interacting dipoles and let 

N{d,t)m\6d6 

be the number of dipoles per unit volume in a range dd near 6 
at time t. The function N{6,t) in general varies with time because 
individual dipoles continuously change their directions. These 
changes are either due to thermal fluctuations or to the action 
of the field. The latter, according to 11.19, causes all dipoles 
to be displaced by an amount 80, 

he=^~U 'fnineU, (]J.2(i) 

(it g 

ill a short time interval 8t. Tlierefore within 8^ seconds a number 
iV(0, ^)sin 9 80 will pass through the surface of a cone with angle 0. 
Owing to the action of the field the rate of change of the number 
of dipoles within an interval dd near 0 is, therefore, using 1 1 .20, 
given by 

( 11 . 21 ) 

This equation follows because a number A^(0,^)sin0 80 wiU 
leave the range dd within St seconds, whereas the number of 
dipoles entering it is given by the same expression replacing 0 
by d—dd. 
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From 11.21 it follows that the rate of change of the function 
N(6, t) itself due to the action of the field is given by 

In the case of large fluctuations with which we are dealing at 
present this equation is more suitable than 1 1.11> to express the 
effect of the field on dipoles. 

Equation 11.22 can be used to find the rate of change of the 
dipole moment il/y in the field direction, of the assembly of dipoles 
under the influence of the field. Assuming the dipoles to be rigid, 

TT 

Mf fjL J cohO N(6,t)^md dd (11.23) 

0 

per unit volume. Therefore using 11.22 the rate of change of 
due to the action of the field is given by 

0 

TT 

= ~ J Gosd-^{N(e,t)fimW}de, 

0 

or integrating in parts, 

rr 

= Y j if)sin*0}8in 6 dd = mipg, 

** (11.24) 

where Nq is the total number of dipoles per unit volume and 
sin^ is the average value of sinW. As discussed above, the field 
influences the angular distribution of dipoles only very slightly. 
Therefore the value of sin?5 in equilibrium in the absence of a 
field can be inserted in 11.24, i.e. 


Hence 


TT .IT 

sia^d = j ^inW sind d6 / j &md dd ~ 

0 f Q 



2 

3 ^ ■ 


/ 

(11.25) 

(11.26) 
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The total rate of change of Mf contains in addition to 1 1.26 a 
term due to thermal motion. This term will tend to restore the 
equilibrium distribution in the absence of a field, in which case 
Mf — 0. On the assumptions which we have made, namely (i) 
that there is no interaction between dipoles, (ii) that in a short 
time interval 8^ the 0- value of a dipole is altered only very slightly, 
this second term must be a linear function of the deviation from 
equilibrium for / = 0, i.e. it must be proportional to —Mf. 
Therefore introducing 1/t as proportionality factor, the total 
rate of change of Mf is given, using 1 1.26, by 


^ 2 ^ 
r ^ ‘ 


(11.27) 


In view of the linear relationship between Mf and dMfjdt 
equation 11.27 leads, of course, to exponential approach to 
equilibrium and hence, according to § 10, to the Debye equations. 
This result in fact should always be expected when the following 
conditions hold: 

(а) Absence of interaction between dipoles. 

(б) Only one process leading to equilibrium (e.g. either transi- 
tion over a potential barrier, or frictional rotation). 

(c) All dipoles can be considered as in equivalent positions, 
i.e. on an average they all behave in a similar way. 

The value of r can be derived from 1 1 .27 by making use of the 
equilibrium value of Mf in the presence of a field. Thus, using 

M, — (11.28) 

This expression must be equal to the equilibrium value for 
Mf resulting from 11.27 with dMfjdt — 0, i.e. 

Mf^^Tl^. (11.29) 

Hence comparing 11.28 and 11.29 and using 11.18, 

This important formula due to Debye [D2] is often used 
to discuss the relaxation time of dipolar liquids. It seems 


(11.29) 


(11.30) 
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appropriate, therefore, to recapitulate the assumptions under 
which it should hold. First of all the above-mentioned assump- 
tions (a), (b), and (c) require 

(a) dilute solutions of dipolar molecules in a non-polar liquid; 

(j8) axially symmetric molecules; 

(y) isotropy of the liquid — even on an atomic scale in the 
time average over an interval small compared with r. 

Here (a) follows from (a), and (j3) and (y) are both connected 
with (b) and (c). All three assumptions are necessary for the 
Debye equations to hold ; but they are not sufficient to obtain 
the value 1 1 .30 for the relaxation time. This relation, as has been 
discussed above, is based on the further assumption that the 
dipolar molecule is fixed fairly rigidly relative to its neighbours 
so that large jumps of the di|)ole direction are unlikely. Now 
the empirical temperature -dependence of the viscosity, 

(11.31) 

suggests that jumps over a potential barrier of height are 
carried out by the molecules of the liquid in processes connected 
with viscous flow. Therefore if H is the height of the potential 
barrier related to jumps of the dipolar molecule, the condition 
that these jumps happen only very rarely suggests that 

( 8 ) 

This condition, however, is correct only if the cf)efficients 
A (cf. 11.3) have the same order of magnitude for both types of 
transition. In general it should be expected [F12, 5^3] that low 
viscosity liquids have small values of and they are more likely, 
therefore, to satisfy (8) than are liigh viscosity liquids. 

Formula 11.30 can thus be used to find the temperature 
dependence of r which is essentially the same as of 07, i.e. 

TOce^V^^’, (11.32) 

because the other terms vary only slowly with T. It follows that 
this temperature-dependence is independent of the nature of 
the dissolved molecules and is a function of the viscosity of the 
solvent only. To obtain the absolute value of r as well would 
require knowledge of the effective radius a which is the radius of 
a solid sphere having the same frictional constant as the dipolar 



90 DYNAMIC PKOPEBTIES HI, § 1 1 

molecule. No theoretical investigations on the correlation of 
this radius to the molecular radius normally used have yet been 
carried out. One should expect, however, that a depends on 
both solvent and solute and, therefore, cannot be assumed to be 
a molecular constant. Thus at the present stage of development 
equation 11.30 camiot be used to obtain absolute values for the 
relaxation time. 

12. Generalizations 

The derivation of Debye’s relation 1 1 .30 — connecting the re- 
laxation time of a dipole molecule with the viscosity of the 
liquid in which it is dissolved — required the assumption tliat the 
dipolar molecule is bound so strongly to the surrounding mole- 
cules that large jumps of the dipole direction are very unlikely. 
This may be true for a number of cases, but others may exist in 
which the opposite is more likely. A dipolar molecule will then 
make many jumps over the potential barrier separating it from 
another dipole direction during the time required for an appreci- 
able change in direction by viscous flow. Clearly this holds for 
solids where flow may be considered as entirely absent; but some- 
times it may also be expected in liquids, and in particular in 
amorphous substances for which the viscosity is so high that 
flow is practically negligible. In liquids it might also happen that 
the process which prevails is different for different kinds of 
dissolved molecules. A further possibility mentioned by Schalla- 
mach [S3] is the coexistence of both types of transitions. 

As an example of the type of substance exhibiting the second 
type of behaviour mentioned above we may consider a dilute 
solution of dipolar molecules in a liquid or in an amorphous solid. 
In this case, in contrast to § 11, we assume that the dominant 
process for changing the dipole direction is that involving many 
large jumps. A single dipole will then behave similarly to a 
dipole in a crystalline solid; as in 1 1.3, the probability for a jump 
over the potential barrier (height H) will be proportional to 
exp(--HlkT), In contrast to conditions in crystalline solids, 
however, the arrangement of the nearest neighbours is not 
exactly the same for all dipoles, and hence the heights H of the 
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potential barriers will also differ, thus producing different 
values of the transition probabilities. The dipolar mole- 
cules in such a substance may then be classified according 
to the height H of their respective potential barriers. If the 
substance has been polarized, and the external field has been 
removed, the contribution of the dipoles in a small range of 
energies near H will decay exponentially — as shown in § 11 — 
with a relaxation time r related to H by equation 11.3. There- 
fore, instead of H the individual relaxation time r may be used 
to classify the molecule. Let ^(t) dr be the contribution to the 
static dielectric constant of the group of dipoles having individual 
relaxation times in a range dr near r. Since interaction between 
dipoles can be considered to be absent (dilute solution), the 
contributions of the various groups superpose linearly. Their 
total contribution to the static dielectric constant is, therefore, 
given by ^ 

^ J y{'r)dT. ( 12 . 1 ) 

0 

The function y{r) wliich describes the distribution of relaxation 
times will be called the distribiition function. 

To obtain the complex dielectric constant we first consider 
the decay function oi{t) (cf. §§ 2 and 10). The dipoles with relaxa- 
tion times in a range dr near r make a contribution to a(^) which 
is proportional to exp( —tjr) and to y{r)drjr, which corresponds 
to the coefficient of the exponential term in equation 10.8. 
Therefore the total contribution of all the dipoles is given by 

00 

«(/)- ( 12 . 2 ) 

0 

Using relation 0.1, the complex dielectric constant is now 
obtained from 12.2, 

00 GO CO 

J* 0 L{x)e}^*^ dx == j dx J* — e~^l'^y{r) 

0 0 0 

dx (12.3) 

0 0 
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The last step is due to an interchange of x and r integration. 
The -integral is then identical with the integral in 1 0. 1 4. There- 


e(co)—€^ 


f y{^) dr 

J 1 — icor’ 


(12.4) 


or separating real and imaginary parts according to 2.8 similarly 
to 10.15 and 10.16, 


i \ „ r 

eiH 

0 


(12.5) 


( 12 . 6 ) 


For a substance of the type described in this section these equa- 
tions will replace the Debye equations (10.15, 10.16, 10.17), and 
the latter will not be satisfied. A clearer insight into the meaning 
of equations 12.5,12.6 will be obtained by a discussion of the way 
in which their solutions deviate from those of the Debye formulae. 

These deviations are best considered in connexion with the 


shape of the power loss — frequency curves e^io}). For a detailed 
discussion, knowledge of the distribution function y{r) is, of 
course, essential. As a preliminary step, however, we may note 
that y{T) is always positive and that therefore consists of a 
superposition of Debye curves a>T/(l (cf. Fig. 13) with 
different positions of their respective maxima. The resulting 
curve € 2 ( 0 )), therefore — supposing it has a single maximum — 
has a larger half-width than that of a single Debye curve whose 
maximum coincides with that of € 2 ( 00 ), As a simple example 
[F7] consider a model in which each molecule has two equilibrium 
positions with opposite dipole directions, and with equal energy 
in the ground level, as in the case of disordered solids; but in 
contrast to the latter the potential banier between the two 
positions has a different height for each molecule. It will be 
assumed that the heights H of the potential barriers are equally 
distributed over a range between Hq and i.e. 
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Thus if iVJ) is the total number of dipoles per unit volume, 

iV„- (12.8) 

is the fraction with fl^-values in a range dv near //o+^* 

For dilute solutions interaction between dipoles can be 
neglected. As shown in §§ 4 and 6, the contribution of a dipolar 
molecule to is then independent of i/, and therefore is the same 
for all molecules, i.e. per molecule. The individual 

relaxation time r, however, depends on H according to 11.3 and 
11.8. Hence, using 12.7 and considering A as constant, 

T -= T„ = (12,9) 

The individual relaxation times r therefore cover the range 

To T < Tj where Tj (12.10) 

To determine the distribution function yir) we note that 
7 /(t) — 0 outside the range 12.10. Now consider r as a function 
of V so that from 12.1, using 12.9, 

Ti tJo Vo 

= J ?/{t) = J ^ J 

To 0 0 

( 12 . 11 ) 

This means that y{T)r dvjkT is the contribution of the molecules 
in the range dv to the static dielectric constant. On the other 
hand, we have just seen that this contribution per molecule is 
Coo)/A^o independent of v. Therefore, since 12.8 represents 
the number of molecules in dv, 

y{r)T dv _ „ dv 

~ kT -' No %o’ 
or 

y{r) = (es-e^)~-, if T„ r < = Toe'’./*:y 

y(r) = 0, if T < To and r > 

The above distribution function — corresponding to an equal 
distribution of potential barriers over a range Vq — depends on 


|. ( 12 . 13 ) 
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temperature. In particular its relative width (tj— To)/tq de- 
creases with temperature since 

TlUl? = f>rolkT___ I (12.14) 

To 

The dielectric constants and eg now be obtained by 



Fig. 14. Dependence of dielectric loss e 2 (o>) on angular frequency a> accord- 
ing to equation 12.19 for the three values 1, 5, and 10 of the parameter 
They correspond to a range of heights of potential barriers of width 
Vg = k7Uog € 2 {o))l€ 2 {o)^) is plotted against a>/a>^ on a logarithmic scale; 
represents the value of a> where eg reaches its maximum. 


insertihg ^(t) from 12.13 into 12.5 and 12.6. All integrations can 
be carried out in an elementary manner leading to 


€i(a))- 




.) 1- 


l+coVge2««/*^\ 


■ 


(12.15) 


— {€g — €oo) — (tan ^(cutq — tan“^coTo). (12.16) 
^0 

These equations replace the Debye formulae 10.15 and 10.16 
for our present model. In this case ei(cu) and € 2 ( 0 )), considered 
as functions of frequency, depend on two parameters — a relaxa- 
tion time Tq, and the factor vJkT determining the width of the 
range of relaxation times by 12.14. The Debye formulae are 
obtained for vJkT — 0; they contain, therefore, one para- 
meter only. The shape of eg can be seen from Fig. 14, where 




in, §12 GENERALIZATIONS 95 

€^{o))/€ 2 (co^) is plotted against co^ being the frequency at 

which €2 has its maximum. 

It can be seen that with increasing values of the second para- 
meter vJkT the curves become more and more flattened out. 
The position of tlie maximum is obtained from 12.16, using 

^^2 -r 

7— ^ = 0 II o) “ oj^. 

dct) 

This leads to == i-c-.-i’WfcT „ _J: — ^ (12.17) 

1 

€<;-e"oo 
0-5 


0 5 10 

Vo/kT 

Fiq. 16. Dependence of the maximum value of the loss, €a(ct)^)» on the 
width t’o of the range of heights of j^otential barriers. 

which for v^) — 0, of course, becomes identical with the corre- 
sponding value 10.22 for the Debye formulae. Inserting 12.17 
into 12.16 yields the maximum value of eg? 

hT 

€2(00^) (c^,— £00) — [tan“M^«/^^’— tan'^e-^^'o/*'^]. 

^'0 

(12.18) 

Again for ^0 — ^ quantity becomes equal to 10.23, but it 
decreases with increasing vJkT, as shown in Fig. 15. 

With the help of 12.14, 12.17, and 12.18 equation 12.16 can 
be written in the form 

which shows clearly the dependence of € 2 ( 0 )) on the parameter 
Tq/ti < 1. For tJti = 1, equation 12.19 is identical with the 
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corresponding equation of § 10, and the smaller is, that is 
the larger the larger is the deviation from a Debye curve. 

The above equations aie based on the assumption of equal 
distribution of the height of potential barriers over a range 
leading to a range of relaxation times given by equation 12.14 
and to the distribution function 12.13. This may seem to be 
rather a specialized assumption. It can be shoAvn quite readily, 
however, that the behaviour of in the main absorption range 

is altered but httle, if the distribution function 12.13 is replaced 
by any other smooth function which is large between and 
and small outside this range. Quite generally, therefore, in the 
neighbourhood of the maximum of € 2 ( 0 ;) this function can be 
considered to be determined by two parametei’s, namely the 
angular frequency at which the maximum occurs, and the 
width of the range of relaxation times It should be 

remembered that the absolute value of is then determined 
by €00 with the help of the phenomenological relation 2.18. 

In contrast to the behaviour of € 2 ( 0 ;) in the main absorption 
region the behaviour outside this region may be very sensitive 
to small changes of the distribution function y{r) for r-values 
outside the range tq < t < which contains the bulk of the 
relaxation times. For the latter, though very numerous, make 
only small contributions to if cotq 1 or cdt^ ^ 1. These 
may well be smaller than the contributions of a few relaxation 
times near r ~ l/co which make their maximum contribution 
at the frequencies outside of the main absorption band. 

For many substances the value of € 2 (a») does not approach 
zero outside the main absorption range, but remains at a very 
small value which changes very slowly with frequency [02], 
This may be explained if a sufficiently wide range of relaxation 
times is used. Garton [01] has suggested that this residual 
absorption may be due to the temporary formation of additional 
equilibrium positions of the molecules due to thermal fluctua- 
tions. The probability for the occurrence of such a potential 
well of depth in a range dw (cf. Fig. 16) is assumed to be pro- 
portional to 




( 12 . 20 ) 
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This is a plausible assumption, though no proof for it has been 
given. We shall proceed to sliow that such an assumption does, 
in fact, lead to a value of € 2 ( 0 )) which is substantially independent 
of frequency for frequencies well outside of the main absorption 
range. 

Let H be the height of the potential barrier measured from a 



Fig. 10. Potential well of depth w duo to thermal fluctuations. 

normal (permanent) position and consider a molecule near which 
a temporary position of depth w has been formed. According to 
11.3 the time it spends in the temporary or in the permanent 
position is given by 

T == and Tq = (12.21) 

respectively on the further assumption that the temporary 
position exists for a time which is long compared with t+To- 
The coefficients B and Bq are considered to be independent of w. 
Assuming r ^ the relative probability of finding this mole- 
cule in the temporary position is given by T/(T-f Tq) t/tq. The 
average number of molecules in such positions, using 12.20 and 
12.21, is thus proportional to 

L^-wlkTdw = - - ~dr = -kT— oc — . (12.22) 

Tq Tq t dr Tq t r 

Each molecule makes a contribution to proportional to 
a>T/(l+^^T^)* Also since > 0 the smallest value that r can 
take is T = J5. For large values near r = Tq the developments 
become invalid, and tq will be used as the upper limit for r; this 

4980.11 H 
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should not be an essential limitation, however, because near 
r = tq the main contribution to is due to the transitions 
between permanent positions which are much more numerous. 
Thus 

T« 

€2(co) oc J ~ = tan-iwTo-tan-i^w (12.23) 

B 

for frequencies higher than the main absorption range, ^ 1. 
If, in particular, 

1/t, < CO < 1/5, (12.24) 

then is practically independent of co. Here IjB will 

probably be a very high frequency. 

The preceding discussions w^ere all based on the assumption 
that interaction between dipoles can be neglected. It seems 
likely that the influence of interaction wall always tend to 
broaden the Debye absorption curve, because the energy of a 
molecule in an equilibrium position may have a whole range of 
values depending on the positions of its neighbours. One would 
thus expect the heights of the potential barriers to vary accord- 
ingly. In the present case, however, the time which a molecule 
spends in a given equilibrium position is of the same order as 
the time during which the height of the potential barrier remains 
constant. If an attempt is made to include interaction under 
these conditions, great mathematical difficulties are encountered 
and no solution for this case has yet been obtained. 

13. Resonance absorption 

In § 9 it was pointed out that two types of power loss should 
be expected to occur in dielectrics: (i) The loss due to displace- 
ment of charges bound elastically to an equilibrium position. 
Such charges have a proper frequency of oscillation, say ojqI27t, 
and the power loss (and hence Cg) will be expected to have 
a maximum near this frequency (the resonance frequency), 
(ii) The loss due to transitions of charges or dipoles between 
equilibrium positions separated by a potential barrier. Such 
transitions are described by a relaxation time r, and in this case 
the power loss will have a maximum near a frequency 1/27 tt as 
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discussed in §§ 10-12. This frequency is usually strongly depen- 
dent on temperature, in contrast to the resonance frequency in 
case (i). The chief aim of the present section is to derive expres- 
sions for the complex dielectric constant when there is resonance 
absorption (case (i) above). 

It is our intention to consider only the simplest case, and it 
will be seen presently that the required formula can be derived 
from very general considerations without choosing a particular 
model. We shall employ the same method that has been used in 
§ 10 to derive the Debye equations. The polarization of the sub- 
stances considered there was assumed to approach equilibrium 
exponentially with time. This cannot be exj^ected to hold in the 
present case of elastic displacement. One would rather expect 
damped oscillations with frequency about the equilibrium 

polarization, as has already been discussed in § 9. Therefore 
instead of 10.1 we shall now assume a decay function 

(x{t) == ye~^l'^coB(ojQt-\-ip)y (13.1) 

and then, making use of equation 9.1, derive an expression for 
the complex dielectric constant €. Equation 13.1 contains two 
constants y and ip which will be determined below. 

The reader is reminded that a(t) is a macroscopic quantity 
referring to the behaviour of the polarization of the whole sub- 
stance. It would be wrong to describe the amplitude of the 
oscillations of a single molecule (in the absence of a field) by a 
function similar to 13.1, because in equilibrium this would lead 
to vanishing amplitude and thus would not take account of the 
thermal motion. 

By inserting 13.1 into 9.1 w^e now find for the complex 
dielectric constant €, 

00 

^—^00 = y J e~^^'^cos(coQX+i/j)e^'^^ dx 

0 

yr / 

2 \1 — i(6t>o+^)'^"^ l+^(^o — ^)' 

l-f^tan0 ^ 1— itan^ \ 2) 

1 — i(a>Q + Cu)T ] +|(cOq — 6u)t/ 
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As the next step, the two constants y and ^ will be determined 
by considering the limiting cases of very low and very high 
frequencies. Low frequencies, w coq, are considered to be 
outside the main range of absorption, i.e. in a region in which 
€2 is very small and, therefore, the dielectric constant is prac- 
tically real. In this range e will approach a constant value, say 


€ = €ao~i-A€, CO OJq, (13.3) 

where Ae is a real (piantity. This value of the dielectric constant 
will hold down to frequencies at which another absorption range 
starts. Fn the absence of further absorption regions 13.3 repre- 
sents the static dielectric constant. The high-frequency dielectric 
constant on the other hand, is the value of e at the high- 
frequency side of the al)sor|)tion peak, again in a region in which 
€2 is very small. 

Comparing 13.3 with ecjuation 1 3.2 and neglecting co compared 
with coq, we find 


Ae 


cos l/j—C 


yr 


TsinJj , 1— coArtanj/f 

^ == yr cos ifj ^ ^ 


l+colr^ 




(13.4) 


as a first condition for the determination of y and i/j. 

The second condition is of a more subtle nature and concerns 

the behaviour of e at angular frequencies oj ^ co^y At these 

frequencies cuq can be neglected compared with co, so that 13.2 

becomes , 

cos lU . / 1 o ^ \ 

€—€^ = yT- — r— , to > COq. (13.5) 

1 — 2 COT 


On the other hand, for co ^ coq the influence of the restoring force 
should be negligible during one period 1/co. Therefore a be- 
haviour according to the Debye theory should be expected with 
a static dielectric constant given by 13.3. Thus by replacing 
^00 by Ae in 10.13, 

g— goo = ~~ , CO>COo. (13.6) 

1 — tOJT 

Combining 13.5 and 13.6 we find 

Ae^yrCOSip (13.7) 

as the second condition for the determination of y and i/j. 
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From 13.7 and 13.4 we now conclude that 

ojqT = — tant/r. (13.8) 

The insertion of this expression and 13.7 ipto 13.2 finally leads 
tof (cf. references V 5 and F9) 


\1 — ^(a>o+ai)T l+i(ojQ~-aj)rJ 
or separating real and imaginary parts according to 2.8, 


I+topT \ 
aj)T/’ 


(13.9) 






l + (a>-a>o)2 


(13.10) 


l + («>+‘"o)^'^^ l+(w— a»o)V 


(13.11) 


These equations represent the dielectric constants in the 
simplest case of resonance absorption. In many practical cases 
Ac is a small quantity so that is approximately equal to e^)- 
In these oases the loss angle according to 2.5 is given by 


tan^ = ^c=!^, ifA€<l. (13.12) 

^00 

If Ae is small it will be difficult to determine this quantity 
experimentally with reasonable accuracy by direct measure- 
ment. Following 2.18, however, the relation 

00 

Ae — - I = — I €o^/(log6L>) (13.13) 

TT J “ CD TT J “ 

0 

can be used instead, "i'his equation holds for tlie present case, 
as can be checked by inserting eg from 13.11. 

We shall now proceed to discuss power loss (oc eg) in more 
detail. At constant temperature eg as a function of cd has a maxi- 
mum when o) — cD^, 

+ (13.14) 

T 

as can be found from the condition dejdo) — 0. The maximum 
value of eg is then given by (inserting 13.14 into 13.11) 

^2(^m) = iAe^(l-f togr^) = JAew^^r. (13.15) 

t Another derivation of 13.0 is given in the Appendix A 4» 
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Consider now that the time of relaxation t depends on tem- 
perature and may be expected to decrease with increasing 
temperature. 

Then 

a>^ l/r and €^{w^) if ojqT < 1, (13.16) 

i.e. at very high temperatures. But 

and 62 (co,J iA€a>oT if > 1, (13.17) 

i.e. at lower temperatures. Comparison with 10.23 shows that 
case 13.16 leads to behaviour similar to that for Debye absorp- 
tion. Without detailed consideration of a special model it is 
impossible, however, to say whether this region can be reached 
at temperatures corresponding to the solid state. The low- 
temperature limiting case 13.17 leads to typical resonance 
absorption. In contrast to Debye absorption, the angular 
frequency cjq of maximum absorption is independent of tem- 
perature, but the resonance peak becomes narrower and higher 
the lower the temperature, i.e. the larger ojq t. Figures 17 and 18 
give a comparison of Debye absorption and resonance absorp- 
tion for various temperatures. 

A number of models have been considered in greater detail, 
and equations 13.9-13.11 have been found to hold. These deriva- 
tions are much more complicated than the one given above, 
but they may be considered to be more rigid (cf. Appendix A 4). 
Transitions between rotational levels of gases have been con- 
sidered by van Vleck and Weisskopf [F5] and by van Vleck [F4J. 
In gases the relaxation process is due to collisions between mole- 
cules. Therefore r will be expected to decrease with increasing 
pressure and with increasing temperature. At temperatures or 
pressures at which cug r ^ 1 the absorption peak according to 
13.16 is near the frequency = 1 /t. In this region, therefore, 
o) > (jJq, so that coq can be neglected. 

In the main absorption region, therefore, the Debye equations 
are fulfilled whenever coqt 1. At small pressures, however, 
coqT ^ 1 and the main absorption peak according to 13.17 will 
be found near the resonance frequency coq. 

Resonance absorption may also be of importance in solids 
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containing dipolar molecules with large moments of inertia 
\_F9, H3, Sll]. The resonance in this case is due to rotational 



Fig. 17. Frequency dependence of dielectric loss for three dift'erent tempera- 
tures Ti a T 2 <i T 3 (schematically), (a) Debye case ; {b) resonance case. 



Fig. 18 . Dependence of the angular frequency a>^ at which dielectric loss has 
its maximum on the relaxation time. D refers to the Debye case, B to the 

resonance case. 


oscillations of the molecules about their equilibrium positions. 
The behaviour of the relaxation time in this case, however, has 
not yet been investigated. 




CHAPTER IV 


APPLICATIONS 

The present chapter is designed to show by examples how the 
general theory developed in previous chapters can be applied 
in the discussion of the properties of dielectric materials. It is 
the intention to cover some important types of materials (e.g. 
non-polar solids, ionic crystals, dipolar liquids, etc.) ; this should, 
however, not be considered as an attempt to discuss systemati- 
cally the properties of dielectric materials, but rather as a 
selection of typical examples. 

14. Structure and dielectric properties 

Atoms 

In the present section we shall commence by discussing the 
structure and the dielectric properties of atoms and molecules 
leading up to a classification of dielectric materials according 
to their structure. An atom, as the reader will know, consists 
of a positively charged nucleus surrounded by electrons which 
just compensate its charge and whose combined mass is very 
small compared with that of the nucleus . By adding or removing 
one or several electrons one obtains negative or positive ions. 
In a stationary state the electric dipole moment of an atom (or 
of an ion relative to its nucleus) vanishes. 

A very general quantum mechanical theorem says that under 
the influence of external (static or alternating) electric fields the 
electrons of an atom behave like an assembly of classical har- 
monic oscillators whose frequencies and other properties can be 
specified, but will not be required by us. Usually the frequencies 
are equal to or higher than those of visible light, and, therefore, 
are much higher than the frequencies of electric fields contem- 
plated in the present book. Under the influence of an external 
electric field f an electric dipole will, therefore, be induced 
which has all the characteristics of a dipole produced by elastic 
displacement of electric charges (electrons) as discussed in case 
(i) of §§ 4 and 9. In the ranges of field-strengths and frequencies 



IV,§ J4 


STRUCTURE AND DIELECTRIC PROPERTIES 


106 


to be considered here, the strength of the induced dipole is 
proportional to the field-strength /, but independent of the 
frequency, and it shows no phase-shift relative to the field. These 
properties of an atom can be described by a constant polariza- 


bility 


= <xj. 


(14.1) 


The resulting polarization will be called optical polarization, 
to indicate that the characteristic frequencies are usually in the 
optical region. Often it is also described as electronic polariza- 
tion because it is due to the displacement of the electrons; this 
description may, however, be misleading as will be seen presently. 


Molecules 

Consider now a diatomic molecule A B composed of two atoms 
A and B. In view of the interaction between the atoms — leading 
to the chemical bond — the distribution of electrons in the mole- 
cule is different from a mere superposition of the electron distri- 
butions of the free atoms A and JS. The distribution should, 
however, be axially symmetrical around the line joining the 
nuclei of the two atoms. It must, therefore, be expected that a 
diatomic molecule will j)ossess a dipole moment in the direction 
A-B. An exception occurs in the case of two equal atoms, 
A = B, where the dipole moment must vanish for reasons of 
symmetry. Thus, for instance, the molecules HCl or CO have a 
dipole moment in contrast to Hg, Og, or Cl^- The magnitude of 
the dipole moment gives valuable information about the distri- 
bution of electrons. Thus the large dipole moment of HCl 
(/X 1 X l()~^%.s.u. — 1 Debye unit) comy^ared with that of 

(JO (/X ~ 0*1 Debye units) gives some justification for considering 
the HCl molecule (but not the CO' molecule) as composed of a 
positive and a negative ion, H++C1“. 

Under the influence of an external field f an average dipole 
moment m is induced in a molecule which, according to §§ 4 and 
9, can be considered as composed of contributions of two types: 
(i) those due to elastic displacement of charges, and (ii) those due 
to a change in the average orientation of the permanent dipole 
of the molecule. The former again are composed of the contri- 
butions of the various normal modes. In view of the great 
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difference in the mass of electrons and nuclei there is one group 
of normal modes which is practically entirely due to displace- 
ment of electrons relative to the nuclei. Their proper frequencies 
are usually in the optical region, visible or ultra-violet. The 
proper frequencies of the normal modes connected with 
oscillations of the nuclei, on the other hand, are usually in the 
infra-red region. They must contain displacements of electrons 
relative to their respective nuclei as well, because of the variation 
of the interaction between electrons and nuclei with the distance 
between the latter. Thus 

m = (14.2) 

where the first two terms are due to elastic displacements with 
proper frequencies in the optical and in the infra-red region, and 
is the moment due to dipole orientation. The quantities 
and are often denoted as electronic and atomic polarization, 
respectively, which may bo considered as somewhat misleading 
because the latter also contains terms due to electronic dis- 
placements. Since the frequencies contemplated in this book 
are much smaller than all proper frequencies, the relations 

m„ = ocj, m,., = (14.3) 

should hold with constant values of the polarizabilities and 
independent of the frequency of the field or of temperature. 
Polyatomic molecules behave in a similar way to diatomic 
molecules in that the average dipole moment induced by an 
external field can be considered as a superposition of three terms 
according to 14.2: two terms representing the elastic displace- 
ments with proper frequencies in the optical and infra-red 
regions respectively, and a dipolar contribution. The former 
can be obtained from the optical and the infra-red polarizabilities 
and according to 14.3; the value of the latter depends on 
the dipole moment of the molecule (e.g. § 6). The determination 
of the dipole moment is of great interest in investigations on the 
structure of molecules. For instance, the fact that COg does not 
possess a permanent dipole moment leads to the conclusion that 
the three atoms must be arranged in a straight line with the car- 
bon atom half-way between the two oxygen atoms. In contrast 
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to this the HgO molecule has a permanent dipole moment, which 
means that HgO is a triangular molecule. 

A great number of interesting applications to stereochemistry 
have been made, but a detailed discussion of this subject is 
beyond the scope of the present book. The reader will find 
examples and reference to the literature in the fairly recent book 
by Le Fevre [1^2] or in the article by Sutton [SIO]. The earher 
developments are to be found in the books by Debye [D2] and 
by Smyth [/Sf^]. 

Two quahtative rules concerning large molecules should be 
mentioned. Such molecules often contain a number of dipolar 
molecular groups such as the hydroxyl group 0 — H, or the ketone 
group C=0. To a fair approximation these groups contribute 
the same dipole moment in whatever molecule they occur. The 
total moment of such a molecule is then the vector sum of the 
moments of all groups present in the molecule. Clearly this rule 
cannot hold exactly for several reasons; one of these is the in- 
fluence of the interaction between the various groups on their 
individual moments. 

The second rule refers to the optical polarizabihty of large 
molecules. As mentioned above, optical polarization is due to 
displacement of electrons on the assumption that all nuclei are 
kept in fixed positions. Clearly the largest contributions will 
arise from the electrons with the smallest binding energies, i.e. 
from valency electrons. The behaviour of these electrons is 
different in a compound than it is in the atomic state of the atoms 
contained in the compound. Thus the optical polarizability of 
CO 2 is different from the sum of the polarizabilities of a carbon 
and two oxygen atoms. On the other hand, the distribution of 
electrons in certain molecular groups is to a fair approximation 
independent of the type of molecule in which they occur, as in 
the above case of the dipole moment of such groups. In applying 
this rule it should be remembered that frequently the polariza- 
bility of such a group depends on direction; for the ketone group 
C=0, for instance, the polarizability in the direction of the 
line connecting the two atoms would be different from that in a 
direction perpendicular to it. 
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Classification of dielectrics 

The division of polarization into three types which was intro- 
duced above is of a yery general nature. It leads at once to the 
following division of dielectric materials into three classes: 

(i) Non-polar substances showing optical polarization only. 

(ii) Polar substances having optical as well as infra-red 
polarization. 

(iii) Dipolar substances which in addition show also polariza- 
tion due to dipolar orientation. 

In the first class of materials, non-polar ones, an electrical 
field produces elastic displacement of electrons only. This class 
contains all dielectrics consisting of a single type of atom, 
whether they form gases, liquids, or solids. Examples are 
diamond, oxygen (solid, liquid, or as vapour), the inert gases, 
and many others. The dielectrics of the second class, polar 
materials, are capable of infra-red polarization as well as optical 
polarization. Substances of this type may contain dipolar 
groups of atoms, but these groups must show only elastic dis- 
placement. If, on the other hand, there are several equilibrium 
positions for the dipole, then the substance belongs to the third 
class of materials. 

The second class contains, first of all, substances consisting 
of molecules whose total dipole moment vanishes, though they 
contain dipolar groups of atoms. Examples are CDg, paraffins 
CH3 — (CHa)^^ — CH3 (cf. § 1 5 ), benzene CgHg, carbon tetrachloride 
CCI4, and many others in the solid, liquid, or vapour phase. In 
most of these substances the infra-red polarizability is only a 
small fraction of the optical polarizability. From a practical 
point of view, therefore, their behaviour is very similar to that of 
non-polar substances. At sufficiently low temperatures many 
solids consisting of dipolar molecules fall into the same category 
because the dipoles in these solids freeze in, i.e. the thermal 
energy is insufficient to turn them in a reasonable time into 
other equilibrium positions. The most representative substances 
of the polar type are, however, ionic crystals which may show 
very large infra-red polarizability. Examples are rock-salt 
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NaOl, other alkali halide crystals, TiOg crystals, and most otlier 
crystals of salts. In contrast to molecular lattices, which have a 
whole molecule at each lattice-point, ionic lattices contain one 
ion at each lattice -point. Thus rock-salt, for instance, forms a 
simple cubic lattice of and Cl“ ions (cf. § 18). In contrast 
to other dielectrics, most salts on melting become (ionic) con- 
ductors. 

All materials consisting of dipolar molecules belong to the 
third class, except in the very low temperature range men- 
tioned above where the dipoles are frozen in. Often in solids 
this process starts just below a critical temperature at which the 
substance undergoes a j)hase transition. In nearly fill of these 
materials tlie turning of a dipole into another equilibrium posi- 
tion is connected witli a turning of the whole molecule (e.g. for 
ketones, § 17). Exce[)tions are ice and some other crystals where 
a turning of the direction of a dipole may be obtained by a 
transfer of a ion from one equilibrium position to another. 

15. Non-polar substances 

The simplest type of dielectric substances show elastic dis- 
placement of electrons only, and in the classification of § 14 are 
designated as non-polar. For such substances the lowest fre- 
quency vq at which appreciable absorption occurs is usually in 
the visible or int he ultra-violet region. For all frequencies v 
which are less than by a sufficient amount, the dielectric con- 
stant should be independent of frequency. Thus for v the 
dielectric constant € should be equal to the static dielectric 
constant and should satisfy the Maxwell relation e ~ That is, 

(15.1) 

should hold between the static dielectric constant and the refrac- 
tive index n at frequencies v v^, Whitehead and Hackett [ W6] 
have recently checked this relation on diamond and found that 
it holds within the range of accuracy of the measurements. They 
measured the dielectric constant at frequencies between 500 
and 3,000 cycles per second and obtained the value S'OSiO'OS. 
The refractive index n was obtained by extrapolating measure- 
ments at various wave-lengths in the optical region to long 
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waves where n is independent of wave-length. The accuracy of 
this procedure can be seen from Fig. 19 and leads to the value 
z=z 5-66. This agreement means that absorption of electric 
waves at frequenciej^ between 3,000 cycles per second and optical 
frequencies must be so weak that (cf. 2.18) (2/7?) J e^dvjv < 0-03 
if integrated over this range. Actually diamond absorbs at 
infra-red frequencies, but the reason for this absorption has not 



Fig. 19. Dependence of the refractive index n of diamond on the wave num- 
ber according to measurements collected by Whitehead and Hackott [W6]. 

yet been found. The above condition, by the way, need not 
imply that eg is small provided the absorption bands are suffi- 
ciently narrow. 

Gases of non-polar molecules fulfil fairly closely the conditions 
required for the vahdity of the Clausius-Mossotti formula, 
(cf. §§ 6 and 8 and Appendix A 3). These conditions require 
(i) elastic displacement only, (ii) absence of non-dipolar (short- 
range) interaction between molecules, (iii) isotropy of the 
polarizability of a molecule, (iv) isotropy or cubic symmetry in 
the arrangement of the molecules. Condition (i) is fulfilled for all 
non-polar molecules; (ii) holds so long as the distance between 
molecules is sufficiently large; (iii) holds for spherical molecules 
only, and (iv) holds always in gases. Hence deviations may be 
expected for non-spherical molecules (condition (iii) violated) 
and at very high pressure (condition (ii) violated). All conditions 
are thus fulfilled for the rare gases. The main check of the 
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Clausius-Mossotti formula consists in measurement of the 
dielectric constant of a gas as a function of its density. For, 
since the number iVJ) of molecules per c.c. is obtained with the 
help of the Avogadro number (6*02 x 

iy;, = £6-02x102® (15.2) 

(d = density, W ~ molecular weight). 

It follows with 6.34 and 6.32 that 


e ^+2 


2-52x1024.^ a. 

It 


Often the quantity 

]) 2-52 X 1024 a: 


(15.3) 


(15.4) 


is defined as molecular polarizability. Using 15.3 and the Max- 
well relation (cf. 15.1), 


W€-l 
d e+2 


(15.5) 


should be a constant independent of density, temperature, or 
frequency. For the rare gases the following values for p have 
been found: 

He Ne Ar Kr X 

p = 0-5 1-0 4-2 6-3 10 c.c. 

The increase of polarizability with atomic weight is mainly due 
to the increase in the number of electrons per atom. Also, 
relation 15.1 has been found to hold very accurately. 

Agreement of equation 15.5 with experiment holds for many 
other non-polar substances over a very wide range of densities. 
Thus we notice with van Vleck f V3] that for Og the right-hand 
side of 15.5 has the magnitude 3*869 for the gas and an almost 
identical value 3*878 for the liquid despite the fact that the den- 
sities differ by a factor of over a thousand. Similar agreement is 
found in many other cases, and since the Maxwell relation holds, 
€ in 15.5 can be replaced by n^. In nitrogen, for instance, it was 
found from measurements of the refractive index that the right- 
hand side of 15.5 varied by less than 1 per cent, in a range of 
pressui’cs between 1 and 2,000 atmospheres. 
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In judging the significance of such a striking agreement it 
should be noted that in spite of tlie large variation in the magni- 
tude of €—1, the largest values of this quantity are usually 
smaller than unity. In this case, even with the assumption of 
complete absence of interaction between molecules, there is 
fair agreement with experiment. This assumption (cf. 6.14, 
using 15.4 and 15.2) leads to 

(15.6) 


In order to compare 15.5 and 1 5.(i we note that if we introduce 
the rnoleculai’ volume v given by 



(15.7) 


equation 1 5.5 is ecjnivalent to 


Thus 



(15.8) 

(16.9) 


The first term in this series is identical with 15.6. The higher 
terms which are characteristic of the Clausius-Mossotti formula 
are usually small. Thus in the above-mentioned case of nitrogen 
pjv is only slightly larger than 0*1 at 2,000 atmospheres. Com- 
pared with equation 15.6 the Clausius-Mossotti formula intro- 
duces a correction of only about 10 per cent. 

Kirkwood has investigated the influence of the anisotropy 

of the polarizability of molecules and found that 15.9 has to be 
replaced by 


where a is a quantity measuring the anisotropy of the polariza- 
bility. The effect of this and other corrections on the magnitude 
of € are usually small, since they affect the small second-order 
term only. 
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There are a great number of materials which according to the 
classification of § 14 are polar, but which, for all practical pur- 
poses behave like non-polar substances. They comprise mole- 
cules which contain a number of polar groups in such a way that 
their resultant dipole moment vanishes. The infra-red polariza- 
bility of these molecules is very often so small that its contri- 
bution to the dielectric constant is negligible. Important 
examples are the paraffins, and in the following we shall show 
why the dipole moment of a paraffin molecule vanishes. 

A molecule of a normal paraffin CH3 — (CHg)^^ — CH3 consists 
of a chain of CH 2 groups with a CH3 group at each end . The bond 
angle C — C — C is very nearly equal to the tetrahedral angle 
6 “ 2 cos~^(l/V 3 ) 109 °. In solids the chain forms a plane 

zigzag (e.g. § lf 5 , Fig. 28 ), but in liquids and gases this is not the 
case in general. Here rotation around the ( ■ — C bonds as axes 
may occur, but it will be assumed that this does not alter the 
bond angle. Each CH2 and CH3 grou]) has a di})olo moment 
whose direction will be discussed on the assumption that all 
four bonds of a C-atom make angles of 2 cos~^(l/V 3 ) with each 
other, t Thus, if a C^atoni is considered to occupy the centre of a 
regular tetrahedron, the bonds are directed towards the vertices. 
Now, according to § 14 , the C — C bonds should have no dipole 
moment in contrast to C — H bonds which will have a moment. 
The resultant moment of a CHg group is then always in the 
C — C — C plane (cf. Fig. 20) bisecting the C — C — C angle, where- 
as the CH3 moment has the direction of the adjoining C — C bond. 
To show this we consider a cube with a C-atom at its centre 
(Fig. 21). Then the four bond directions point towards four of 
its corners which are not neighbours. These corner points form 
the vertices of a regular tetrahedron. Thus, if /x is the C — H 
moment, the CH2 moment is 2/xcos^/2 — 2 jLt/V 3 . This may be 
considered as composed of two dipoles of strength /x in the direc- 
tions of the C — C bonds (dotted arrows in Fig. 20). Furthermore, 
the moment of a CH3 end-group is also equivalent to a dipole /x 
in the direction of its adjoining C — C bond. For the resultant of 

t It should be noted that in Fig. 28 the distance refers to CHg groups as 
found from X-ray measurements and not to the C — C distance. 

4980.11 J 
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any pair of H-atoms of the CHg group is a dipole of strength 
2/X/V3 in the direction perpendicular to a surface of the cube, 
as has just been shown (this assumes that C — H dipoles in a 
CHg and in a CHg group are equal). The resultant thus forms an 
angle cos”^ 1/V3 wdth the diagonal in which the residual dipole 
lies, and thus contributes (1/V3) (2/x/V3) = 2/x/3 to the CHg 
dipole. Now three CHg groups can be selected from CHg, but 
then each CH dipole is counted twice. Thus the moment of the 


I 

I 



Fm. 20. Fig 21. 

Fig. 20. Tho full ai row represents the dipole of a CHg group seen to lie in the 
C — C — C piano of a paraflin molecule. It can be decomposed into the two 
dotted arrows in the C — C directions. 

Fig. 21. A C-atom of a paraffin molecule with its four neighbours, two C- 
aiid two H-atoms. The former are indicated by rings, the latter by dots ; the 
arrow indicates the dipole of the group. The lines joining the central atoms 
with its neighboiu’s indicate the directions of the latter rather than their 
distance ; actually the C — H distance is smaller than the C — C distance. 

CHg group is -|(2/i/3) = p,. It follows that the paraffin chain 
can be considered as composed of C — C rods with a dipole /x 
on each end. These end dipoles point in opposite directions, so 
that the dipole moment of each rod, and hence of the whole 
chain, vanishes exactly. 

As an example we mention pentane CgHjg which at 30® C. 
under a pressure of 1 atmosphere has a dielectric constant 
€ = 1*82. The refractive index isn = 1*36, i.e. = 1*85 which 
shows that the Maxwell relation 15.1 is fulfilled — probably 
within experimental error. Hence the contribution of the 
infra-red polarization must bo very small. The following table 



IV, § 15 


NON-POLAR SUBSTANCES 


116 


of measurements by Danforth [Di] shows that the Clausius- 
Mossotti formula also holds to a good approximation. 


Pressure 

€ 

(€- 1 )/( c + 2 ) 

1 jv 

P ^ t’(c-l)/(€+2) 

1 atm. 

12000 atm. 

1-82 

1 2-38 

0-216 

0-.308 


0-356 

0-339 


Again it can be seen that (cf. 15.9) 1 so that the zero 

approximation of 15.9 accounts for more than two-thirds of 
e — 1 . More accurate comparison with equation 15.10 shows that 
a value a —0*05 is required to account for the experimental 
data. This shows that the next approximation — i.e. the first one 
which is significant for the Clausius-Mossotti formula — deviates 
by only 5 per cent, from 15.9. 

For an example exhibiting a larger contribution of infra-red 
polarization we turn to the measurements of Michels and 
Hamers \^M2\ and Michels and Kleerekoper [MS] on OOg. This 
molecule is rectilinear and has, therefore, no dipole moment. 
Dielectric measurements give the value 7*5 for ^(c— l)/(6-f-2); 
from refractive index measurements, on the other hand, the 
value 6*7 is found for v{n^— 1 )l{n^-j-2). Hence about 10 per cent, 
of the static polarization is due to infra-red contributions in this 
case. 

Relatively large infra-red contributions are obtained in dipolar 
solids at low temperatures at which dipoles are frozen in, so that 
these substances behave in a similar manner to non -polar 
materials (temperature-independent dielectric constant in a 
certain range). As an example, the case of ketones will be dis- 
cussed in more detail in § 17. 


16. Dipolar substances 

Oases and dilute solutions 

Provided the pressure is not too high, the static dielectric 
constant e^ofa dipolar gas satisfies the relation e^—1 ^ 1. There- 
fore should satisfy equation 6.15, which, with the help of 15.2 
can also be written as 


= t6.02xl0-4i|, 


(16.1) 
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or, inserting numerical values and expressing fi,,, in Debye units 
(cf. § 6) and T in degree Kelvin, 

1-83x1 (10.2) 

This equation, according to § G, has been derived on the assump- 
tion that the interaction between dipoles is negligible. This 



Fig. 22. Temperature-dependence of the dielectric constant of 
various gases according to Sanger [SI]. 


involves an error of the order of 1)^, as can be seen by com- 
parison with formulae where the interaction has not been 
neglected (e.g. 6.36). 

Since no assumptions about the interaction have to be made for 
the derivation of 16.1 or 16.2, measurement of the temperature- 
dependence of the static dielectric constant of gases should make 
possible a reliable determination of the dipole moment of free 
molecules. As an example Fig. 22 shows the results of measure- 
ments of the dielectric constants of CH4, CH3CI, CHgClg, CHCI3, 
and CCI4 gas plotted against the inverse temperature IjT (from 
Sanger, SI). Since the highest value of 1 is only about 10“^, 
the above equations should hold very well. As required, the 
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experimental points lie on a straight line whose slope allows us 
to calculate the dipole moment with the help of 16.1 or 16.2 if 
the density is also measured. 

The results are also of interest in the field of stereochemistry. 
Thus the fact that the dipole moments of methane (CH4) and of 
carbon tetrachloride {CCI4) both vanish (because is indepen- 
dent of T) is evidence in favour of a molecular structure in which 
the four hydrogen (or chlorine) atoms form the corners of a 
regular tetrahedron with the carbon atom at the centre. 

Figure 22 also permits a determination which is obtained 
by extrapolating the straight lines to l/T — 0. The value 
obtained should be slightly larger than the square of the optical 
refractive index n. But, as already mentioned in § 15, the differ- 
ence €^—n^ is very small for most molecules; in gases it is prob- 
ably too small to be measurable using present experimental 
techniques. 

There should be no appreciable frequency dependence of the 
dielectric constant, and hence no loss, up to the absorption 
frequencies of the molecules, which often lie in the infra-red 
region. For a number of molecules absorption begins at the still 
longer wave-lengths of the ultra-short electric wave region. 
Such gases will then show resonance absorption in this region 
following the laws derived in § 13. Relevant experiments will be 
discussed in detail in the book by Jackson and Saxton [J 4 ], 
Here we merely wish to discuss briefly the shape of the absorption 
spectrum of ammonia which, following Cleeton and Williams 
\C 2 ], has been investigated in the centimetre region in great 
detail by Bleaney and Penrose [B 3 ]. Ammonia, NH3, forms a 
pyramidal molecule and the absorption is connected with the 
swinging of the nitrogen atom through the plane of the three 
hydrogen atoms. The spectrum shows considerable structure, 
which is resolved at pressures below 5 cm. Hg, as shown in Fig. 23. 
The shape of a single absorption line should be given by equation 
13.11, provided the concept of a single relaxation time r can be 
applied. Assuming l/r to be proportional to the gas pressure, 
the shape of the absorption spectrum at a pressure of 10 cm. of 
mercury has been calculated by Bleaney and Penrose with the 
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help of 13.11 using data from measurements at 0*5 mm. Figure 
24 shows excellent agreement with the measured absorption. 
At a pressure of 60 cm., however, agreement is no longer so 
good. 

We shall now turn to dilute solutions of dipolar molecules in 
non-polar substances. The concentration will be assumed to be 



Fic. 23. The ammonia absorption spectrmn at 1, 2, 5, and 10 cm. pressure 
according to Bloaney and Penrose {B3]. 


sufficiently small so that interaction between dipoles may be 
neglected. It then follows from equation 6.19 that the static 
dielectric constant satisfies an equation similar to those of 
gases, but with an effective dipole moment which depends on the 
nature of the solvent and on the structure of the dipolar mole- 
cule. This effective dipole moment can be calculated in a simple 
way only when the molecule can be approximated by a point 
dipole at the centre of a sphere with refractive index n. In this 
case equation 6.24 should hold. In general, we may put (eg = 
dielectric constant of solvent; Nq = number of dipole molecules 
per unit volume) 


ZkT \ 3 / 


{1-y)^ 


■8 


■00 


(16.3) 
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where y is unknown except in the case of spherical molecules for 
which, according to 6.24, 

y = ^ for spherical molecules. (16.4) 

(2€o+^^)(fo+2) 

It follows that \y\ is usually much smaller than unity, but, in 
view of the factor €q— may be either positive or negative. 



Fig. 24. Calcxilated and measured absorption spectrum of ammonia at 
10 cm. pressure according to Bleaney and Penrose [B5]. 

The main conclusion to be drawn is that by using experimental 
values for the temperature-dependence of it is not possible to 
calculate the dipole moment of the free molecule but only the 
quantity /x^(l— y) which differs from though usually only by 
a small amount. 

Besides the solvent effect, one of the main differences between 
gases and dilute solutions is their behaviour in alternating fields. 
The interaction between the dipolar molecules and the mole- 
cules of the solvent leads to an energy loss whose maximum at 
most temperatures falls within the frequency region of electric 
waves. This affords a new way of obtaining the important 
quantity without measuring the static dielectric constant 
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€/, for if the dielectric loss, and hence the imaginary part 
of the complex dielectric constant e(co) (of. § 2) is known as a 
function of frequency, €qo is obtained from equation 2.18, 
which may also be written as 

Cg— ^ J e^{oi)d{logoj), (16.5) 

and which holds independently of the special shape of the func- 
tion € 2 ( 0 )). The integration has to be carried out over the whole 
range of frequencies in which € 2 ( 0 )) has appreciable values. The 
present method can be expected to yield more accurate values for 
€^—€ao than measurement of the static dielectric constant 
because is nearly equal to the dielectric constant €q of the 
solvent and is small compared with in the case of dilute 

solutions. Very high accuracy in the measurement of would 
therefore be required to give reasonably accurate values of 

The investigation of dilute solutions permits a study of the 
behaviour of single molecules independent of dipolar interaction; 
and the results obtained will be of interest in attempts to under- 
stand the properties of less dilute solutions, or of pure dipolar 
substances. To obtain full information the dielectric loss should 
be measured over a large range of frequencies at different 
temperatures. Absorption can usually be described by the 
.Debye formulae (§ 10) or by a generalization (§ 12). Frequency 
measurements will establish whether the Debye equations hold; 
if they do, they will lead to a value for the relaxation time, but 
otherwise they will permit an estimate of the width of a band of 
relaxation times. Temperature measurements will then yield 
the temperature -dependence of the relaxation time — or of the 
width and position of the band of relaxation times. No complete 
system of measurements is available at present. 

As an example, we shall now discuss the experiments by Jack- 
son and Powles \JS] who have measured at one temperature 
(19^^ C.) the frequency dependence of dielectric loss of dilute 
solutions of dipolar molecules in benzene and in paraffin. In 
Fig. 25 the experimental values of tan (<^ = loss angle) for 
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benzophenone solutions (1 gm./lOO c.c. at 19^^ C,) are plotted 
against frequency. For the solution in benzene these values 
lie on a curve which represents the Debye formula. According 
to 10.17, 10.22, and 10.23, making use of 10.29, this formula can 
be written in terms of the angular frequency at which the 
maximum loss-angle cf)^ occurs: 

tan^ = (10.6) 

whore t&n<f>„^ = < 1. (16.7) 



Fig. 25. Loss angle of dilute solutions of benzophenone in benzene (inaximuni 
near cycles) and paraffin ; experimental values according to Jackson and 
Powles The full line represents a Debye curve (equation 16.6, or 16.8 

with p — 1); the dotted line represents equation 16.8 with — 5. 


From 10.22 the relaxation time is then given by t = l/oi^. 

In contrast to benzene the solution of benzophenone in 
paraffin leads to a broader absorption curve than could be satis- 
fied by a Debye curve. According to § 12 this can be described 
by introducing a whole range of relaxation times. If, in par- 
ticular, we introduce a model where the relaxation times are 
distributed between the two values tq and according to the 
distribution function 12.13, then tanc;^ can be obtained after 
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division by €i(a>) e., (of. 10.29) from 12.16 or from 12.18 and 

12.19. Thus with the use of 12.17, 


where 

and 


tan (/) = tan <f)^ 


tan“^jS— tan“^( l/)3) 


tan — 


1 €,,—€00 tan“ij8~tan-i(l/j3) 

2 log^ 

P = VKK) > 1- 


(16.8) 

(16.9) 

(16.10) 


For j8 ~ 1, equations 16.8 and 16.9 are identical with the Debye 
equations 16.6 and 16.7. For jS > 1, however, the curves are 
broader and flatter. Figure 25 shows that the curve with j3 = 5 
satisfies the experimental points fairly well; appreciable devia- 
tions occur only at frequencies very far from a>^ where tan^ 
is relatively small. This is rather to be expected, as follows from 
the discussion in § 12 following equation 12.19, because the 
distribution function 12.13 can be considered to represent 
position and width of the distribution of relaxation times but 
not finer details. 

It will be of interest to see how the time which in the 
Debye case represents the relaxation time is now connected 
with the two limits tq and r^, between which the relaxation 
times are distributed. According to 12.17, 

= VK'^i) == ^"^0 = (16.11) 

Thus, in the present case l/co^^ = 5to ti/5. 

Since no experiments at other temperatures are available at 
present, it is not yet possible to make use of the temperature- 
dependence of the relaxation times to discuss the relaxation 
mechanism. On the assumption of a model in which the relaxa- 
tion time is determined by the frequency of jumps of a dipole 
over potential barriers (cf. § 11), the factor allows us to calcu- 
late the maximum difference Vq in height of these potential 
barriers on the assumption that they are equally distributed 
between the values Hq and Hq+Vq (cf. 12.7). 

Thus from 12.17, using 16.10, 


vJkT = log(Ti/To) = 21og^. 
In our case, therefore, vJkT — log 25 ci 3. 


(16,12) 
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It would be of interest to investigate the variation with tem- 
perature of the ratio as well as that of (or tq) itself. 
Measurements at various temperatures were carried out by 
Whiffen and Thompson \ W2\ though unfortunately, up to the 
time of writing, the range of frequencies has been rather re- 
stricted. The dielectric loss in solutions of chloroform in hep- 
tane seems to follow a Debye curve as indicated in Fig. 26, but it 



Fig. 26. Loss angle of a dilute solution of chloroform in heptane according to 
Whiffen and Thompson [TV2]; the curve represents the Debye formula 16.6. 


must be realized that in the available frequency range there is 
only a variation of a factor of three in tan^. The temperature 
was varied in a range between — 70° C. and 80° C. It was found 
that the logarithm of the relaxation times is proportional to l/T. 
A similar result was obtained for solutes, as shown in Fig. 27. 
To test the validity of the Debye model (§ 11), according to which 
the temperature variation of r should be the same as that of the 
viscosity rj of the solvent, log?; is also plotted in the same figure. 
Putting (16.13) 


Debye ’s model according to 1 1 .32 requires == //^. The experi- 

mental values are as follows: 


cx-Bromo-naphthalene 
Methyl benzoate 
Camphor 
Chloroform 


in k,cal. 

. 1-8 
. 1-8 
. 1*7 

• 1*5 


On the other hand, for heptane = 2-0, which thus is larger 
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than for all the above molecules. This indicates that in the 
above cases the Debye model probably does not apply exactly. 



2-5 3-0 3-5 4-0 4-5 5*0 5*5 

i/txio^ 


Fig. 27. Temperature-dopendence of the relaxation time r of various solutions 
in heptane, and of the viscosity t} of heptane, according to Whiffen and 

Thompson IW2], 

In all experiments on dilute solutions the value of €oo can 
be obtained from dielectric loss measurements by carrying out 
the integration indicated in equation 16.5. It should be remem- 
bered that for dilute solutions both quantities 6i(a>), are nearly 
equal to the dielectric constant of the solvent, so that (cf. 
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2.5) €^{cji)) can be replaced by cotan^. Prom -Coo, then, the 
dipole moment can be obtained according to 16.3 if y can be 
neglected. This method was first used by Sillars [aS 6] in the case 
of dilute solutions of dipolar long-chain molecules in solid 
paraffin -wax. Although equation 16.4 for y is certainly not 
applicable in this case, it nevertheless indicates that y should be 
particularly small if the square of the refractive index of the 
pure solute is nearly equal to the dielectric constant of the sol- 
vent. For long-cliain molecules this should be very well fulfilled. 

- Chain axis 
CH3 

Fia. 28. ParafYiii dmin. Tlie indicated distances refer to X-ray scattering 
centres of (/Hg gronpH ac(;ording to Muller [M6'] ; they are sligVitly different 
from C — 0 distances. 

The investigation of the solid solutions mentioned above 
are of great interest for the elucidation of the mechanism of 
dielectric loss and will, therefore, be discussed in greater detail. 
We shall start with a discussion of the structure of substances 
built up of long-chain molecules. Most of these structures can 
be derived from that of paraffins, which was investigated by A. 
Muller \^M6\ A paraffin molecule, in a crystal, consists of a 
plane zigzag with a CHg group at each corner and a CHg group at 
each end (cf. Pig. 28). The distance between neighbouring 
CHg groups is about 2 A (1 A ™ 10~®cm.), while the distance of 
their projections to the chain axis is about 1*25 A. In a paraffin 
crystal the chains are arranged in layers whose thickness is 
approximately equal to the chain length. Within such a layer 
the molecules form rectangular cells with side length a, 6, c, 
where a 5 A, 6 7*5 A, and c is slightly larger tlian the chain 

axis and is parallel to it. Figure 29 shows how the chain planes 
cross the a-b plane which is perpendicular to the chain axes. It 
is of importance to notice that in the subsequent layer the whole 
arrangement is shifted by about 1 A in the direction of the 6-axis. 
This is indicated in Fig. 29 by the dotted lines. 

A paraffin molecule carries no dipole moment, as has already 




126 


APPLICATIONS 


IV, § 16 


been mentioned in § 15. Now suppose that some dipolar long- 
chain molecules, c.g. esters or ketones, are dissolved in a paraffin 
crystal. This means that some of the paraffin molecules are 
replaced by dipolar long-chain molecules. These latter can be 
usually derived from paraffin molecules. A ketone molecule, for 
instance, is obtained from a paraffin molecule by replacing one 
CHg group by a CO group. This latter carries a dipole moment 
whose direction is perpendicular to the chain axis and probably 



Fig. 29. Fig. 30. 

Fig. 29. Section of paraffin crystal perpendicular to chain axis indicating the 
chain planes, according to Muller [MS'], The dotted lines indicate the section 
of chain planes in a subsequent layer. 

Fig. 30. The two equilibrium positions of a ketone molecule — indicating its 
dipole directions — ^replacing a paraffin molecule of longer chain length. The 
size of the rectangular cell is the same as in Fig. 29. 

lies in the chain plane. Also assume that the length of the ketone 
molecule is shorter than that of a paraffin molecule so that it 
can easily replace a paraffin molecule. There should then be two 
possible directions for the dipole which differ by 180°, as indi- 
cated in Fig. 30. For, if the ketone molecule is shorter than the 
paraffin molecule, the former should still fit into the structure if 
it is turned through 180° and shifted by the length of one link 
along the chain axis (cf. Fig. 31). In a dilute solid solution such 
a dipolar molecule has thus two equilibrium positions with 
opposite dipole directions. From the similarity of the two 
positions it can be assumed that the energy of the molecule is 
the same in both. 

In general, a molecule will oscillate about either of the two 
equilibrium positions. Occasionally, due to a thermal fluctua- 
tion, it will acquire sufficient energy to turn into the other 
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equilibrium position. Such asubstance thus may well bedescribed 
by the high -temperature model of dipolar solids discussed in 
§11, where a dipole was supposed to have two equilibrium posi- 
tions separated by a potential barrier; the minimum energy 
required to lift it over the barrier was denoted by H. The energy 
loss of such a substance is described by the Debye formulae 
(e.g. 16.6 and § 10) and the relaxation time is given by equation 
11.3 (cf. also 11.8). 

Loss measurements on dilute solutions of the type described 
above have been carried out by W. Jackson [Ji, 2], Sillars [S6], 



Fio. 31. Result of turning a ketone molecule by 180® and shifting it by one 
chain length. The arrow indicates the dipolar CO group. 

and Pelmore \P2]. These authors find that the dielectric loss 
curves satisfy the Debyeequations fairly well, though notexactly. 
This slight deviation from the Debye equations may be due to 
the fact that a dipolar molecule which is shorter than a paraffin 
molecule, say by z links, has 2 possible locations (apart from the 
two dipole directions) corresponding to different positions on 
the c-axis. Its energy in these positions may not be exactly the 
same, and the same may be true of the energy H required to lift 
a molecule over the potential barrier. This would lead to a distri- 
bution of relaxation times. It should be realized that in calcu- 
lating the relaxation time t according to the Debye formula 
1 0.22 some error is introduced (cf. 16.11), though it may be negli- 
gible if the variation of logr with temperature, and not the 
absolute value of r, is considered. 

The most interesting result of these experiments is the 
measurement of the dependence of the relaxation time on the 
chain length of the dipolar molecule. If the molecule is con- 
sidered to be rigid, and if therefore each link has to be lifted over 
the barrier at the same time, the total energy required to lift 
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a whole molecule of chain length m over the barrier should 
increase in proportion to m. Thus, since according to 11.3 

logr “ constant (16.14) 

the experimental values of logr plotted against m should lie on 
a straight line. This is not the case. Actually molecules are not 
rigid but have a certain flexibility, and, if this is taken into 
account (cf. reference F5), satisfactory agreement with experi- 
ment is obtained. The influence of the flexibility on the value of 
can easily be discussed qualitatively. Suppose for a moment 
that the chain is rigid. Then the whole chain has to be lifted over 
the potential barrier at the same time. Thus if is the energy 
required to lift a single link over the barrier, the total energy 
would be equal to mH^ where m is tlie number of carbon 
atoms in the ohain. Actually, however, the chain has a certain 
flexibility. This allows the chain to be gradually lifted over the 
potential barrier so that a total energy less than is requii*ed. 
On the other hand, a twisting of the molecule requires a certain 
energy; but this energy decreases with increasing chain length 
for a given angle of twist of the two ends of the chain. Thus 
for short chains oc m, but for long chains will tend to a 
constant value independent of the chain length. A simple 
calculation [i^5] yieldsf 

07) 

= i/imQtanh— , (16.15) 

m 

where m is a constant which separates short chains (m <,m, 
little torsion) from long chains (m > m, appreciable torsion). 
Comparison with experiments shows that 1/30 e-volt, 

while m = 26, both reasonable values. To obtain t a further 
constant is required according to equation 16.14. Thus three 
experimentally determined constants are included in the final 
formula (r in seconds) 

log27rT = -60-4+i^^tanh|^. (16.16) 

Here the relaxation time t depends on two parameters, the tem- 
perature T and the chain-length m. The formula holds, however. 


t Note that tan h® = (e®-e-*)/(e*+e-*). 
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only if m is smaller than the chain length of the paraffin in which 
the dipolar chains are dissolved. Figure 32 shows that the results 
obtained from equation Hi. If) agree well with tlie experimental 
results. The expression tanh(m/26) is plotted as a function 
of m/26 and is compared with the experimental values of 
(log 27rT+50-4)T/13800 for four different chain-lengths (m “ 20, 
22, 24, 32). Two of the constants in equation 16.16 have 
been determined from the experimental values for m = 20 and 



Fig. 32. Comparison of the theoretical dependence of the relaxation time on 
chain length with experimental values, according to [EJ]. 


m = 22, while the third was determined from the temperature- 
dependence of m = 20. For m ~ 24 and 32 no more adjustable 
constants are available, which shows that formula 16.16 leads 
to a correct dependence of the relaxation time on chain lengths. 

From the numerical value —50-4 for the constant in equation 
16.14 (comparing with 16.16) it might appear that the order of 
magnitude of the constant ^ of equation 11.3 (cf. also 11.8) could 
be estimated. For assuming for the frequency of oscillation 
o)J 27 t 10^^ one finds log^ ~ -”25. Not much significance 
should be attributed to this value, however, in view of the un- 
certainties connected with the experimental determination of 
the absolute value of t; for, as mentioned above, the deviation of 
the loss curve from a true Debye curve brings with it some un- 
certainties in determining t which cannot be adjusted by apply- 
ing equation 10.22 as in the Debye case. This would probably 

4980.11 17 
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influence the value of ^ to a much greater extent than the less 
sensitive 

In all examples discussed above the concentration of dipolar 
molecules was assumed to be so small that dipolar interaction 
would be neglected. In this range, therefore, the quantity 
Cg— £oo well as dielectric loss sliould increase proportionally 
with the concentration. For further increases of concentration 
the effects of dipolar interaction should start to appear, and an 
investigation of the deviations of from equation 16.3 

would be of great theoretical interest. In particular, it would 
be of interest to see to what extent and in which range of con- 
centration the Onsager formula 6.38 can account for such devia- 
tions. This would be of particular interest for a solute which in 
the pure state does not satisfy the Onsager formula for the 
temperature range contemplated. f It will be remembered (§ 8) 
that the Onsager formula takes into account the long-range 
dipolar forces, but neglects short-range forces. Then with 
increasing concentration we should expect three stages: (i) all 
interaction can be neglected, leading to equation 16.3 for (e^— ; 
(ii) long-range dipolar forces have to be considered, but short- 
range forces can still be neglected, leading to the Onsager formula 
6.38; (iii) short-range forces have to be considered as well. This 
leads to the Kirkwood formula 8.14. 

Of very great importance would be the investigation of the 
dependence on concentration of the shape of the dielectric loss- 
frequency curve. The deviations from a Debye curve would 
probably increase with increasing concentration owing to the 
interaction between dipoles. As yet no theoretical investiga- 
tions on this subject have been carried out. 


17. Dipolar solids and liquids 

Survey 

From the developments of §§ 7 and 8 it follows that no formula 
exists which represents in a simple manner the dielectric 

t Measurements on the dependence of €g on concentration have been 
carried out, of course. They do not extend, however, to the region of very low 
concentration where €g—€cx> is best obtained from loss measurements with the 
help of 16.5. 
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behaviour of dipolar solids and liquids. Qualitatively, however, 
most of these substances behave in a very similar way. At very 
low temperatures in the solid state all dipoles are frozen in, and 
if the solid has not become permanently polarized the dipoles will 
not contribute to the dielectric constant. In this region (T 0), 
therefore, the dielectric constant e^(0) should be nearly inde- 
pendent of temperature, and its magnitude should be approxi- 
mately equal to that of the dielectric constant e^oiT) at a higher 
temperature but at a frequency which is so high that the dipoles 
do not take part in the polarization. This equivalence holds 
only, of course, if the volume of the substance is the same at both 
temperatures, and it includes the assumption that upon turning 
a molecule from one equilibrium position into another its con- 
tribution to the high-frequency dielectric constant is not 
altei'ed. Then 

€g{0) €^(T) = constant volume. (17.1) 

In this formula it has been indicated that is composed of 
the square of the optical refractive index n (cf. § 15), and of a 
term Ae due to elastic displacement of nuclei or of dipoles. The 
corresponding polarization was called infra-red polarization in 
§ 14. As pointed out in § 15, this type of polarization is frequently 
negligibly small in non-dipolar substances. In dipolar substances 
Ac usually has small but noticeable values, mainly due to a dis- 
placement of the equilibrium directions of the dipoles by the 
field. The frequency of the rotationary oscillations which the 
dipoles carry out about their equilibrium positions is normally 
in the far infra-red region, but there are indications that in large 
molecules it may approach the centimetre region (cf. reference 
F9), Such substances would thus show resonance absorption 
(§ 13) at ultra -high frequencies. Substances which are expected 
to show such absorption have been discussed by Szigeti [5^17] 
and a further theoretical investigation on the expected fre- 
quency-dependence of the absorption has been carried out by 
Huby [7?3] leading to a generalization of equation 13.11. In 
this connexion it should be mentioned that Girard and Abadie 
[03] found that normal liquid long-chain alcohols absorb in the 
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centimetre region besides showing Debye absorption at longer 
wave-lengths. It was suggested by Magat \M1] that this is an 
example of resonance absorption. Further experiments would 
be required to establish the validity of this interpretation. 

Let us now turn to the regions of higher temperatures and 
lower frequencies wliere dipole orientation plays an important 
role. Within creasing temperature an increasing number of 



Fig. 33. Temperature -dependence of the dielectric constant of solid hydrogen 
sulphide according to Smyth and Hitchcock [<5,9] showing an order-disorder 
transition near 103° abs. A further transition and melting are of little influence. 

dipoles can turn into other equilibrium directions, as pointed 
out in § 8. The static dielectric constant will thus increase 
with temperature, at first according to the relation 

e^~~e^az (17.2) 

as follows from 8.17 (F(0) is a constant energy discussed in § 8). 
At higher temperatures the rate of rise of gets increasingly 
larger until at the temperature the solid carries out an order- 
disorder transition beyond which (T > T^) the dielectric con- 
stant €s falls with increasing temperature. The order-disorder 
transition may occur in the solid state; a number of examples are 
shown in Pigs. 33 and 34. In this case there is usually no appre- 
ciable change of at melting. In other substances melting occurs 
before the order-disorder transition is completed. In these 
cases the decrease of starts at the melting-point. 
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As shown in the discussion of entropy in § 3, a positive value 
of dcJdT indicates that the degree of order is decreased when 
the field is applied, while a negative value of deJdT indicates the 
opposite, namely, that the degree of order is increased by the 
application of the field. It follows that in a substance in which 
the dipole directions are completely disordered deJdT must be 
negative, whereas if the dipole directions are completely ordered 



Fig. 34. Temperature-dependence of the dielectric constant of solid ethylene 
cyanide according to White and Morgan [W4] showing an order-disorder 
transition near — 40® C. without a discontinuity at molting. 

deJdT must be positive. Tims from purely thermodynamical 
j easoning it follows that the transition from positive^ to negative 
values of deJdT must be connected with a change in the order 
of dipole directions. A large number of examples of such transi- 
tions has been reviewed in a recent article by Smyth [S7], 
Dielectric loss due to dipole orientation is expected to occur 
at T > as Well as for T <Tc (Figs. 35 a and b give examples 
for the latter). Except very close to the latter will, of course, 
be much smaller than the former because of the much smaller 
value of € 5 — Thus (cf. Fig. 36) at any temperature we expect 
resonance loss in the far infra-red or ultra-high frequency region, 
and loss of the Debye type (but requiring a distribution of 
relaxation times) at longer wave-lengths. For the latter, no 
satisfactory theoretical treatment of the shape of dielectric 






- 150 - 100 ~ 50 0 + 20*C 

Fig. 35. (a) Dielectrio constant and .(5) dielectric, loss oc of di-isopropyl 
ketone according to Schallamach [6*^]; — O — af ft frequency of 1-2 Mc./s., 
— 3 — at 4-4 Mc./s., — • — at 20 Mc./s. Melting occurs near — 73° C. and there is 
no previous transition. 
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loss -frequency curves exists at present. Compared with dilute 
solutions the difficulty consists in accounting for the influence 



Fia. 36. Frequency -depen den CO of dielectric loss, above and 

below the transition temperature (schematically). cool2n represents the 
lowest resonance frequency which is usually in the far infra-red or in the ultra- 
high electric frequency region. 



Fig. 37. Frequency-dependence of loss angle according to Hartshorn, Megson, 
and Rushton [HJ] for various materials. 

(1) Phenolic resins. (2) Benzyl alcohol resins. (3) Rubber- sulphiu* compoimds 
(Scott, McPherson, and Curtis, S5). (4) Cetyl palmitato in j)araffin wax 
(Sillars, S7). (5) Chlorinated diphenyl (Jackson, J2). (6) Represents a Debye 
curve; comparison with Fig. 14 sliows that cui’ves (l)-(5) can be represented 
by equation 12.19 with suitable values of the parameter ti/tq. 

of interaction. This, one should expect, would lead to a broaden- 
ing of the Debye curves similar to that obtained through the 
existence of a whole range of relaxation times (cf. examples in 
Fig. 37). The width of the latter ought to be larger in amorphous 
substances than in crystalline solids. It would be of interest, 
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therefore, to investigate substances which at the same tempera- 
ture exist in various states of order; e.g. amorphous substances 
which have been quenched from different temperatures, or 
substances existing at the same temperature as super-cooled 
liquids and as crystalline solids. 


Analysis of curves 

Without any knowledge of the structure of a dielectric, in- 
teresting conclusions about the behaviour of its dipoles can be 
drawn from the experimental T curves with the help of our 
general theory. For from equation 7.39 it follows that we can 
obtain the quantity 


— -~mm* 
3 k 




3e, 


T 


(17.3) 


provided that the optical refractive index nis known in the correct 
region (cf. § 15). If this is not the case, we can make use of the 
fact that the infra-red contribution Ae = is usually small 

(except in ionic crystals) compared with and will replace 
by and hence, according to 17.1, by ^^(O). Thus, if 


B{T) 


3 k 


mm’' 


(17.4) 


we find approximately 

B(T) ~ T. (17.5) 

This value 17.5 for B[T) differs from the correct one by 
y^^ g8(7’)+2e^(0)-Ae ^ 

which is of the order TAe. Using approximation 1 7.5, it is thus 
possible to obtain B{T) from measurement of the static dielec- 
tric constant at various temperatures. Theoretically B{T), 
according to 17.4, is proportional to mm*, and the value given 
by 17.5 is thus constant (subject to the small correction 17.6) 
independent of temperature, if Onsager’s formula holds, i.e. if the 
dipole moment of a molecule can be considered as a constant and 
if there are no short-range forces leading to mutual orientation 
of dipoles* Since for T ->co all interaction can be neglected 
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and the dipole orientation is completely random, a function 
B{T) which increases with decreasing temperature points to an 
increasing tendency towards parallel orientation of dipoles. 
Similarly if B{T) decreases with decreasing temperature the 
dipoles tend towards anti-parallel orientation. It should be 
pointed out again (cf. §§ 7 and 8) that this orientation is due to 
non-dipolar short-range forces. Long-range dipolar forces too 
lead to a certain tendency for dipole orientation. But averaging 
over all angles between dipole direction and radius vector, this 
orientation cancels in an isotropic substance. Since m* is the 
average moment of a sphere if one of its units (molecules) has 
the moment m, dipolar interaction alone leads to m* = m, 
i.e. to a value of mm* which is independent of temperature if 
m is a constant. 

If the assumption is no longer made that the dipole moment 
m of a single unit is constant, then similar conclusions to those 
above can be drawn from the temperature-dependence of B{T). 
A decrease of B{T) as the temperature is lowered can now mean 
either an increasing tendency towards anti -parallel orientation 
between the dipoles of neighbouring units (molecules or unit 
cells) or a decrease of the dipole moment of a single unit, or 
vice versa if B{T) increases. 

In Figs. 38 (i) and (ii) a few examples of this analysis are 
given. We shall now show how in the case of a known structure 
the general theory can be applied. 

Water [K4, 02, K6] 

Water no doubt is one of the most important dielectrics, and 
the correct calculation of its dielectric constant must be con- 
sidered as a great success of Kirkwood’s [K4\ formula 8.5, 8.14 
which applies to dipolar liquids in general. A water molecule can 
be considered as a negative O"” ion with the two positive ions 
attached in such a way that the lines connecting the latter with 
the centre of the oxygen ion form an angle of 105°, as can be con- 
cluded from the infra-red absorption bands of HgO vapour. The 
dipole moment of the HgO molecule is a vector directed along 
the line dividing the H — 0 — H bond angle into two equal angles 
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of 52|®. In the liquid state the position of an individual molecule 
is strongly correlated with those of its neighbours, and X-ray 



Fig. 38. Examples of — T analysis, (ii) Temperature -dependence of the 
function B(T) obtained from eg(T) (cf. i) with the help of equation 17.5. 

(o) Oichloro -propane, (CHgluCCL, solid, using c^(0) ^ 2-23. (6) Tertiary- 
butyl -chloride, (CH8)3CC1, liquid, using €^(0) = 2-45. (c) Penta methyl- 
chlor-benzene, (CH 8 ) 5 C 1 , solid, using €g{0) = 2-8. Measurements: (a), 

(6) Turgewich and Smyth [7’2]; (c) White, Biggs, and Morgan [W3]. 

With decreasing temperature, (a) indicates increasing tendency to 
anti-parallel — (6) to parallel — orientation of dipoles, (c) points to 
random orientation, and to validity of Onsager’s formula. 

investigations show that the average number of nearest neigh- 
bours is close to four. Also, X-ray investigations of liquids show 
that in general the nearest neighbours of a molecule are arranged 
in a fairly ordered way corresponding to a certain type of crystal 
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structure. In contrast to solids, however, this order gradually 
disappears for more distant molecules. In water, according to 
the model of Bernal and Fowler [jB 2], the four nearest neigh- 
bours of a given molecule form a nearly regular tetrahedron with 
that molecule at its centre. The bond between neighbouring 
molecules, called the hydrogen bond by 
chemists, is supposed to be directed 
along the O — H bond of one molecule 
towards the oxygen ion of the other. 

The arrangement of a molecule and its 
four neighbours may be represented 
schematically as in the diagram, (Fig. 

39) in which the hydrogen bonds are 
indicated by dotted lines. 

To simplify the following considera- 
tions we shall assume the bond angle 
to be 2cos“"^ 1/V3 cr! 109° instead of 
105°. This should make little differ- 
ence to the final result. It enables us to 
consider the tetrahedron formed by the 
four neighbours of a given molecule as regular, with that mole- 
cule at the centre. A simple way to picture the arrangement of 
molecules and their dipoles is to start with a cube with four 
of its vertices forming a regular tetrahedron, and the centre 
occupied by HgO molecules (Fig. 40). The four bonds emanating 
from the centre, namely two — O — H* • •0<( bonds and two 
)0* • -H — O — bonds, follow the lines joining the centre with 
the vertices. A dipole can have any one of the six directions 
perpendicular to the faces of the cube. Once the direction of 
the central dipole is fixed, however, each of the neighbouring 
dipoles can have one of only three possible directions: for the 
two molecules bound by — O — H • • • 0< bond the dipole direc- 
tions are away from our cube, and for the other two towards 
it, as can be seen from the figure. The other three directions 
are excluded because they would require — 0 — H • • • H — 0 — 
bonds. It will now be further assumed that correlation with 
more distant neighbours need not be considered. Then the throe 



Fig. 39. A HgO molecule 
with its four neighbours, 
schematically. Full lines 
represent ordinary chemical 
bond, dotted lines indicate 
hydrogen bond. 
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possible dipole directions of the neighbours will be equally 
probable. Hence the average value cosy of cosy, where y is 
the angle between neighbouring dipoles, is 

(17.7) 

because in two positions the dipoles are perpendicular to each 



Fig. 40. An HgO rnoleonle with its four neighbours; dotted lines represent 
hydrogen bond. The arrows at the comers show the possible dipole directions 
of the neighbours for a fixed dipole direction of the central molecule. 


other (cosy = 0) and in one they are parallel (cosy ^1). The 
same result can be obtained if the O — H bond outside the centre 
is permitted to have any direction which does not destroy the 
O — H- • 0 bond; in other words, if free rotation around the 
O — H • -0 bond is assumed (cf. Fig. 41). For this leads to an 
average value ji in the direction of the ()• • H — 0 

bond, if is the dipole moment of a molecule, and hence to a 
value /I/V3 in the direction of the fixed dipole. 

Upon inserting the value of cosy from 17.7 into 8.14, using 
2 ; = 4 for the number of neighbours, we find since 
1+^co^ 1+4/3 = 7/3, 


36, / 7i^+2 \M^/igjyo7 

2€,+7i2\ 3 / 3 AjT 3’ 


(17.8) 


or since the refractive index n == 1*33, i.e. 


3c,/(2e,+n2) - 3/2, 
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Now substituting the following numerical values: Avogadro’s 
number -- fiXlO^^, density 1, molecular weight 18, 
k r- 1-4 X l()“i®, /Xy l-9x 1*77, we find, using 15.2, 

c:::i — — j.e. err:: 55 for T 500 abs., 

(17.10) 

which compares favourably with the experimental value of 78. 

Fig. 41. Direction of the dipole relative to an O — H...O bond. 

At 80” (J,, i.e. at about 353” abs., the experimental value of 

is 60 and the theoretical one is 53. This shows that the ex- 
perimentally determined temperature-dependence is also in fair 
agreement with the IjT law given by 17.10. 

In judging the approximations made above it must be remem- 
bered that the rigid structure assumed for a molecule and its 
neighbours does give only an approximate picture of the correla- 
tion between neighbours. X-ray investigations [M4\ show that 
the average number of neighbours is slightly larger than four ; 
this may be the reason for the larger experimental value of 
Also, with increasing temperature the restriction placed on the 
relative orientations of neighbouring dipoles by the hydrogen 
bond may occasionally be overcome. This would decrease the 
dielectric constant and account for a stronger temperature- 
dependence than that given by equation 17.10. 

It thus seems that on the whole the Kirkwood formula 
accounts in a satisfactory way for the static dielectric constant 
of water. 

No satisfactory theoretical treatment of the frequency - 
dependence of the dielectric properties of water is available at 
present. This would require a generalization of the Kirkwood 
formula applicable to time -dependent fields. The experimental 
results by Collie, Hasted, and Ritson [C3] and by Saxton and 
Lane [^2] point to very simple properties. It seems that the 
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Debye formulae 10.15-10.17 are satisfied, and moreover that 
the relaxation time r is connected with the viscosity r] by Debye’s 
relation 11.30. This means that r and have the same tem- 
perature-dependence; hence T must be a linear function of rj/T. 
Figure 42 shows that this holds for water as well as for heavy 
water. 


3 

T 

2 

1 

0 

Fig. 42, Dependence of relaxation time t on rjjT {rj — viscosity) for water 
and heavy water, according to Collie, Hasted, and Ritson [CS], 

Ketones 

As an example of long-chain substances we shall consider solid 
ketones, and we shall begin with a discussion of their structure 
in the lowest energy states which is realized at sufficiently low 
temperatures. A ketone molecule is obtained from a paraffin 
molecule by replacing one (or several) CHg group(s) by one (or 
several) CO group(s). In a solid a ketone molecule forms a plane 
zigzag, as shown in Fig. 28, § 16. The dipole lies in the chain 
plane and is directed perpendicular to the chain axis. In the 
simplest case the crystal structure is similar to that of paraffins 
which has been investigated by Muller [^M6]. The molecular 
chains are arranged in layers whose thickness is approximately 
equal to the chain length. Within such a layer the molecules 
form rectangular cells with side length 

a Ci:: 5*10“® cm., b c:::i 10“®cm., 

and c slightly larger than the chain length. The dipoles of such 
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a layer are all in one plane, the dipolar plane, and they are 
arranged as shown in Fig. 43. These dipolar planes are thus 
polarized in the 6-direction, and it is of great importance in 
studying the behaviour of the whole crystal to know the relative 
directions of the polarization of successive dipolar planes. First 
then consider tlie position of successive layers in paraffins. 
According to Muller, distinction must be made between two 
cases, depending upon whether the number of C-atoms in a 
chain is even or odd . Fig. 44 shows the positions of two successive 



Fig. 43, Dipole directions in the ground state of 
a dipolar plane. 


layers for the two cases. Now let us attach a dipole to each 
molecule at distances and Cg from the two ends so that c^-fCg is 
equal to the chain length. For both even and odd chains there 
now exist two possibilities. Except for the gap between the 
chains the distance between successive dipolar planes is (a) 
Ci-fCg, or (6) it is alternately 2ci and 2c2. Figure 45 shows that 
for odd chains both cases lead to opposite directions of the 
polarization of successive dipolar planes. For even chains, 
however, this is only so in case (6), whereas in case (a) successive 
layers have the same direction of polarization. This latter case 
will thus lead to a strong polarization of the crystal with all its 
dipoles nearly parallel, while in the other cases the polarization 
of successive layers cancel. 

For paraffins both cases, (a) and (6), are identical. Their 
energies, therefore, differ only by the contributions of the 
dipoles. 

The interaction between dipoles can be considered as composed 
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of the interaction between the dipoles of each single layer and 
of the interaction between dipolar layers. Since the distance 
between neighbouring dipolar layers is large compared with the 
distance between neighbouring dipoles within a layer, it is possible 
to approximate a dipolar plane by a continuously polarized plane. 
The interaction of a given dipolar plane with all tlie others is 



Fig. 45. The two possible positions of 
Fig. 44. Position of paraffin chains dipoles in the ground state of odd and 
in successive layers (schematically), even ketones (schematically). 

then equal to its interaction with the surface charge which the 
other planes produce at the surface of the specimen. For a 
sufficiently large specimen this interaction energy should be 
nearly equal to the self-energy of a continuously polarized 
specimen of the same size and with the same total moment. 
According to the appendix (A 2.ili) this energy is positive, but 
its value depends on the shape. It thus follows that the energy 
of the spontaneously polarized structure is higher than that of 
the other structures, so that we cannot expect the spontaneous 
formation of permanently polarized crystals, corresponding to 
ferro-magnetic substances in the magnetic case. On the other 
hand, if such a structure could be obtained in an even ketone, 
it would probably persist for a very long time, because to reach 
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an unpolarized state the chains of some dipolar planes would be 
required to reverse the direction of their chain axis; this should 
hardly Tie feasible in the solid state. It is possible that the 
required structure can be obtained by solidifying an even ketone 
in the presence of a strong electric field. The strength of the field 
required depends on the shape of the specimen; the most favour- 
able case would be that of a needle-shaped specimen with its 
axis parallel to the field, because such a shape leads to the 
smallest self-energy. 

Now let us consider the dielectric properties of ketones which 
are not pt^rmanently polarized. At low temperature the only 
effect of the field will be to displace sliglitly the ecpiilibrium 
direction of the dipoles , aj)art from the displacement of electrons. 
By comparing, at these temperatures, the static dielectric con- 
stant of a ketone with that of a paraffin of equal chain length, 
Muller \M7]^ found that the contribution of dipole displacement 
is of the order A, ~ 0-1. (17.11) 

This value can be used to estimate the frequency io^j^TT of the 
rotational oscillations of the dipoles. Suppose that the field E 
makes an angle 6 with the equilibrium direction of a dipole^. 
It then exerts a couple iiEmid on the dipole. The restoring 
couple is if I is the moment of inertia and <f> the average 
angular displacement of the dipole direction towards the field 
direction; hence 


4> = 


fjiE sin 0 


(17.12) 


For weak fields, ^ 1, so that fjL(/> sin 6 is the projection of the 

induced dipole in the field direction. Averaging over all direc- 
tion between p. and E, i.e. replacing sin^tf by leads to an in- 
duced moment per unit volume 

2fjL^ENQ 


-34 = 3 


/tog 


(17.13) 


if Nq is the number of dipoles per unit volume. Hence using 1.9 
the contribution to the dielectric constant is 

Sir 


A6=: 


3 /CU§* 


(17.14) 


For the ketone used in Muller’s experiments 2 x lO^^ per 

4»80.n 
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c.c,, (JL 2*5x 10“^® e.s.u. To obtain I we first assume that the 
C-atom remains in its initial position and the 0-atom can be 
displaced. Then / — Ar^, where A 16x 2 x 10“^^ gramme is 
the atomic mass of oxygen, and r cri 10“® cm. Hence using 
17.11 and 17 . 14 , the resonance wave-length Aq is of the order 



( 17 . 15 ) 


On the other hand, if we were to assume that the chain oscillates 
rigidly, then I would be larger by a factor of the order of the 
chain length (~ 20), and hence A^ would be larger by about V20. 
Thus we expect a resonance wave-length between 1/10 and 
1/100 cm. This may well lead to a measurable absorption in the 
centimetre region. 

It was shown in § 10 (cf. Fig. 30 ) that for dilute solutions of 
ketones in a paraffin crystal with longer chain length there are 
two equilibrium positions for each molecule with opposite 
dipole direction and about equal energy. In a pure ketone crystal, 
too, it seems likely that an equilibrium position can be obtained 
by turning a chain plane by 180 ° around the chain axis. This new 
position must, however, have a higher energy than the former 
one, partly due to the alteration of the interaction between the 
dipoles and partly due to non-dipolar interaction. The inter- 
action between dipoles is mainly restricted to those of the same 
dipolar plane. Hence, in view of the strong anisotropy of ketone 
crystals, calculation of this interaction is reduced to a two- 
dimensional problem. In this case the interaction with nearest 
neighbours is a good approximation to the total interaction. 
Since the non-dipolar interaction is also restricted to nearest 
neighbours it is not necessary to take into account any long- 
range interaction if we consider non-polarized states of the 
crystal. It follows from these considerations that with increasing 
temperature an increasing number of dipoles will turn in the 
opposite directio^i, leading to an order-disorder transition in the 
way discussed in § 8 and in the first section of the present para- 
graph. The dielectric constant should, therefore, rise with tem- 
perature, as indicated in Fig. 10, until the transition temperature 
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has been reached. Actually Muller’s experiments show that the 
crystal melts before the transition temperature has been reached . f 
This possibility was envisaged in our previous discussion. 

A transition of the type considered here wg-s found, however, 
by Ubbelohde [ [/ / ] in a paraffin by measurement of the specific 



Fig. 46. Temperature -dependence of the specific heat Cj, of a paraffin (m — 16) 
according to Ubbelohde [U7], showing a transition at about and 

melting. Extrapolations of ‘normal’ specific heat is indicated. 

heat Cp at constant pressure. The turning of molecules into 
different equilibrium positions makes a contribution Ac^ to 
which is characteristic of order-disorder transitions. It leads to 
a sharp increase of near the transition temperature followed 
by an abrupt fall beyond it. Figure 46 shows that Ac^ can be 
separated from the normal specific heat with fair accuracy. 
The transition is connected with a change A/S in entropy which 
is given by the well-knovm thermodynamic relation 

A-S == J ^dT (17.16) 

integrated over all temperatures. If A/S is the entropy increase 
per molecule, then in the disordered state a molecule has 
positions for each one in the ordered state. We should thus 

^S = k\og2. (17.17) 

t This trcuisitioa has now been found by Dr. Vera Daniel (Nature, 1949) in a 
paraffin-ketone solution. Added in proof. 
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Actually, however, the experimental value of AS is^k cmi k log 20. 
Tliis means that above the transition temperature a chain must 
have many more than two positions. It is likely that this is due 
to twisting of cliains [P6]. This may be expected in view of 
the flexibility of chains which has been demonstrated previously 
in the discussion of the relaxation times of dilute solutions of 
ketones (§ 16). 

Additional evidence of the twisting of chains near the transi- 
tion temperature has been given by Muller [ATS]. He measured 



Fig. 47. Temperature-dependence of the dielectric constant of the two di- 
ketones CxoHi80a(l) and CiiHaoOa(2) according to Muller [M8]. The mole- 
cular chains in the imtwisted state are shown as zigzags, the arrows indicating 
the directions of the dipolar ketone groups. 

the dielectric constant of ketones containing two ketone groups, 
i.e. two equal dipoles, for two different cases: (i) the number of 
links between the ketone groups is odd, and (ii) it is even. In 
case (i) the dipoles of the ketone groups are anti-parallel in the 
untwisted chain, so that the total dipole moment vanishes where- 
as in case (ii) they are parallel. Nevertheless, in this case the 
dielectric constant showed a similar increase below the melting- 
point to that of case (i) (of. Fig. 47). This is possible only if chain 
twisting occurs in this temperature region, for otherwise the 
ketone in case (i) would behave as a non-polar molecule and 
have a dielectric constant independent of temperature. 
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18. Ionic crystals 

Ionic crystals are materials in which each lattice point is 
occupied by an ion. Examples are the alkali halides such as 
rock-salt, NaCl, which forms a simple cubic lattice consisting of 
two interlocking face-centred cubic lattices of Na+ and of Cl“ 
ions. Normally the negative ions are much larger than the posi- 
tive ones and often the negative charge of the former overlaps 
into the region of the neighbouring positive ions so that it is 
not always possible to refer to the charge of an ion as an integral 
multiple of In the case of alkali halides the success of Born’s 
[B4] lattice theory suggests, however, that the charges of posi- 
tive and negative ions are well separated. For in this theory the 
binding energy and other properties of the crystals are calculated 
on the assumption that the interaction between ions is composed 
of (i) the attraction between the ions of charge ±e, (ii) the repul- 
sion between nearest neighbours, and (hi) some corrective terms 
due to van der Waals forces. 

The polarization of ionic crystals is entirely due to elastic 
displacements, in both the optical and the infra-red regions. 
Optical polarization, as indicated in § 14, is due to displacement 
of electrons relative to the nuclei and the corresponding reso- 
nance frequencies are in the visible or ultraviolet region. Infra- 
red polarization is connected with the displacement of nuclei, 
and usually is accompanied by a displacement of electrons 
relative to the nuclei. In contrast to its relatively minor part in 
other substances, infra-red polarization plays an essential role 
in the case of ionic crystals. We shall, therefore, have to consider 
it more closely than has been done hitherto. The most general 
(not necessarily homogeneous) polarization can bo considered as 
a superposition of plane polarization waves whose wave-lengths 
range from the order of the lattice distance to infinity, the latter 
corresponding to homogeneous polarization. The frequencies 
of polarization waves fail into two main regions corresponding 
to the frequencies of ultra-violet and of infra-red electromagnetic 
radiation, but the wave-lengths in both regions may range from 
values of the order of the lattice distance to infinity. Polariza- 
tion waves whose wave-lengths are smaller than the size of the 
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specimen can be divided into longitudinal and transverse 
waves, assuming the crystal to be electrically isotropic. For a 
plane wave moving in , the k direction with wave-length ‘iTrjk 
and frequency 

P oc E oc D -= E+47rP x (Ig 

where D and E are the electric displacement and the electric 
field-strength due to the polarization wave. If a magnetic field 
is also assumed to be connected with the polarization wave, then 
with the use of 18.1 the Maxwell equations lead to the well- 
known formula k — nwjc (n — refractive index). Thus if we 
exclude polarization waves with this wave number, i.e. if we 
exclude the presence of electromagnetic radiation, then the 
magnetic field and the magnetic displacement must vanish. 
Assuming also the absence of free electric charges and of con- 
duction currents, we have p == 0, j — 0, JS — 0, and hence from 
the Maxwell equations conclude that (cf. A 1.1 and A 1.3) for 
polarization waves 

div D = 0, curl E = 0, (18.2) 

or using D = E-|-47rP and 18.1, 

div E = —477 div P = — 47rikP. (18.3) 

It follows that the three vectors E» P, and D are parallel. Now 
for longitudinal waves k is parallel to P or D, so that 
div D — ikD 0 

requires that D ~ 0. For transverse waves, liowever, k is 
perpendicular to D or P, so that 18.3 yields div E = 0, which 
together with the second condition 18.2 leads to E = 0. We 
have thus shown that 

D = 0, i.e. E = — 477 P for longitudinal waves; (18.4) 

E = 0, i.e. D = 477 P for transverse waves. (18.6) 

Consider now a spherical specimen whose radius is large com- 
pared with the lattice distance, but small compared with 
Clve\ where v is the frequency of the applied electric field. Prac- 
tically we are thus interested in homogeneous — or nearly homo- 
geneous — ^polarization of the sphere. In this case there is no 
difference between longitudinal and transverse waves. Moreover, 
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it will be shown presently that the frequencies of these long 
polarization waves in a sphere are different from those of a 
specimen which is large compared with the wave-length. For 
this case of homogeneous polarization, if and are, at any 
instant, the optical and infra-red polarizations respectively 
(cf. § 14), the total polarization is given by 

P = 1^,+!^.,. (18.6) 

As shown in § 14, the optical polarization is entirely due to 
the displacement of electrons relative to the atomic nuclei 
(electronic polarization). The infra-red polarization, however, is 
due in part to the displacement of whole ions (considered as 
rigid; atomic polarization) and in part to a displacement of 
electrons which results directly from this change in position of 
the ions. This electronic displacement 'induced’ by the motion 
of the ions vanishes only in exceptional cases. There are thus 
two displacements of the electrons superimposed upon each 
other. However, since the second one is considered produced by 
the change in position of the ions it is essentially an infra-red 
frequency phenomenon and should be lumped with the 'atomic’ 
polarization to form The two components of the polariza- 
tion Pq and P^y can then be considered as independent of each 
other, i.e. they superpose linearly without perturbing each 
other, f If the total polarization were divided into components 
in a different manner — for instance 'atomic’ and 'electronic’ 
components — these could not be considered independently since 
the atomic displacement is usually accompanied by a certain 
amount of electronic displacement. 

Assume now that the sphere is brought into a homogeneous 
field Eq due to sources outside the sphere. Then in equilibrium, 
according to A 2.16, the field inside the sphere is 

E = static case, (18.7) 

if Eq is static. Hence, using 1.9, 

P = -g - static case. (18.8) 

4ireg+2 

t Vibrations are called normal vibrations if they superpose linearly. 
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If the external field is not static, but has a frequency which 
is large compared with the infra-red resonance frequency but 
smaller than the optical resonance frequency, then only optical 
polarization will be excited, and the dielectric constant is e = n^. 
Thus corresponding to 18.8 we find 


Po = 


3 n^-l^ 


(18.9) 


This equation holds for all frequencies smaller than those 
within the optical resonance band, including frequencies smaller 
than the infra-red resonance frequency, although in the latter 
case the relation € = is no longer valid. 

Now let a, olq and be the total, the optical, and the infra- 
red polarizability of the sphere, respectively. Then by definition 
(§ 14), if V is the volume of the sphere, 


P = ^Eo, Po = ^Eo, Pi, = ^^E„, static case, 

(18.10) 

holds for the equilibrium polarizations in a static field. Hence, 
using 18.9, 18,8, and 18.6, if is the radius of the sphere, i.e. 
V = 47ra^/3, 

'^ir ^ ^ ^ n s I n 

n*+2“(e,+ 2)(nH2)- ^ ^ 

We shall now consider vibrations of homogeneous polariza- 
tion (i.e. polarization waves of infinite wave-length) of a sphere. 
We shall assume the material to have one infra-red frequency 
only for this type of polarization and denote by its value 

in the absence of an external field. Assuming the vibrations to be 
harmonic, JFJy must thus satisfy the equation = 0 in 

the absence of an external field. Apart from a constant factor, 
the first term in this equation represents the rate of change of 
momentum of the ions, and the second term is the restoring 
force except for its sign. In the presence of an external field J5 q, 
which may depend on time, the right-hand side of the above equa- 
tion must be equal to the force exerted by this field (apart from 
the constant factor), i.e. it must be proportional to where the 
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proportionality factor is independent of the frequency. This 
factor can, therefore, be obtained by considering the static case 
where = 0 and where the last equation of 18.10 must be ful- 
filled. Thus 


— P -^P 

2 ^ir 




(18.12) 


We shall now return to the study of polarization waves in a 
specimen which is large compared with the wave-length. To 
calculate the infra-red frequency in this case we shall consider, 
within this specimen, a spherical region whose radius is small 
compared with the wave-length but large compared with the 
lattice distance. The polarization within such a sphere is thus 
nearly homogeneous, so that the equations derived above can 
be applied. In particular, equation 18.12 will hold, if Eq now 
represents the field inside the sphere produced by the polariza- 
tion of the surroundings. This means that Eq must be periodic 
with the same period as the polarization wave, and it may be 
expected to be proportional to below, equation 18.15). 

Now within the sphere such a field Eq gives rise to an induced 
optical polarization (proportional to Eq) which can be calcu- 
lated from 18.9 as well as the infra-red polarization given by 
18.12. Thus, assuming that 18.12 gives the same value for 
as was previously considered to exist in the small isolated sphere 
with no applied field, there will be superimposed upon this an 
optical polarization due to the field Eq produced by the remainder 
of the specimen when the small sphere is merely a region in that 
specimen. Hence in a large (compared with the wave-length) 
specimen an infra-red polarization wave is connected with a 
polarization in contrast to the case of a small (compared 

with the wave-length) sphere where a polarization wave is con- 
nected with a polarization only (this was the definition of 
Pfj). This means that the polarization I^>+Po to an infra- 
red polarization wave in a large specimen is composed of atomic 
and electronic polarization in a different way than in a small 
sphere. t We shall show presently that this composition is 

t Pir, which was a normal vibration of the sphere, is no longer a normal 
vibration of a wave in a largo specimen. 
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different for longitudinal and for transverse waves, and that the 
frequencies of these two types of waves are different as well. 

The external field Eo,may be expected to be proportional to 
the total polarization P, say 

477 , 




477 ^ 




(18.13) 


where the constant q will be determined later. Inserting 18.13 
into 18.9, and solving for P^, yields 
„ q(n^— 1 ) 


" n^-j-2—q(n^—l} 

Introducing this again into 18.13 leads to 




(18.14) 


477 


Eo = #*7 1- 


q(n^-l) 


n^-i-2—q(n^—l) 


h 


^Pi 


n ^+2 


irZ 


(18.15) 

Hence with the use of 18.11 and 18.15, equation 18.12 becomes 


— P 4-P 
2 


3g 


This can be wiitten as 


€^+2 n^-\-2—q(n^—l) 






0, 


leading to vibrations with a frequency ajj 27r given by 


= 1- 




3q 


(18.16) 


(18.17) 


(18,18) 


€^.+2 n^+2—q(n'^-~\y 

To determine q we notice that the macroscopic electric field 
E inside the sphere can be considered as composed of the contri- 
bution Eq of the region outside the sphere, and of the self-field 
Eg of the sphere. According to A 2.21 the latter is equal to 
- 477 P/ 3 , so that 

E = En-^P, (18.19) 


E.-fp. 


or inserting E from 18.4 and 18.5, 

477 M 


En 


E+yP 


3 


for transverse waves 


-877 


( 18 . 20 ) 


P for longitudinal waves. 
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Comparison with 18.13 shows that 

1 for transverse waves , ^ , 

(18.21) 

— 2 for longitudinal waves. 


Denoting the angular frequencies of longitudinal and transverse 
waves by and respectively, we find by inserting 18.21 into 


__ ojf _ 71^+2 

col ^s+2 ’ col ^s+2 * 


(18.22) 


The ratio of longitudinal to transverse frequency is thus given 
by (of. references Fll, L4, L5, K2) 



(18.23) 


We thus see that the frequency of longitudinal polarization 
waves is greater than that of transverse waves with an inter- 
mediate value for the frequency of a sphere whose radius is 
small compared with the wave-length. We also find, inserting 
18.21 into 18.14, that the induced optical polarization has 
opposite sign for the two types of waves. Finally it should be 
emphasized again that these conclusions hold only for wave- 
lengths which are large compared with the lattice distance. 

We shall now proceed to calculate the static dielectric constant 
in terms of the infra-red frequencies. We shall again consider a 
sphere liomogeneously polarized by a constant external field 
^0 leading to a static polarization P ^ The change 

F—Fq in free energy due to the infra-red type of polarization 
can be considered as composed of a self-energy proportional 
to P|. due to the elastic displacement, and the interaction 
energy — P^,.Eo V. Thus following a procedure similar to that 
in Appendix A 2.iii we obtain 


F-Fo = -P,,Eo V+^CWPl, (18.24) 


where the constant y in that appendix has now been denoted as 
In equilibrium F must be a minimum if P^^ is treated as 
a vector parameter, because K and Pq are independent of each 
other in the case of a sphere. Hence, as in the development in 
the appendix, ^ 

p ^0 


( 18 . 25 ) 
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Inserting this expression into the last equation of 18.10 and 
making use of 18.11 we find 


6^— __ w.2+2 1 


(18.26) 


Now for elastic displacement a generalized displacement co- 
ordinate O referring to a single unit cell can always be found so 
that the self-energy of the sphere is given by 

(18.27) 


where is the reduced mass of the ions and is the number of 

unit cells per unit volume. This expression must be equal to 
JC'FP?,, so that 


PI 


(18.28) 


Furthermore, the polarization must be proportional to the 
displacement 0 and to the number of cells per unit volume. We 
thus can introduce an effective ionic charge e* byf 

P,, = e*N,Q. (18.29) 

Inserting the value of from 18.28 into 18.26 and making use 
of 18.29 we thus findj (cf. Szigeti, S13) 


€g—n^ __ 477 


(18.30) 


If instead of co^ we introduce the transverse angular frequency 
oj/ from 18.22 the above equation becomes 


e,-n^ = 


/ to 2 + 2\2 e*mo 

[ 3 ) M^cf 


(18.31) 


t It should be remembered here that I^y was defined as polarization due to 
the infra-red normal vibrations of a small sphere. The effective charge e* is 
therefore a quantity which specifically refers to a sphere. 

if An alternative way of deriving equation 18.30 would be to make use of 
the general theory of § 7 according to which (7.21, 7.44) 

c,-I^47rM^ 

€s+2 3F BkT* 

where is the average square of the spontaneous fluctuation of the moment 

of a dielectric sphere of volume V in vacuum. is in our case of harmonic 
vibrations the sum of oi^tical and of infra-rod terms and is proportional to kT^ 
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In deriving this formula we proceeded on macroscopic lines, 
except when introducing by means of equations 

18.27 and 18.29. Of these two constants the reduced mass can 
always be obtained from the mass of ions and from a knowledge 
of the structure with the help of Bom’s \^B4\ lattice theory. 
Calculation of the effective charge 6*, however, would require a 
detailed knowledge of short-range interactions and of the charge 
distribution in the lattice. The reason for this different behaviour 
of mass and charge is due to the fact that the mass is concentrated 
in the nuclei and contributions by the electrons are negligible, 
whereas the charge cannot be considered as concentrated at 
points except for long-distance interaction. 

As an example consider a crystal of the NaCl type. Let 
and M~ be the masses of the positive and negative ions and let 
r+ and r“ be their displacements. In an oscillation of the sphere 
corresponding to a polarization wave with infinite wave-length 
the polarization is homogeneous at any instant. This means 
that the displacements r+ and r~ have the same value in each 
unit cell, and that r+ is opposite in direction to r~. The restoring 
force must then be proportional to |r+— r“| and act in opposite 
direction on the positive and the negative ions. This suggests 
the relation 


o r+— r - 


(18.32) 


because then if we put the restoring force as equal to zb-^red^l 0 
we obtain 18.27 as the self-energy of the sphere. This assump- 
tion for the restoring force requires 

0 == O = 0 . 

(18.33) 


Dividing by and M~ respectively and subtracting the two 
equations we obtain, using 18.32, 

which leads to oscillations of the required angular frequency 
£0, if 


1 
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We have thus calculated the reduced mass and now require the 
effective charge which IiUkS to be obtained by calculating the 
polarization of the sphere for a given displacement Q according 
to 18.29. On the assumption that no non-dipolar interaction 
exists, and that the charges of neighbouring ions do not overlap, 
it is easy to show that c* is equal to e, where ±6 is the charge of 
an ion. The assumption of the absence of non~dipoIar forces is 
entirely fictitious, however. In fact the restoring force is mainly 
due to short-range repulsion between neighbours. A satisfactory 
calculation of c* has not yet been carried out. 

We have thus shown that at the present stage of development 
of the theory all quantities in equation 18.31 except c* can be, 
obtained from experiment. Use of this equation, therefore, 
allows a semi-empirical determination of c*. 'The following table 
(Szigeti, S13) shows that for alkali halides e*/e is smaller than 
unity. This need not be taken as an indication of pronounced 
overlapping of charges, however. It may also indicate that on 
setting up a homogeneous polarization by a displacement of 
the nuclei of ions, an electronic polarization in the opposite 
direction is induced by the short-range forces. 


Table 

(A; = 27rcla)i) 




n* 

A{X 10^ cm. 

g»/c 

LiF 

9-3 

1*92 

32*6 

, 0-83 

NaF 

6-0 

1-74 

40-6 

0-94 

NaCl 

5-6 

2-25 

6M 

0-76 

NaBr 

60 

2.62 

74*7 

0-85 

Nal 

6-6 

2-91 

85-5 

0-71 

KCl 

4-7 j 

213 

70*7 

0-80 

KBr 

4-8 

2-33 

88'3 

0-76 

KI 

4*9 1 

2-69 

10-2 

0-69 

RbCl 

50 

2-19 

84-8 

0-86 

RbBr 

5-0 

2-33 

114 

0-88 

Rbl 

50 

2-63 

129-6 

0-78 

Csa 

7-2 

2-60 

102 

0-88 

CsBr 

6-5 

i 2-78 

134 

0-81 

TlCl 

32 

510 

! 117 

Ml 

CuCl 

10 

3-57 

63 

1 MO 

CuBr 

8 

4*08 

67 

1 10 

MgO 

10 

2*95 

17'3 

1 2x0-88 

CaO 

12 

3*28 

27*4 

2x0-76 

SrO 

13 

3*31 

47 

1 2x0-60 
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For the oxides a factor 2 is separated from the e^je value to 
indicate that an ideal oxygen ion is doubly charged. We see that 
'|e*/e in this case is somewhat smaller but of the same order as 
e*/e for alkali halides. It is also noticeable that the high dielectric 
constant of TlCl does not lead to an excessive value of e*/e, but is 
largely due to the high refractive index. In fact, even in the case 
of Ti 02 , whose dielectric constant is larger than 100, a value 
0-7 is obtained (referring to the oxygen ions), as has 
been shown by Szigeti The high refractive index, together 

with the high charges of the ions, are mainly responsible for the 
high dielectric constant. 


In the above discussion we have made use of equation 18.31 
and not of 18.30 because and not can be obtained from 
experiment. Equation 18.31 is very similar to an equation 
derived by Born [B4] by an approximate method. This equation 


does not contain the factor 


1 ^ 2 ^ 2 


and it replaces e* by the 


actual ionic charge. 

It should bo remembered that, in view of the relation 18.22 
between and ca^, equations 18.31 and 18.30 are equivalent. 
This raises the question of the possibility of permanent polariza- 
tion of ionic crystals. For on solving 18.30 for it is found that 


Cg 00 as 


47r 71^+2 e*‘Wo 
3 3 


(18.36) 


On the other hand, equation 18.22, for ->oo, requires ca^ -> 0 
(if we exclude ca^-^^oo), which Avith 18.31 would also lead to 
€,'-> 00 . Investigations on these lines should be of importance 
in view of the properties of crystals like barium titanate. They 
have not been developed far enough, however, to be included in 
the present book. 
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(i) Conservation of energy 

Macroscopic oloctromagnetic theory is based on the Maxwell equa- 
tions, 

curl E — — - ^ (A ].l), div D - ^np (A 1.3), 

c ot 

ciirl H = - ^+— j (A 1.2), iliv B -= 0 (A 1.4), 

c ot c 


where H and B are the magnetic field strength and induction respec- 
tively, p is the density of the true charge, j the density of the conduction 
current; for E and D cf. § 1. To allow a unique calculation of the field 
vectors E, D, H, B from p and j, the Maxwell equations must be supple- 
mented by two further relations of the typo 

E-E(D) and H H(B). (A 1.5) 

These relations are not included in the fundamental equations for the 
electromagnetic field, but are characteristic of the type of material 
employed. For substances of interest in this book 

H = B (A 1.6) 

can be assumed to hold to a good approximation. 

To introduce the energy law, multiply equation A 1.2 by E, and sub- 
tract from it equation A 1.1 multiplied by H. Use 

H curl E-E curl H = div[E x H]. 

Then with equation A 1.6, after multiplication by c/47r, 

Books on electromagnetic theory discuss in detail the fact that the 
last term of this equation represents, per unit volume, the rate of con- 
version of electromagnetic energy into other typos of energy (heat, 
kinetic energy of particles, etc.) in so far as this is connected with con- 
duction currents. c[E x H]/47r is the Poynting vector representing the 
rate of flow of energy; the second term, therefore, gives per unit volume 
the rate of efflux of electromagnetic energy. Then from conservation of 
energy it follows that the first term must represent the rate of change 
of energy content per imit volume, provided that there is no flow of energy 
of another type (e.g. heat currents). 


Now 


JLh— - ~ — 

4^r dt dt Bit* 


so that H^ISrr is readily recognized as the density of the magnetic energy. 
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It is, then, often suggested that 

S / 

should represent the density of the electric energy. The integral is con- 
sidered to (3xtend from E -- 0 to the actual value. In fact, if E — D/e^, 
one finds, for a field-strengtli 

(A 1.8) 

0 

if can be considered as constant. Such a procedure is, however, possible 
only if E is a univalent function of D, for otherwise the integral has no 
unique meaning. Thus only in this case can the electric energy be defined 
without further investigation. 

In general, the conclusion to bo drawn is that 

:^EdD (A 1.9) 

477 

represents the change of energy density (not necessarily electric energy) 
connected with a variation of D by dD, if ther(3 is no flow of energy of any 
other type. It has been shown in § 3 how the total energy con be calcu- 
lated from this expression. 


(ii) Conduction current and energy loss in periodic fields 

As in § 2, assume a homogeneous field with E -- E^ cos wt, and 




€2 ^E 
oj dt ’ 


Then using d^Ejdt^ = —w^E, and assuming the conduction current y to 
vanish, the right-hand side of equation A 1.2 becomes 


1 ^ 

c dt c dt^ c 


(A 1.10) 


If, on the other hand,y 0, and if Ohm’s law holds, i.e. if (cr is the con- 

ductivity) 


and if furthermore €2 
reads 


0, then the right-hand side of equation A 1.2 
c dt c 


Comparing equations A 1.10 and A 1.11 it follows that in periodic 
fields the introduction of two dielectric constants €i(a>) and € 2 (a>) (cf. § 2, 
equation 2.8) is equivalent to using a single €i(a>), putting Cg == 0, but 
introducing instead a frequency dependent conductivity cr(to). The two 
representations are connected by 

(x((o) = (A 1.12) 

477 


4980.11 


M 
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J oule’s law for the rate L of loss of energy per unit volume then becomes 

\ctEI, (A 1.13) 

where the bar has been used to represent the average over one period. 
Inserting a from A 1.12, this expression is found to bo equivalent to equa- 
tion 3.15. 


(iii) Relation between €i(co) and CsCco) 

To derive equations 2.16 and 2.17 apply the theorems of Fourier 
transformations to equations 2.14 and 2.15. 

Then from 2.14, 


00 

I f djx, (A 1.14) 


00 
2 r 

and from 2.15, ol(x) = — I € 2 (fi)^in fjLX d/x. 

0 

Introducing A 1.15 into 2.14 yields 

00 . 00 . 
€i(co)‘-€^ ~ ^ f J' £a(/x)sin^a7 dfij 

6 ' 0 ' 

<x> . R . 

= - lim ,/ dJe^in) J COS ojx sin fix dxj 
0 ' 0 ' 

_?lim 

\ /x+ta 


(A 1.15) 


-cos(/x— a>)J? 


fl — W 




(A 1.16) 

Now the integrals containing the cos-terms vanish if B-> oo. The re- 
maining terms lead immediately to equation 2.16. 

Finally, introducing A 1.14 into 2.15 gives 


00 . 00 . 
ej(a>) = ^ / da;(sina>a; J {€i(ii)—€a^}co&fixdfi\ (A 1.17) 
0 ' 0 

which yields equation 2.17 by a similar method. 


(iv) Relations between dielectric constants and optical constants 

From the Maxwell equations A 1.1-A 1.4 the wave equation can be 
derived. For this purpose apply the operator curl to A 1.1 and 1/c Bjdt 
to A 1.2. Then, using A 1.6, H can be eliminated. Assuming periodic 
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solutions, we have D == €(a^)E, and hence in the absence of free charges 
(i.e. p ~ 0) and of conduction currents (j — 0) we find 




(A 1.18) 


if use is made of the relation (since from A 1.3 and p ~ 0, div E — 0) 
curl curl E = graddivE— V^E ” — V^E, 

Now assume a solution representing a wave penetrating the dielectric 
in the ir-direction, ^ 1.19) 

where A is a constant vector perpendicular to the a;-direction. Here by 
the usual definition of the optical constants n is the refractive index and 
K is the absorption coefficient. They can be expressed in terms of the 
complex dielectric constant e (cf. 2.8), for after inserting A 1.19 into 
A 1.18 we find + + 


Hence 

Cl = n^—K^ 

(A 1.20) 


€2 = 2n/f. 

(A 1.21) 

Alternatively, using A 1.12, 

2tt 

riK ~ — or. 

OJ 

(A 1.22) 


A 2. DIPOLE MOMENTS AND OTHER ELECTROSTATIC 
PROBLEMS 

(1) The basic problem 

We consider an infinite homogeneous dielectric with static diolectricf 
constant which contains a spherical region of radius a and dielectric 
constant € 2 . We wish to calculate the electric field duo to any of the 
following sources: 

(а) Outside sources leading to a constant field Eq^ (say in the z-direction) 
at a large distance from the sphere. 

(б) A point dipole p. (say in the z-direction) at the centre of the sphere, 
(c) An extended dipole M (say in the z-direction) inside the sphere 

arising from a homogeneously polarized sphere of radius a, 

M==ya»P„ (A 2.1) 

Pq being the polarization. 

Let <D be the electrostatic potential so that 

E == —grad <E> (A 2.2) 

is the field -strength. O satisfies the Laplace equation 

V* O 0 (A 2.3) 

f In contrast to the denotation in the main sections, and now represent 
static dielectric constants. 

4980.11 3^2 
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subject to the following conditions: if r is the distance from the centre of 


the sphere and $ the angle between r and the s-axis, then in 
case (a) E — E^^, i.e. <E> = —E^r cob 6 if r ^ a. 

(A 2.4) 

case (b) 

_ u cos^ ^ 

(I> if 0, 

(A 2.5) 

case (c) 

D = €2 if 

(A 2.6) 


D = E if r > a. 

(A 2.7) 


wb(^ro D is the electric displacement. In all three cases the normal com- 
ponent of D and the tangential compound Eq of E must be continuous 
at r = a. 

The calculations are simplified by treating all three cases together 
because they all lead to the same angular dependence of O. In general 
theory O is developed into a series of spherical harmonics. Owing to our 
three conditions, and to the boundary conditions, only terms proportional 
to cos^ appear. For this angular dependence the general solution of 
A 2.3 is given by 


containing two arbitrary constants A, B. They will have different v'alues 
outside (Ai, B-^) and inside {A 2 , B^) the sphere and have to be determined 
from the boundary conditions. Using A 2.4 it follows that -■ and 
hence 

(A 2.8) 


<1, = 


Also using A 2.5, we have —Ai ~ i*<^‘ 

^ = {^-B2r)coB0, 
According to the boundary conditions 


Eq ™ 


1 ^ 
r d$ 


r < a. 


(A 2.9) 


must be continuous oX r ~ a, so that 

- -. 77 .+®- 

Similarly, using A 2.6 and A 2.7, continuity of Dcosd at r ~ a 
yields . ^ 

(A 2.11) 

From these two equations one finds 


2 fj. iirPc 

O ® 2^1 €2 ^ 2 € i - j -€2 2 Ci “|-€2 

R _ p .2 Ci-c, fi ^ttPc 

“ 2€i+€2 Zci+e* 2ei+«,' 


(A 2.12) 
(A 2.13) 
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Inserting A 2.12 and A 2.13 into A 2.8 and A 2.9 gives the solution of the 
general problem. 

For the discussion we shall now separate the three cases. 

Case (a): No internal sources, /i, = 0, = 0. In this case the field 

inside the sphere is equal to and will bo denoted as cavity field 0\ 
Thus from A 2.13 with A 2.2 


is the field inside the sphere. If in particular € 2 = 1 , i.e. for an empty 
sphere o 

G'(e, ^ 1 ) = G = 2 ^ (A 2.15) 

If, on the other hand, the sphere is in vacuum, then = 1, and the field 
inside becomes 

(A 2.16) 


The field outside the sphere according to A 2.8 is composed of the field 
Ego at infinity and of a dipolar field with potential (using A 2.12) 


^1-^2 
’261 + 62 




(A 2.17) 


Case (6): Point dipole, no external field, Eg^ ~ 0, “ 0. In this case, 

according to A 2.9 and A 2.2, Pg represents the deviation of the field 
inside the sphere from a purely dipolar field, i.e. is the reaction-field 
P' acting on fx. Hence with A 2.13 


R'=: 


2 €i--C2 , 


(A 2.18) 


if in particular 


€2 2€i-f€2 

1, then the reaction field will be denoted by R, and 


Outside the sphere we have a dipolar field with potential 
^ AiCOS^__ 3 IXCOS0 


(A 2.19) 


(A 2.20) 


r* 261-l-ea 

Case (c): Homogeneously polarized sphere, no external field, Eg^ = 0, 
fx — 0, Consider first the sphere in vacuum, i.e. = eg = 1. Then Pg 
represents the self-field inside the sphere which using A 2.13 and A 2.1 is 
thus given by 

(A 2.21) 




M 


Now surround this sphere by a medium with dielectric constant cj. Then 
the increafle of the field inside the sphere, Pj— is the reaction-field E, 
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Thus using A 2.13 with C 2 = 1» 

4ttP„ 47rp 2(t|— I) 47r 2(ei— 1) M 

2€i+1'^3 2€i+ 1 3 " 2ci+l a» 

(A 2.22) 


in formal agreement with A 2.19. 

Thus, if the sphere contains a point dipole (x surrounded by a homo- 
geneously polarized sphere with moment M, then its reaction field is given 

2^-1 

R = -3^^(M+|x)-9(M+h), (A 2.23) 

where g is defined in A 2.19 and M+p. is the total dipole moment of the 
sphere. 

For regions outside of the sphere also case (6) turns into case (c) if fi 
is replaced by M. This holds even if # 1. Thus as far as the field in the 
region external to the sphere is concerned, cases (6) and (c) lead to the 
same results. 


(ii) Dipole moments 

The potential of a rigid (non-polarizable) dipole with moment Pg in an 
infinite medium of dielectric constant ej is given by 

0 = (A 2.24) 

€i r® 

Now let us consider a model of a molecule consisting of a rigid dipole p 
at the centre of a sphere with dielectric constant eg. (The same model wiJl 
of course apply to an extended dipole with the same moment correspond- 
ing to case (c).) If this molecule is embedded in a medium of dielectric 
constant Ci, then according to A 2.20 (and to the remarks at the end of the 
discussion of case (c)) the potential outside of the sphere can be repre- 
sented by equation A 2.24 if 

Therefore Pg is defined as the external moment of the molecule in the 
medium Ci. 

It should be realized that the moment of the molecule in vacuum, p^, 
is different from p. For p^ is defined as the external moment in a medium 
for which Ci == 1. Thus the vacuum moment is given by 

(1, = . (A 2.26) 

It follows that the external moment pg in a medium can be expressed 
in terms of p^,. Using A 2.25 and A 2.26, 


€2 4" 2 Sci 
3 2€i-j“€2 




(A 2.27) 


In contrast to the external moment pg — ^which is the moment a point 
dipole must have to produce the same field a& the molecule — ^the internal 
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moment is the actual moment of the molecule if embedded in the 
medium Cj. In vacuum — 1) the external and internal moment are, 
of course, equal. In another medium, the internal moment differs from 
the vacuum moment /x^ by the polarization of the molecule through the 
reaction field. Hence equation 6.18 for is obtained and this, in turn, 
with the help of A 2. 1 9, leadi> to equation 6 . 20 . Equation 6.20 can also be 
derived by using directly the properties of the external moment. As we 
have seen above, the field produced outside the sj)here does not depend 
on whether the moment is produced by a point dipole or by imiform 
polarization of the sphere if this leads to the same moment. Therefore 
the molecule (as represented by a sphere with dielectric constant €3 and 
a dipole p, at its centre) must have the same moment as a sphere with 
dielectric constant Cj which is in a medium with the same dielectric con- 
stant €1 and has a moment /Xg at its centre, because by definition of 
both produce the same field outside the sphere. Thus 

= |Jie+ J P (A 2.28) 

sphere 


showing that |X{ is composed of the moment of the rigid dipole pig and the 
moment contained in a sphere around it. Now if points in the 2 :-direc- 
tion, the integral must be a vector with the same direction for reasons of 
symmetry. Hence using 

-(€,- 1 )-^ (A2.29) 


and 


cos 6 z 

l( 

11 

^l 

11 

dz\ 


we find with the help of A 2.24 that 

/££)-= ISO 

and J V‘>(^)dT = - 477 . 

Inserting A 2.30 into A 2.28 yields 


(A 2.30) 


which is identical with 6.20. 

The above calculation also shows that the dipole moment contained 
in a sphere surrounding a dipole is independent of the size of this sphere. 
Hence the moment contained between two spherical surfaces inside the 
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dielectric vanishes. This holds even if the two spheres are not con- 
centric, f 

Thus, if we consider a large sphere inside an infinite dielectric of dielec- 
tric constant and this sphere contains a molecule with internal moment 
then the moment of the large sphere is also 

This is, however, no longer the case if the large sphere is not embedded 
in its own medium because then the field outside the molecule is different 
from the field of a point dipole. By analogy with A 2.26 one should expect 
the moment in this case to be 

(A 2.32) 

as is actually confirmed by calculation if the radius of the sphere is large 
compared with that of the molecule. 

Wo have thus introduced five different dipole moments, /x, 
jig. To prevent confusion we give a brief summary. 

fjL has a meaning only in terms of the special model used. 

= moment of the molecule in vacuiim. 

fjLi — moment of the molecule in a medium Cj 

= moment of a sphere (containing the molecule) within an infinite 
medium €i. 

— moment of a rigid point dipole producing in a medium €i the same 
dipolar field as the molecule. 

fig = moment of a dielectric sphere in vacuum containing the molecule. 

To summarize the formulae: 


3 2€i+e2 


ei(Cj + 2)'*‘ 


glL+i-v -- J ». 

€j + 2 2e. + l‘*' €2+2** 


(A 2.33) 


f The moment in any direction, say with radius vector s, is proportional to 
the integral J dQ>ld8 dr extended over the space between the two spheres with 
radii and This integral can be transformed into the difference of two 
integrals over the surfaces of the two spheres (which are designated as external 
surface e and internal surface i): 

rj J Og cos 0 dO, J cos ^ dCl. 

e i 

Here dQ is the element of the solid angle for each sphere, <De and are the 
values of O at the two surfaces, ip is the angle between s and the normal to 
the surfaces pointing away from the centre. Both and can be developed 
into spherical harmonics of which only the first representing a dipolar poten- 
tial of a dipole at the centre of the sphere concerned gives a contribution to 
the integral. These terms are proportional to cos 0/rJ and cos respec- 
tively, which means that the two surface integrations just cancel. 
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(Ui) Self-Energy 

The free energy F of a dielectric polarized by a field can be considered 
as composed of three terms: ( 1 ) the free energy Fq of the field and of the 
dielectric before the latter is brouglit into the field; (2) the energy of 
interaction between the polarized dielectric and the field; and (3) the free 
energy required to polarize the dielectric. The latter will be called the 
self-energy Fg. If the dielectric is polarized homogeneously, then it will 
be assumed that Fg — yM^, where M is the moment of the dielectric. 
The energy of interaction is equal to — (MEq) if Eg is the field (supposed 
homogeneous) before the dielectric is inserted. Thus 

F-Fo - ~(MEo)-f (A 2.34) 

In equilibrium F must be a minimum where M ~ (M^, M^) is treated 

as a parameter. Hence from 


we find 


dF __ ^ ^ 

dM^ ” dMy ~ dM^ 

2yM = Eg or y ~ 


- 0 , 

A 

2M‘ 


(A 2.35) 
(A 2.36) 


Hence our assumption Fg -- yM^ implies that M is parallel to Eg. The 
value of y depends on the shape of the dielectric. Thus a sphere of radius 
a brought into a homogeneous field Eg has a moment 
M = a3Eg(€,- !)/(€, + 2). 


Hence 


_ I ^ j_ 1 m 

2 eg- 1 a^’ ’ 2 e,-I a® ■ 


(A. 2.37) 


A slab of volume V with surface perpendicular to Eg, on the other hand, 
has a moment M =“ FEo(e^— l)/47r€j,. In this case 


27r 


F, 


27tM^ 


(A 2.38) 


6,-1 V 

In both cases, of course, Fg is proportional to the volume if the polariza- 
tion ikf/F is constant. 


A 3. THE CLAUSIUS-MOSSOTTI FORMULA 
In § 5 an equation for the static dielectric constant e, was derived 
(5.13) which holds exactly for the model used in that section. This is 
usually called the Clausius-Mossotti formula. Later in § 6 it was found 
that a similar formula holds approximately for liquids of non -polar 
molecules, but that the approximation required negligible short-range 
interaction. In § 8 it was pointed out that by neglecting short-range 
interaction the general theory of § 7 will lead to either the Clausius- 
Mossotti or to the Onsager formula, depending on whether one considers 
non -polar or polar spherical molecules. The proof of this will be given at 
the end of this section. 

Discussions as to whether or not the Clausius-Mossotti formula holds 
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exactly have formed part of the literature on dielectrics over a period of 
many years. Controversial conclusions have been drawn, mainly due to 
a misunderstanding of the significance of this formula. In fact one should 
distinguish between a macroscopic and a molecular formula. The same 
mathematical symbols are normally used in the two cases and, as a 
result, they are often confused with each other. The macroscopic formula 
is valid exactly, but the molecular formula holds only subject to the 
conditions mentioned in § 15. 

To derive the macroscopic formula, consider a sphere of a continuous 
isotropic dielectric to be brought into a constant electric field f. Accord- 
ing to equation A 2. 1 6, the field inside the sphere is tPien given by 



if €g is the dielectric constant. Let Mjg be the dipole moment induced in 
the sphere. Then the polarizability ot^ of the sphere will be defined by 
(the suffix m indicates ‘macroscopic’) 

(A 3.2) 

On the other hand, in view of 1.9, 

= (A 3.3) 

47r u 


if V is the volume and is the radius of the sphere. Equating A 3.3 and 
A 3.2 and introducing / from A 3.1 loads to the Clausius-Mossotti formula 


€s + 2 


(A 3.4) 


An alternative way of deriving this formula is to consider a spherical 
region within a homogeneous dielectric. If a constant field E is produced 
in the dielectric, the part of the field inside the spherical region which is 
due to sources outside this region (the inner field) is given by expression 
5.9 and is therefore identical with our field f (A 3.1). The remaining 
development is then identical with the one given above. 

Equation A 3.4 is always correct if the spherical region is sufficiently 
large, so that the material contained in it may be considered from a 
macroscopic point of view. From a molecular point of view equation 
A 3.4 has, however, no significance. To give it such a significance the 
polarizability a must be expressed in terms of other quantities whose 
values will not be found experimentally by measurement of the static 
dielectric constant but by experiments of a different nature. 

Equation A 3.4 finds its main application in the case of dielectrics in 
which polarization is connected with elastic displacement of charges 
(case (i) of § 4). The most general displacement of the dielectric material 
in a sphere may then be developed in terms of its normal vibrations 
(cf. van Vleck, V 3) ^ and it is then found that the polarizability a is a con- 
stant independent of temperature if the density of the substance is kept 
constant. This procedure is of importance because it shows in a very 
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general way that the assumption of elastic displacement leads to a 
temperature-independent dielectric constant. 

We shall now discuss the ‘molecular’ Clausius-Mossotti formula. In 
this case it is assumed that the macroscopic polarizability can be 
expressed in terms of a molecular polarizability oc — ^ which is a constant 
independent of temperature or density — by 

— N(x (A 3.5) 

if iV” is the number of molecules in the macroscopic sphere of radius 
Furthermore, if 


(A 3.6) 


IZT s — YL — ^ 

3 “ ~'N~ S'N 

is the volume occupied per molecule, insertion of equation A 3.5 and 
A 3.6 into A 3.4 yields 

g: — 1 rv 

(A 3.7) 


6^ + 2 

This molecular formula has the same mathematical structxiro as formula 
A 3.4, but its significance is different. In equation A 3.7 the polariza- 
bility oc is a property of a single molecule independent of macroscopic 
parameters. The quantity Zj^Tra^ is equal to the number of molecules 
per unit volume which is determined by the molecular weight W and the 
density d, o ^ 


47Ta» ■“ 


A, 


(A 3.8) 


where A is Avogadro ’s number. Hence aja^ is proportional to the density 
of the dielectric, which may bo varied by altering the external pressure. 
No such conclusion can bo drawn from equation A 3.4, where the macro- 
scopic polarizability cx„,, may depend on the density in an unspecified way. 

The crucial step leading from equation A 3.4 to A 3.7 is contained in 
the hypothesis A 3.5. One should expect A 3.5 to be correct only in the 
absence of short-range interaction between molecules, because such an 
interaction would influence the reaction of a given molecule to the 
field. Actually in the two cases in which we have derived the molecular 
Clausius-Mossotti formula (5.13 and 6.34) short-range interaction was 
absent. 

Finally let us derive A 3.7 and the Onsager formula by applying the 
general theory of § 7 to a model in which short-range forces are entirely 
absent. Also it will be assumed that the unit cell is spherical and that the 
average polarization outside the imit cell is equivalent to that of an 
isotropic continuum polarized by the moment m of the unit cell. In 
this case the average moment m* of a sphere (embedded in a larger speci- 
men) containing the unit cell is equal to the moment m of the unit cell, 
as was shown in § 7. Hence the spherical region introduced in § 7 may be 
taken as equivalent to the unit cell, and equations 7.11 and 7.33 become 
identical if 7.34 with N = 1 is taken into accotmt. Thus, using A 3.8, 

3c, ^ 

2c, -fl kfa^* 




(A 3.9) 
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where is the average square of the spontaneous dipole moment of the 
spherical cell embedded in its own medium. Assuming first that the unit 
cell contains a rigid dipole with moment p, m — p, and — p,*, the 
Onsager formula follows from A 3.9. Consider, on the other hand, a 
spherical polarizable molecule without a permanent moment. Such a 
molecule can be represented by an elastically bound charge e. Then if r 
is the displacement, the internal energy Ui must be quadratic in r, i.e. 

(A 3.10) 


where c is a molecular constant. To express this constant in terms of the 
polarizability, assume the molecule to be brought into an external field f. 
Then its potential energy is given by 


efr. 


(A 3.11) 


The equilibrium value r of r is obtained by making this expression a 
minimum. Hence 


er : 


(A 3.12) 


is the average moment of the molecule. Since this must be equal to of, 

we obtain ^ 

«=-, (A 3.13) 

c 

and hence with A 3.10 ^ 

= (A 3.14) 

To obtain we require the energy U == Ui-^U^ because, as in 7.12, 
using m = er. 


00 .00 

^ = e* I r’‘e-VI’^Tridr j e-t^Wdr. (A 3.16) 

0 0 

Now the external energy is given by — JmR where R is the reaction 
field, as can be shown by the argument which led to 7.18. Therefore 
making use of 5.10 and the fact that ^Tva^j^ = F, 

V = iizJir* 

‘ 2€,+ 1 a3 2e,+ l 

Introducing U ~ into A 3.15 and integrating yields 

^ 3 kT 2a^(2€,^\)oL 


(A 3.16) 


2a3(2e,-fl)-2(€,-l)a- 


(A 3.17) 


Inserting this expression into A 3.9 leads at once to equation A 3.7. 

In the above derivation the importance of the external energy 
including the interaction with the reaction field should be stressed 
It is interesting to compare this with the case of a non -polarizable 
molecule with dipole moment p (Onsager case), where the interaction 
with the reaction field is irrelevant. In this latter case the only variable 
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is the direction of {i., but its amount, /x, is a constant. Thus since m = |x 
in this case, == — constant, so that Ug according to A 3.16 is a 

constant as well. Hence in A 3.16 the terms containing Ug can be taken 
outside the integrals and, therefore, cancel. 


A 4. SHAPE OF ABSORPTION LINES {F9) 

We shall consider a set of linear harmonic oscillators of mass Af, 
charge e, proper frequency a)o/27r. Their number N per unit volume is 
assumed to bo so small that interaction between them can be neglected. 
Then in an external field E, which may depend on time, 

x~ —a)lx-\-^Ecosd (A 4,1) 

if X is the displacement of the charge e, 6 the angle of E with the x-direc- 
tion, and a dot represents differentiation with respect to time t. The 
oscillators are supposed to make frequent collisions with a medium which 
is in thermal equilibrium so that they also tend to reach equilibrium. To 
describe their behaviour we introduce a distribution function /(ir,:c) so 
that f(XfX) dxdx represents per unit volume the number of oscillators 
with a displacement between x and x-\-dx, and a velocity between x and 
x-\-dx. Hence 

00 

J J f(x,x) dxdx = N. (A 4.2) 


With the help of the distribution function the polarization P in the field 
direction is obtained as 

CO 

P — ecosd J J xf(XfX) dxdx. (A 4.3) 

— 00 


Let us first consider the static case for which the field is independent of 
time, E ~ E^. If U is the energy of an oscillator, i.e. 

JJ = UQ—eExcos6, (A 4.4) 

Uo = iMwlx^+iMx\ (A 4.5) 

then, in equilibrium, according to the Boltzmann theorem, 

f{x,x) = (A 4.6) 


where C is independent of x and x and can be determined with the use 
of equation A 4.2. Now using A 4.4 and E — Eq, 


^-U{x.x)lkT ^ Q-UalkT 


eEx cos 0 


I (eEx cos BY \ 

■*■21 ¥ t ~) 


(A 4.7) 


and hence we can write, considering A 4.6, 

eEcosOBfa 
Meal a* 


(A 4.8) 


where the dots represent terms proportional to higher powers in 
In the calculation of the polarization P we shall be interested in linear 
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terms in E only. Therefore we need not consider these higher terms; to 
be exact we should also prove that the series for P, in terms of Eq con- 
verges, but we shall omit this. In A 4.8 

= Ce-VolkT (A 4.9) 


represents the equilibrium distribution in the absence of a field. The fac- 
tor C is the same as in A 4.6 if in calculating the inte^gral A 4.2 for the 
latter case only terms which are linear in, or independent of, Eq are con- 
sidered. 

Now inserting A 4.8 into A 4.3 we find for the polarization Pq in a con- 
stant field 

00 00 

Pq ~ ecosdl j j xfQdxdx~^^^^j~ J |* x~ dxdxj. (A 4.10) 
— 00 —00 

ISj&tQ the first term vanishes because is an even function in x. Inte- 
grating the second integral by parts we find, using A 4.2, 


e^Ecos^ON 

Mwl 


(A 4.11) 


Averaging over all directions of the field loads, of course, to a replacement 
of cos^^ by 1/3. Hence the contribution Ac to the static dielectric con- 
stant is given by 


A. -5 


eW 


E ZMojX 


(A 4.12) 


Conside^r now the case of a periodic field E described by 

E E^e-^^K (A 4.13) 


Again we shall introduce the distribution function/ satisfying A 4.2, but 
now depending on time t. We shall calculate / by considering its rate of 
change dfjdt for a given value of x and x. This is composed of two terms, 
one due to collisions of the oscillators with the surrounding mediiun and 
the other due to their motion. Denoting by (^//^Ocoii rate of change of 
/ due to collisions the simplest assumption we can make is to put 

(I) (A 4.14) 


in which the relaxation time t is assumed to be independent of x and x. 
This implies that equilibrium (/equ) is approached exponentially. 

For time -dependent fields we shall assume that A 4.14 still holds if 
/equ is the equilibrium distribution corresponding to the field at the time 
in question. Thus making use of A 4.8 and A 4.13, 


/egu — /o' 


cEq cos 9 dfQ 


Mwl * 

For the distribution function / in this case we shall assume 
f(x, X, t) == /o(a?, x)^g(xy x)Eq cos 6 


(A 4.15) 


(A 4.16) 
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where g{Xj x) is a function of x and Xt independent of L To satisfy A 4,2 
this requires 


J J g(x,x)dxdx — 0. 

— CO 

Inserting A 4.16 and A 4.16 into A 4.14 now leads to 


(A 4.17) 


\dt)n 


(A 4.18) 

^coll 

There is a second contribution to dfjdt due to the motion of the oscilla- 
tors. This contribution will be denoted by so that 


dt \dtl^oll 


(A 4.19) 


To find an expression for (df/dt)jf^wo notice that all the oscillators contained, 
at a time t, in an interval AxAx, the number of such oscillators being 
f(Xf X, t) Ax Ax, were at the time < — in an interval of the same sizef Ax Ax 
with coordinates x—xht, x—xht provided the number has not been 
changed by collisions; all/(a;,:r, ^—8^) Ax Ax oscillators contained in this 
interval at the time have meanwhile passed into another interval. 
Thus 

^ {/(^ ~xhtiX—xht,t—ht)— f{x, X, ^ — 80 } 
is the rate of change of / due to motion. Hence with 8^ -> 0, 


m 

\dtL 


(Jx ' dx * 


or making use of A 4. 1 , 




cEq cos 6 df 


Now from A 4.9 and A 4.6 it follows that 


0. 


(A 4.20) 


(A 4.21) 


(A 4.22) 


Thus inserting from A 4. 16 into A 4.21 and neglecting terms proportional 
btain 

®). = (A 4.23) 


Also from A 4.16, 


dt 


= —ioigEQe-^^^^GosB. 


(A 4.24) 


Inserting now A 4.24, A 4.23, and A 4.18 into A 4.19 we find that g 
satisfies the differential equation 


- -Sx'^+'^Odx'" M 8x 


Mcoir dx* 


(A 4.26) 


t This can be proved with the help of A 4.1. 
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To solve it we put 


dx* 


(A 4.26) 


where a and h are constant. 

This does not represent the inost general solution of A 4.25, but it has 
been shown by Huby [H3] that it is the only solution satisfying all 
requirements. 

The polarization P according to A 4.3, A 4.16, and A 4.26 is thus given 
P = ecoa^6E„e~'‘^* j j 

-■» (A 4.27) 

To find the constant a insert A 4.26 into A 4.25. Making use of 


and 

we obtain 

Hence 


/ .d^^d\dfo .,a/o 

\ dx ^ dx/dx dx * 

\ ^T/ Mcoir 

(_ia.+i)6 = 


e ojoT^+l— 


McOq a>o -|- ( 1 — icor)^ 
__ 6 1/ l“-^a>oT 

Mwq 2\1 — i(co-\-(x}a)T 


1 +^’^>0'^ 
i{w-\-a}Q)T^ I — i(oj—o)o)r^ 


)• 


(A 4.28) 
(A 4.29) 
(A 4.30) 
(A 4.31) 

(A 4.32) 


Together with A 4.27 and A 4.12 this leads to equation 13.9. 
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a-Bromo Naphthalene, 123-4. 
Absorption, 70 ff., 90 ff. 

— coefficient, 98 ffi, 117 ff., 173 ff. 
Alkali halides, 158. 

Ammonia, 118-19. 

Anisotropy of polarization, 112. 
Argon, 111. 

Atomic polarization, 106 ff. 

Atoms, 104 ff. 

Benzene, 108. 

Benzophenone solution, 120—1. 
Benzyl alcohol resin, 135. 

Caesium halide, 158. 

Calcium oxide, 158. 

Camphor, 123, 124. 

Carbon dioxide, 106, 107, 108, 115. 

— monoxide, 105. 

— tetrachloride, 108, 116, 117. 
Cavity field, 25, 34, 39. 

Cetyl palmitato, solution in paraffin 
wax, 1 38. 

Chlorinated diphenyl, 135. 

Chlorine, 105. 

Chloroform, 116, 123, 124. 

Clausius -Mossotti formula, 26, 28, 35, 
110, 169 ft. 

Copper halide, 158. 

Debye equations, 70 f¥., 78 ff., 120 ff. 
Decay fimction, 6 ff. 

Diamond, 108, 109, 110. 

Dichlor methane, 116. 

Dichloral propane, 138. 

Dielectric constant, 1 ff. passim. 

(static), 2, 16 ff. 

general statistical theory, 

36 ff. 

of dilute solutions of polar 

molecules in nonpolar substances, 
30 ff., 118 ff. 

of gases, 28 ff., 116 ff. 

of ionic crystals, 149 ff. 

of mixtures, 47 ff. 

of polar liquids, 49 ff., 130 ff. 

of solids, 63 ff., 130 ff. 

spherical molecules, 31, 33ff. 

temperature -dependence of, 

12-13, 28 ff., 46, 48-61, 116-17, 
186 ff. 

(complex), 6, 62 ff. 

frequency dependence of, 8, 

62 ff., 73, 115 ff., 118 ff., 130 ff. 

temperature-dependence of, 

76 ff. 


Dielectric constant, relation between 
real and imaginary parts, 8, 162. 
Di-isopropyl ketdno, 134. 

Dilute solutions of polar molecules, 
in benzene, 120. 

in non-polar liquids, 

30 ff., 89, 118 ff. 

in paraffin, 120. 

in solids, 118 ff., 125 ff. 

Dipolar gases, 28 ff., 1J6 ff. 

— interaction, 21 ff., 33 ff., 38 ff., 
51 ff. 

— liquids, 49 ff., 83 ff., 130 ff. 

— solids, 53 ff., 79 ff., 130 ff. 

Dipole, 19 passim. 

— molecular, 26, 105 passim. 
Distribution of relaxation times, 91 ff. 

Electric dipole moment, 4, 15 ff. 
27 ff., 105 ff., 166 ^.passim. 

— displacement, 3. 

Energy, static electric field, 9 ff., 160. 

— loss in peidodic electric field, 5 ff., 
13, 161. 

Entropy, 9 ff. 

Ethelyne cyanide (solid), 133. 

Forces between dipoles, 21 ff., 33, 
36 ff., 48 ff. 

Free energy, 9 ff., 40. 

Frequency dependence of dielectric 
constant, see Dielectric constant 
(complex). 

Gases, 28 ff., 110 ff., 115 ff. 

Halides, 158. 

Harmonic oscillator, 16, 63. 

Helium, 111. 

Heptane, 123. 

Hindered rotation, 62. 

Hydrogen chloride, 105. 

— peroxide, 105. 

— sulphide (solid), 132. 

Hydroxyl radical, 107. 

Inert gases, 108, 111. 

Internal field (Lorentz), 22 ff., 163 ff. 
(Onsager), 26 ff., 163 ff. 

Ketones, 107, 126 ff., 142 ff. 
Kirkwood’s formula, 49 ff., 137. 
Krypton, 111. 

Longchain molecules, 113 ff., 126 ff. 
Loss angle, 14, 78 ff., 101, 121 ff. 
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Magnesium oxide, 168. 

Methane, 116, 117. 

Methyl benzoate, 123, 124. 

— chloride, 116. 

Models, for crystalline solids, 54 ff. 

— for Debye equations, 78 ff. 

— for complex dielectric constant, 
63 ff. 

— for static dielectric constant, 16 ff. 
Molecules, 105 ff. passim. 

Neon, 111. 

Nitrogen, 111. 

Non-polar molecules, 26, 105 ff. 

Onsager’s formula, 33 ff., 49 ff., 53, 
130, 171. 

Optical constants, 14, 162. 

— polarization, 105 ff., 149. 
Order-disorder transitions in dipolar 

crystals, 63 ff., 132, 146. 

Oxygen, 108, 111. 

Paraffin, 108, 113 ff., 125 ff. 

Pen ta -methyl chlor- benzene, 138. 
Pentane, 114. 

Phenolic resin, 135. 

Polar molecules, 26 ff., 105 ff. passim* 
Polarizabihty, 28, 105 ff. 

Polarization, 2, 64. 


Polarization waves, 149 ff. 

Power loss of dielectrics in periodic 
field, 13, 73 ff., 92 ff., 98 ff. 

Rate of unimolecular reactions, 81. 
Reaction field, 25, 31, 34, 41, 165. 
Refractive index, 28, 163. 

Relaxation time, 73 ff., 88, 91 ff., 
121 ff., 128 ff. 

Resonance absorption, 98 ff., 173 ff. 
Rubber-sulphur compounds, 135. 
Rubidium halide, 168. 

Self-energy, 144, 145, 165, 169. 
Sodium chloride, 109, 157. 

Strontium oxide, 168. 

Superposition principle, 6. 

Tallium chloride, 158. 

Temperature -dependence of dielectric 
constant, see Dielectric constant. 
Tertiary butyl chloride, 138. 
Titanium oxide, 109. 

Viscosity, 84 ff., 89, 123. 

Water, 107, 137 ff. 

Xenon, 111. 


PRINTED IN GREAT BRITAIN 
AT THE UNIVERSITY PRESS, OXFORD 
BY CHARLES BATEY, PRINTER TO THE UNIVERSITY 




date of issue 

This book must be returned within 3/7/14 
days of its issue. A fine of ONE ANNA per day 
will be charged if the book overdue. 




